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CHAPTER | 


Mathematical methods and 
differential equations 


1.1 Introduction 


In the first two chapters of this course we need to lay the groundwork for 
the more advanced material that comes later, so here the emphasis is on 
consolidation of material that you may have seen before. We also need to 
be sure that you are familiar with the notation that we shall use later. 


In this chapter we concentrate on areas that are essential components of an 


applied mathematician’s tool-kit: functions, complex numbers and di 
tial equations. In addition to reminding you of some useful techniques, we 
also provide you with a number of exercises that are preparatory to material 
that will be discussed in later chapt 


1.2 Some initial definitions and 
results 


1.2.1 Functions 


A function f is a mapping from a set A known as the domain, to a set B 
known as the codomain, in which every element of A maps to a unique 
element of B. The notation 


y=f(z) for ceA 


is used to encapsulate such a relationship, where z is an element of A, y is 
an element of B, and the symbol € is an abbreviation for ‘is a member of” 


Often A and B are sets of real numbers, in which case f is said to be a 
real function. For example, f(x) = a? for x , where R i 
numbers, is a real function. The set R is often referred to as the real line. 
A real function f is even if f(a) = f(—2x). and odd if f(x) = —f(—«). For 


example, y = x? is even, while y = 2° is odd. 


the set of real 


A function g is an inverse of the function f if f(g(x)) = x, and it is common 
practice to denote an inverse of a function f by f~'. For example. if 
f(x) = Vx for x > 0 (where /x denotes the positive square root of x), then 


f*@) 


because Vs 
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Later we shall meet functions for which the domain is a set of real numbers 
and the codomain is a set of complex numbers. Such functions are known 
as complex functions of a real variable. z 


We shall also meet functions for which the domain is itself a set of functions. 
Such functions are often known as operators or transforms. This may 
at first seem a rather strange concept, but you have in fact met the idea 
before. Whenever you differentiate a function, you are actually applying 
the differential operator to a function to obtain another function, namely 
its derivative. So differentiation is in fact an operator that maps functions 
to their derivatives. We shall not discuss operators in any great depth, but 
the concept will occur again when we introduce the Fourier transform later 
in the course. 


The domains and codomains of specific functions are often defined implicitly. 
For example, for the real function 


f(z) = 


it is implicitly assumed that the domain is the set —1 < x < 1 since this is 
the largest set of real numbers on which the definition is sensible. Strictly 
speaking, we should distinguish between a function f and a function value 
f(x), but generally we shall make no such distinction and, for example, refer 


vi-2x? 


Functions themselves are sometimes defined implicitly. For example, the 
equation (of a circle) x? + y? = 1 determines two functions, y = V1 — 22 
and y = —/1—<%, and we require more information in order to discover 
which of these functions is intended. In either case the domain would be the 
set -l<a<l. 


i 


to the function y = 


Smooth functions 


Often we shall need a function to be ‘sufficiently smooth’. The minimum 
requirement for this is that the function should be continuous at the point 
in question. In intuitive terms, a function is continuous if its graph has no 
sudden breaks, but we can be a little more precise. A function f is said 
to be continuous at a point « =a if jim f(x) = f(a). In other words, as x 


approaches a from above (a > a) or below (a < a), the function value f(x) 
approaches the value f(a). 


In order for the function to be ‘smooth’, we also require that it should be 


differentiable at the point in question. This means that for a function f It can be shown that a 
: . rene 2 — f(a) function that is 
to be smooth at the point 2 = a, we require that the limit aaa me differentiable at a point 
= 


—a + 
£ 4 CRaiteaeene * ene must be continuous there. 
exists. Again in intuitive terms, a smooth function is one for which the 


graph has a tangent at the point in question. 


The derivative of a function y = f(x) is often denoted by = dy , although we 


dx 
also denote the derived function by f’, so f’(x) = au 


dy 
The second derivative is usually written as f"(x) or —> Fra <Y while the higher 
d"y 
derivatives are written as f(z) or = for n = 3, Aye 


When we require a function to be ‘sufficiently smooth’, we generally mean 
that the function should be differentiable a sufficient number of times at the 
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relevant points of its domain. Sometimes this may imply that the function 
should be differentiable ‘an infinite number of times’. (For example, for the 
sine function f(x) = sina, we have f’(x) = cosz and f(r) = —sinx, thus 
the function is differentiable as many times as we please.) 


A polynomial is a function of the form 
p(w) = ap + ayx + aan? +--+ +anx", 
and the polynomial is of degree n if a, 4 0. Clearly a polynomial is ‘suffi- 
ciently smooth’ at any point on the real line because it can be differentiated 
as many times as we please. 
‘ * ‘ a p( a: 
A rational function is a function of the form a where p and q are 
Qe 


polynomials, and such a function is ‘sufficiently smooth’ at every point on 
the real line except the roots of the equation q(x) = 0. At those points the 
rational function is discontinuous. 


1.2.2 Series 


Convergence of series 


Example 1.1 
You may have already met the geometric series 

S,(a) =1+a+a?+---+a"}, 
where n is an integer and a is some specified number. Such a sum can be 
written in an alternative form if we multiply both sides of the above equation 
by a to obtain 

aS,(a) =a+ ae. a", 
then subtract one equation from the other and rearrange: 

er 


Sn(a) = =—. 


(1.1) 


For large values of n, and for any choice of the number a in the interval 
a" at 


-l<a<\1, the term Ta is small, therefore S,,(a) ~ . More pre- 


cisely, 


im Sn(a) = ' (1.2) 


b= 

provided that -l<a<1l. @ 

Tt is this final statement that leads us to the definition of the convergence 
of an infinite series. We begin by considering a finite sum 


Sp = Ag + Ar + Ag +++ Ana, 


n-1 
which may be conveniently written in the form S, = SS A; We then say 
Es k=0 
that the infinite series Yo Ae is convergent to the sum S if lim S, = S. 
ne 


k=0 
Ifthe limit lim S,, does not exist, then we say that the series is divergent. 
ne 
n=l 
A finite sum such as SS A,, is often known as a partial sum of the infinite 
series. k=0 


The symbol ~ is read as 
‘is approximately equal to’. 
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For the previous example we could then write if-l<a<1l, 


i 
Fr Fe, 
and S,(a)=1+a+a?+---+a"" would be the nth partial sum of the 
infinite series. 

Sometimes it is convenient to write an infinite series in a more informal 
manner, and for example we may write 


ox 
So Ak = Ao + Ar + Ag tert Ante 
k=0 

where the final set of dots indicates that the series continues indefinitely. 


Exercise 1.1 


Show that for -1 <a <1, R,(a) = }~ a* can be written in the form 


a” ken 


R,(a) = 


l-a 


Rate of convergence 


In some circumstances it is not sufficient simply to know that an infinite 
series converges; we may also need to know how rapidly the series con- 


if 


1 
verges. In the previous example we saw that the infinite sum was =r 


—1<a< 1, which means that S,,(a) — is small when n is large. It is 


12 
instructive to rearrange equation (1.2) into the form 


=1+a+a?+---+a"! + R,(a), 


n 


a 


where R,,(a) = is known as the remainder term. The rate at which 


the remainder fuera tare to zero determines the rate of convergence of 
the series, and, loosely speaking, in this case we can see that R,,(a) behaves 
like a”. However, intuitive statements like this can be made more precise 
by introducing some additional notation. 


The O notation 


The O notation is a convenient device for specifying an order of magnitude. 
For example, sin x is of the same order of magnitude as x when x is small. To 
make such a statement more precise, we define f(x) = O(g()) if there exists 
a constant M such that |f(2)/g(a)| < M for all sufficiently small values of x. 
Thus, for example, it is possible to show that sin = O(«) and 2x + 32? = 
O(x) if x is small. In the context of our previous example, we could write 
R,,(a) = O(a") for a sufficiently small, to indicate that |R,,(a)/a"| < M for 
some constant M > 1. So if we are concerned only with orders of magnitude, 
we could write 
re =l+ata?+---+a"? + O(a"). 

More generally, for a function f(x) = g(x) + an2" + @nsi2"*! +--+ with 
constant coefficients an, @n+41,---, we can write f(x) = g(x) + O(x") for suf- 
ficiently small values of «x. 
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Example 1.2 

Given f(x) = O(x?) and g(x) = O(«*) when z is small, establish the orders 
of magnitude of f(x)g(x) and f(a) + g(x). 

Solution 


We know that |f(x)| < kila|? and |g(x)| < ke\a|$, for some constants ky 
and ky. Therefore | f(s)g(:)| < kikalx|°, so f(x)g(x) = O(a). 


Also, we have Here we have used the 
3 triangle inequality 
[F(2) + 9(2)| < [9(e)1 + lg(x)] < halal? + halal? < 2h fa? ta Be tal 
if x is sufficiently small. Hence f(x) + 9(x) = O(z?). © (see Exercise 1.14). 


Taylor series 


It can be shown that a function f that is sufficiently smooth near a point 
x =a can be represented as a series Note that 
(k) ol=1 
f Wo. (1.3) and 
‘ L(a) = fla). 


f(x) = > ax(x—a)* where a, = 
k=0 


Such a series is known as the Taylor series (or Taylor expansion) for f 
near the point «=a, and it is generally convergent for x in an interval 
a—p<x<a+p (which can be written as |x — a| < p) for some positive 
constant p. The largest possible choice for p is known as the radius of 
convergence of the series. 


A sequence of numbers such as a),@2,a@3,... is often abbreviated to {a,}, 
and in this case these numbers are known as the coefficients of the series. 


Example 1.3 
Write down the Taylor series for the function f(«) = TS near x = 0. 
Solution 
- thi son, f(x) = EE so pt (0) = (—1) FA! 
For this function, f\“)(a) = Gea" So f\*)(0) = (—1)*k!. and the co- 


efficients of the Taylor series are a, = (—1)*. Thus 


AY ge 
ee wen 
ne XI 1k, 


Exercise 1.2 


Find the largest interval |x| < p on which the series 
ea 
1—a? +at— 29 +---= S>(—-1)ka* 
k=O, 


is convergent. 
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1.2.3 Some elementary functions and their 
derivatives 


You will probably have previously encountered most of the following func- 
tions and their Taylor series near x = 0. (We include the geometric series 
again for completeness.) We also give the interval of convergence for each 


series. 
1 
Tog Ee he for |2| <1, 
Cie a2 Bok 
expr =e" =1+ er o put a = +-=> > forr ER, 
k=0 
(—1)"22n+1 _ (eit 
Qn+D! eS: Qk+D! for z € R, 
(-1)"22" (—1)ka2* e 
Qn)! a > QR)! for x € R, 
k=0 
2 pk 
See for -l<a<l. 
k=1 


For this course, the natural logarithm of x will be denoted by Ina, although In is the logarithm to 
other authors write log, x or just log. For x > 0 the function In is defined _ base e. 

to be the inverse of the exponential function, so exp(In.x) = x for x > 0. The 

radius of convergence of the Taylor series for In(1 — x) near x = 0 is 1. but 

the series diverges at « = 1 and converges at x = —1. 


Many other series can be obtained by simple manipulations of the above 
series, as the following example shows. 
Example 1.4 


Write down the Taylor series for the following functions near t = 0. 


1 
(a) In(1 + 2¢) (b) a3 


Solution 


(a) Replacing 2 by —2t in the series for In(1 — x), we obtain 


2 (_oy)k 
in +2)=-S>! 2 : 
= 


(b) It can be shown that we are justified in differentiating the geometric 
series term by term to obtain 


aa = é (4) = tae Ss tnx" t+... 
~ Se 
k=l 
Then replacing x by 3t gives 


1 


20 
Gay 14 6t4 270? +---+n(3t)""? +---= Soke)! 


Exercise 1.3 


Show that sinxIn(1 + x) = 2? + O(z*). 
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The remaining trigonometric functions can be defined in terms of sine and 
cosine: 
sinr cost 1 


tang = » cor=— = + sec : ——« 
cosa sing tana cost sing 


The hyperbolic functions 


The hyperbolic sine and cosine functions, sinh and cosh, are defined by Sinh and cosh are 
pronounced ‘shine’ and 


~ 2n+1 
sinha = cee etree = Aras ‘cosh’, respectively. 
2 (2n+1)! 
convergent for x € R, 
e+e* a 
cosh a = =1lt a teeet fos. 
2 2! (2n)! 
» ek 
= + convergent for x € R. 
> (2k)! 2 


The remaining hyperbolic functions are defined by 


sinh cosh x 1 
tanha = he’ cotha = = hea tan Tanh, coth, sech and 
es bee a tanhz cosech are pronounced 
sech z = , cosecha = — A ‘tansh ‘coth’, sesh and 
cosh sinhz ‘co-sesh’, respectively. 


The graphs of sinh, cosh and tanh are shown in Figure 1.1. 


Figure 1.1 ‘The graphs of sinh, cosh and tanh 


Exercise 1.4 


From the definitions of sinh and cosh, show that cosh? x — sinh? x = 1. 


Example 1.5 


Show that a (cosh! x) = 
dx 


Solution 
If w=cosh~! x, then coshw =x. Differentiating with respect to x, using 
the chain rule, we have sinh w =~ =1. Thus 

dw 1 ul 


1 
= = = tf 
dx sinhw /eosh?w—1 va2—-1 
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The Gaussian function 
The Gaussian function 
y = exp[—(a — )?/207), 


where j. and o are constants, arises in statistics. Its graph is shown in 
Figure 1.2 for various values of jz and o. 


Figure 1.2 Graphs of the Gaussian function 


Identities and derivatives 


Various identities for the trigonometric and hyperbolic functions arise quite 
often, and you may find the following lists useful. 


cos? 2 + sin? 2 =1 cosh? x — sinh? ¢ = 1 
sin(x + y) =sinxcosy+cosxsiny | | sinh(« + y) =sinhxcoshy + cosh«sinh y 
cos(2 + y) =cosxcosy —sinasiny | | cosh(x + y) =coshxcoshy +sinhasinhy 
tana + tan tanh x + tanh 
tan(e +4) = thing) 
1—tanztany 1+ tanhztanhy 

Function f(a) | Derivative f’(x) Function f(a) | Derivative f’(x) 
a” (n#0) | nz" 
Ee e Ina 1/x 
sing cos © sinh «x cosh a 
cos. —sing cosh x sinh x 
tana sec? x tanhx sech? x 
cota —cosec? x coth a —cosech? x 
seca sec x tan rc sech x —sech x tanh 
cosec x —cosee x cot t cosech x —cosech x coth x 
sing! & a sinh”! x — 

V1-2? : Va? +1 
cosa se cosh7! x Z 

ie Ve=1 

tana! x 4 tanh"! x : 

1+2? 1-2? 


1.2.4 Integration 


A primitive of a given function f is a function F such that F’= f. The 
term ‘indefinite integral’ is often used as an alternative to ‘primitive’, and 
the notation i {(x) dz is used to indicate an arbitrary indefinite integral 


of f. Thus indefinite integration arises as the ‘inverse of differentiation’. 


1.2 Some initial definitions and results 


On the other hand. definite integration originates from the need to eval- 
uate a limit of a sum, for example arising from the area under a curve. 


These two notions of integration are linked by the fundamental theorem of 
calculus (see below), but it is important to appreciate that they are quite 
distinct concepts. 


In an integral of the form / f(x) dx, the function f(x) being integrated is 


called the integrand, and the variable z is called the integration variable. 
Another common way of expressing such integrals is to write the integrand 


last, i.e. fu f(x). Both notations are used in this course. 


Definite integration is based on the concept of a limit of a sum. Suppose 
that we are given a function f(a) defined on the interval a < x < b. We first 
divide the interval into N subintervals, with end points 


a=2% <a <9 < +++ <ay_-1 < tn =b. 
N 
Then we form the sum S = > f(x,)(- — 2-1). Often it is convenient to 
r=1 
choose subintervals of equal length, which we might denote as dx, and then 
N 


this sum becomes S = y f(a,) 6x, which is often abbreviated to > f(x) dx. 
r=1 

Finally, we take the limit as the number of intervals increases, and the 

maximum length of the subintervals shrinks to zero, giving 


b N 
[[ te)de= jim Y fener 27-2) 


It is important to appreciate that this is how the integral is defined, but to 
evaluate an integral we generally try to employ the fundamental theorem 
of calculus. This tells us that if we are able to find a primitive F() such 
that F'(x) = f(a), then 


b 
Ht f(a) dx = F(b) — F(a). 


Example 1.6 


Find the area bounded by the curve y = Inz, the z-axis, and the lines « = 2 
and x = 3 (see Figure 1.3). 


Solution 


Figure 1.3 


If we approximate the required area A by rectangles, as shown in the figure, 
then A ~ }7In(«) dx. Now we take the limit as the number of intervals 
increases, and the length of the subintervals shrinks to zero, to obtain 


A= [eae 


In an elementary 
introduction to calculus, 
such a sum S is usually 
identified with an estimate 
for the area under the 
graph of f(r) from «=a 
tor=b. 
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To evaluate this integral we may use integration by parts, which gives 


3 F 3 
A=] In(z)dr= {rine dx Try differentiating «In x if 
2 
2 2 you cannot see why this 
= [ena - a3 works. 


=3In3—2In2—1=In(27/4)-1. 


Example 1.7 


Suppose that we are given a wire that does not vary in any way along its 
length: i.e. its cross-section stays the same, and the material from which it 
is made remains fixed. We say that such a wire is uniform, and we expect 
the mass of any segment of such a wire to be proportional to the length of 
the segment: a piece twice as long would weigh twice as much. 


The ratio of the total mass of the wire to its total length is known as its 
linear density, and if we denote this quantity by p, then the mass M of 
a length L of wire is pL. If M is measured in kilograms (kg) and L in 
metres (m), then the units of p are kilograms per metre (kgm™'). 


Now suppose that the wire is non-uniform: how should we extend our def- 
inition of linear density? Imagine that the wire is placed along the «x-axis, 
and concentrate attention on a segment of wire of length dx containing a 


particular point . If the mass of ee segment is dm, then we define the 
. m 
linear density p(x) at « to be lim =. 
br—0 6. 


In order to be more specific, let us consider a wire which is 1 metre long, for 
which p(a) = 2a +1, so that the density increases steadily from the value 
1kgm7! at the origin to 3kgm7! at the other end of the wire. The total 
mass of the wire can then be expressed as a definite integral. 


First we divide the wire into N subintervals, with end points 

O= 29 < t < 29 <-** << oy < ty = 1. 
If the intervals are small, then p(x,) is a good approximation to the linear 
density at all points of the interval z,_) <x < ry, so the mass of this short 
segment of wire is 6m, ~ p(x,)(x, —«,—;). Adding all of these small masses, 


we obtain an approximation for the total mass of the wire: 
If we choose equal 


n n 

f= aA e subintervals, then this 
M= 2 ie = » P(r) (tr — Fr—1) approximation can be 
r= r= 


s written as M ~ S7 p(x) dx. 
= ¥ ala) éx,, where dap = Ly — Ly—1- 
r=1 
This approximation becomes more accurate as we increase the number of 
intervals and the length of each subinterval shrinks to zero. Thus in the 
limit we have 


1 
m= [ p(x) da. 


To evaluate this integral for the particular function under consideration, we 
employ the fundamental theorem of calculus: 


1 1 
= as = 
M i pala [ er+iar 


ie. the wire has a mass of 2kg. 


x? + 2]) = 2, 


Example 1.8 


Calculate the volume of a sphere of radius r. 
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Solution 

First divide the sphere into a number of vertical slices of thickness dx (rather 
like slicing an onion), then approximate a slice at position « by a cylinder 
of depth dx and radius y = Vr? — x? (see the shaded region in Figure 1.4). 


ox 


Figure 1.4 A sphere divided into elemental ‘cylindrical slices’ 


The volume of such a cylinder is wy?éa = a(r? — x?)dx, so the sum 

S=DompPor = 7D (r? — 2? )6x 
provides an approximation to the required volume. In the limit as the num- 
ber of subdivisions increases and dx: tends to zero, the sum becomes a definite 

_ 
integral and we obtain the required volume, V =a | (r? —«*) dx. Now we 
may use the fundamental theorem of calculus to opin 
r 

Yor (r? — a?) dx = m [r?x — 4 iE 

-r 


Integrals over an infinite range 


The integrals discussed so far have been defined over finite intervals, but we 
shall also need integrals defined over infinite ranges. 


Example 1.9 


Suppose that we are told that an object is at the origin at time t = 0, and 
that it moves along the z-axis in such a way that its velocity at time ¢ is 
given by v = cost. What is the maximum displacement of the object from 
the origin? 


Solution 

d. 
We know that = =v =cost. so x(t) is a primitive (or indefinite integral) 
of v. Thus x(t) = | costdt =sint +c, where cis a constant. We choose c so 


that x(t) satisfies the initial condition «(0) = 0, thus ¢ = 0. So x(t) = sin#, 
and it is clear that the maximum value of x(t) is 1. 


Note that this result can also be found using definite integration, by setting 
the lower limit equal to the initial condition x(0) = 0. At any given time T, 


‘a 
we have 2(T) = [ costdt=sinT. @ 
0 
In the previous example the object moved backwards and forwards about 


the origin (consider the graph of x(t) = sin t), but what happens if the object 
moves always away from the origin but at a slower and slower pace? 
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Example 1.10 


Suppose that we are told that an object is at the origin at time t = 0, and 
that it moves along the x-axis in such a way that its velocity at time ¢ is 
given by v = 1/(1+1)?. What is the maximum displacement of the object 
from the origin? 


Solution 


As in the previous example, we can use either definite or indefinite inte- 
gration. We choose definite integration here, but you might like to repeat 
the calculation using indefinite integration, using the initial condition to 
calculate the value of the arbitrary constant. 


At any given time T, we have 
T T 
1 1 1 
= dt = |— =1- i 
aKa) Jo (L+t)? re 14+T 


As time T increases, the displacement 2x(T’) increases from zero, but we can 
see from Figure 1.5 that it never actually achieves a maximum value. Note 
also that the slope of the graph, v(T), decreases as T increases. 


% 


x(T) = 1-W+T) 


Figure 1.5 The graph of x(T) = 1—1/(1+T) 


As T increases, the value of a(T) = 1 — approaches the value 1 from 


below (but more and more slowly as the speed of the object decreases). So 


ie reek 
im ( f aa) = 1, and we normally write [ tan i | 


Thus we are led to our general definition of an integral over an infinite range. 
r 
If the limit jim (/ sat) exists and equals a finite value L, then we 
—-On a 00 
say that the infinite integral is convergent and write [ f(t)dt =L. 
a 


If the limit does not exist, then we say that the integral is divergent. 
oa 


In the context of the previous examples, we would not expect | u(t) dt to 


a 

exist if the object fails to slow down, but you might reasonably expect that 

the condition jim u(t) = 0 would be sufficient to ensure the convergence of 
30 


28 
[ u(t) dt. However, this is not the case, as the following example shows. 
Ja 


Example 1.11 


Suppose that we are told that an object is at the origin at time t = 0, and 
that it moves along the z-axis in such a way that its velocity at time ¢ is 
given by v =1/(1+1). What is the maximum displacement of the object 
from the origin? 


1.2 Some initial definitions and results 


Solution 


As in the previous example, the object moves more and more slowly as it 
moves away from the origin. However, we now have 


T 
1 
2(T) -f pap = +7). 


and 2(T’) increases without bound as T increases. So here we have a cir- 
cumstance where the object is moving ever more slowly, but eventually it 
will move beyond any fixed point that we care to mention. In this case 


2S 

1 

i Tot dt does not exist, and the integral is divergent. 
0 


Sometimes it is possible to establish that an integral is convergent without 
actually calculating the value of the limit. Suppose that v(t) > 0 and that 


T 
2(T) = i v(t) dt <K 
for some fixed constant K and all T > a. In this case, differentiating with 
respect to T gives x'(T) = v(T) > 0, therefore x(T) cannot decrease with 


time T, but we also know that it cannot exceed the value K. So it must 
approach some fixed value (which is less than or equal to K). 


Example 1.12 


dt is convergent. 


erie 
Show that the integral | 
1 


Solution 
We use the notation introduced in the previous examples. In this case, 
t 


u(t) = — > 0ift>1, so 


T t 

2(T) = [ at 
1 

et 


increases with T. But — <e~' ift > 1, so 


t 


T. T T 
ar)= [ : ars [ edt = [-e]) =e -e7 <et. 


So 2(T) is non-decreasing and bounded above, so we can conclude that 
oo tt 
-. dt is convergent. 


jim (T) exists and therefore the integral 
100 i 


The Gaussian integral 


The following integral identity is very useful in many contexts where Gaus- 
sian functions appear: 


0 
f exp(—az*) dx = yi where a > 0. (1.4) 
aan a 


This important relation is obtained by a clever trick, where the square of 
the integral is evaluated using polar coordinates. It will be derived in the 
next chapter in the section on multiple integrals. 
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Exercise 1.5 


(a) The following result is often useful when estimating the size of functions that 
can be expressed as integrals. 


If f(x) is a function such that |f(x)| <M for a <x <b, 


I f(x) dx 


Use this result to show that Ina < a for a > 1. 


then < M(b—a). 


(b) Using the method of integration by parts, or othe: evaluate the following. 


' Ting ve T ing 
of = i) f sa de 


(c) Show the following. 


Ing 

(i) [ a ae eyes 
1 z 
Ina 

(ii) [ AF dr is convergent. 
1 e 


Exercise 1.6 
The following result is often useful when comparing the size of integrals. 


If f(x) and g(a) are two functions such that f(x) < g(x) for a <a < b, 


b 
then [ f(a) dx < [ g(x) dx. 


20 
Use this result to show that [ exp(—2x") der is convergent. 
0 


[im First split the integral as 


oS 1 oS 
[ exp(—2) de = [ exp(—e2) de+ : on(-2*) ae] 
A 0 1 


1.2.5 Partial fractions 


The following process of adding fractions is probably familiar to you: 
2 3 _ A%x+2)+3(~—-3)_ —_—(w—-1) 
a-3 2+2 (x — 3)(x + 2) (x — 3)(a +2)" 
The reverse process is known as expressing a rational function in terms of 


partial fractions, and the procedure is useful in differentiation and integra- 
tion, and when expanding a rational function as a power series. 


Suppose that we are told that the identity 
5(a@—-1) _ A B 

(x — 3)(a + 2) 

holds for some constants A and B, and we want to determine these constants. 


eB. (1.5) The symbol = is often used 
3° r+2 in such a case to emphasise 
the fact that the equality 
spplies far all values ofr, 
There are various methods that can be used to find the values of A and B; but we shall normally 


perhaps the simplest is to substitute various well-chosen values of «. use = unless we wish to 
clarify a point. 
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For example, putting « = 0 in identity (1.5) gives the equation 


eee 

cs 3° 2 
while putting x = 1 gives 

A,B 

j= + =. 

vas ia 
These two equations for A and B are sufficient for us to establish that A = 2 
and B= 3. 


Alternatively, we could multiply both sides of identity (1.5) by (a —3)(x +2) 
to obtain 


5(a@ — 1) = A(x + 2) + B(x — 3). 
Then equating the constant terms and the coefficients of « on both sides of 
this expression gives 

—-5=2A-3B and 5=A+B. 


Again we have sufficient information to establish that A = 2 and B= 3. 
The latter technique is often referred to as equating coefficients of powers 


of a. 
Thi . .. Pia) F 
This method works for a rational function Q) where the polynomial P(x) 


is of lower degree than the polynomial Q(x). 


Example 1.13 


1+ 6x — 2? 


ee ent deheeaabiana: 
Qr—-D(e@+2(a+1) in partial fractions. 


Express 


Solution 


Since the numerator is of lower degree than the denominator, we attempt 
to express the fraction in the form 
1+ 6x — 2? aA ee. Cc 
(Qa—1)(e@+2)(@+1) 22-1 r+2 «+1 
for some coefficients A, B and C. Then we have 


1462 —2? = A(w +2)( +1) + B(2x - 1)(a +1) + C(2x — 1)(a +2). 


Now we either set 2 = 4, 2 = —2, 2 = —1 in turn, or equate coefficients of 
powers of a, to obtain A = 1, B = —3, C = 2. Thus 
1+ 6x — 2? 1 3 2 


(Qc—I)(@+2)@+1) 22-1 24221 


Exercise 1.7 


Express the following in partial fractions. 
lla +5 


() Gr=Gerd 
(b) 2(a + 1)(2a +7) 
(a —1)(a + 2)(a + 3) 
342, x42 P 
(ec) — Hint: Try writing of =2r+ aa 
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Repeated and quadratic factors 


A minor adaptation of the method is required when a factor in the denomi- 


nator is repeated. For example, to express in partial fractions, we 


= 1 
Tae 
write ( 2) 
fet A b B 
(l-2)? 1-2 (1-2)? 
Multiplying by (1 —.2)?, we obtain 22 — 1 = A(1—.r) + B, and then equating 
coefficients of powers of x gives A= —2 and B= 1. 

A similar approach will deal with rational functions in which the denom- 
inator has a quadratic factor that will not factorise into real factors. For 
vtrtl 
(a? + 1)(a — 2) 

v+ct+l _ A Br+C 

(z2+1)(@—-2)  2-2° a? 41° 
The value of A can be determined by multiplying both sides by 2 — 2 and 
then putting x = 2. giving A = Z. Then we multiply by (a? + 1)(@ — 2) and 


example. to express in partial fractions, we write 


equate coefficients, giving B = -2 and C = # Thus 
vt+a+l ae ih plate 
(@?4+1)(e@-2) 5 \r-2° 2241) 


Exercise 1.8 


Express the following in partial fractions. 
r 


©) Booa46 
2 ints Writes == ea HO) 
(b) j-1F Hint: Write G1" 1+ Oa) 
(«=1)? 
©) G377e-D 
(d) cts 


Q+2)0 +2)? 


Example 1.14 illustrates some of the applications of partial fractions, but 
first we remind you of two useful facts. 


if 


=In(x+a)+e forr+a>0. 


S=l+atartetat+--- for |2| <1. 


Example 1.14 


Use the expression obtained in Example 1.13 to do the following. 


1+6z—2 
i dir fe 3. 
(a) Evaluate the inetd | ee lx for x > 3 
1+ 62 —2? 


(b) Express f(x) = G as a Taylor series near x = 0, 


2x —1)(x + 2)(x +1) 
af 


(c) Calculate as 


All quadraties will factorise 
if we allow complex 
numbers in the factors, In 
fact, the technique 
discussed above can readily 
be extended to this case. 


1.2 Some initial definitions and results 


Solution 


(a) Using the fact that f de = Ine for 2 >0, we see that 
foams 1+62 — «2? i 
(2a —1)(a + 2)(x +1) 


2 
-fea* -fz sgee+ {sae 


$In(2x — 1) — 3ln(x + 2) + 2In(x + 1) +e, 


since ¢ > $ 


(b) First we have 
1+ 62-2? eS ee eee 
(Qe-D(@+2)(e+1)° 1-2 1+}e lta 
then we recall that the geometric series can be written as 
1 
aes L+atarteta%+--- for fal <1. 
Replacing xr in this series first by -32, then by 2x, and finally by —2, 
we obtain 
1+ 62-2? 3 
a a eS Se tee a 
Qz— De +2)(@ +1) 3( 22+ (3 
+2(1—2+2?-254.--) 
—(1+ 2a + (2x)? + (2x) +--+) 


for |x| < 


rey 


fa : 
(c) First we notice that if g(a) = (1+azx)7!, then 3 = -3!a¥(1+ar)-4. 


So, choosing appropriate values for a, we obtain 


(sr | (eae ee ON: Sa 
dg \ (1-228 24 442)? (i +2)8 


a ee 
(1—22)4 “ (2+2)4 (1+2)* 


1.2.6 Complex numbers 


A brief history 
Historically, complex numbers arose from attempts to splve equations, such 
as x2 = —1, which have no real solutions. Early mathematicians would 


have said that such equations have no solutions, but in the mid-sixteenth 
century, Cardano and his contemporaries were performing calculations in- 
volving /—1. They discovered that assuming the existence of this so-called 
‘imaginary number’ enabled them to simplify the algebra involved in estab- 
lishing results involving ‘normal’ numbers. Euler introduced the notation 
i = ¥—1 in 1777. and established the relation e'" = —1 (which involves four 
of the most interesting numbers in mathematics). In 1799 Gauss proved his 
fundamental theorem of algebra, which states that every polynomial equa- 
tion has real or complex roots, and early in the nineteenth century Cauchy 
and others extended the notions of calculus to complex functions. However, 
these are generally matters that are better discussed in a course on pure 
mathematics. Here we are mainly concerned with the applications of com- 
plex numbers; we shall meet them again when we discuss Fourier series and 
transforms. 
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Complex algebra 


For our purposes it is adequate to introduce the imaginary number i with In a more formal 
the property that i? = —1, and to regard the set of complex numbers C as definition, of the kind that 
the set of all ‘numbers’ x + iy. where x and y are real, with the normal rules" May find in a course in 


ie pure mathematics, the set 
of algebra. Thus, for example, Gimay be defined 'to be the 


(1 + 2i)8 = (1 + 2i)(1 + 2%)(1 + 2%) set of ordered pairs of real 
: ai 3 numbers (x,y), with 
= 1+ 3(2i) + 3(2é)" + (2i) certain specified rules of 
=1+6i-—12-8i addition and 


multiplication. The 
complex number i is then 
identified with the ordered 
pair (0,1). 


=-ll-2. 


Real and imaginary parts 


It is common practice to write z = x + iy and then refer to the complex 
number 2, with x and y as its real and imaginary parts, respectively. We 
use the form z = 2 + iy so frequently for an arbitrary complex number that 
it is not usually necessary to add that we are assuming that x and y are 
real. 


An expression such as «+ iy is often known as the Cartesian form of 
a complex number 2 (by analogy with the Cartesian coordinates of the 
plane). It is also useful to introduce the notation Re(z) = x and Im(z) = y 
for the real and imaginary parts of the complex number z. You should note 
particularly that the imaginary part of z = x + iy is y. not iy. You should 
also be aware that two complex numbers can be equal only if both their 
real and imaginary parts are equal. Thus, for example, if 2 + iy = 2+7 
then we know that 2 =2 and y = V3. This process is sometimes known as 
equating real and imaginary parts. 


We also identify the set of real numbers with the set of complex numbers 
having zero imaginary parts. In other words, we regard the real number 3 
as being identical to the complex number 3 + 0i. So a real number may be 
regarded as a special case of a complex number. 


The Argand diagram 


Real numbers are often represented as points on a line, whereas complex 
numbers, having two independent components, require two dimensions. It is 
for this reason that they are often represented as points in a two-dimensional 
Cartesian system, and this representation of the complex plane is usually 
known as the Argand diagram (see Figure 1.6). In this representation the 
w-axis is often known as the real line and represents the set of real numbers. 
We may also refer to the «- and y-axes as the real and imaginary axes, 
respectively. 


If z=x-+ iy, then the conjugate of z is defined to be z* = + —iy (see Some authors use = rather 
Figure 1.6). than 2°. 


The modulus of z = x + iy is defined to be |z| = \/x? + y?, and the identity 
|z|? = 2 2* is often useful (as you will see in the exercises). 


It is easy to establish the results 


for any complex number z. Also, 
\z122| = l2i||z2], |21/z2] =|2il/lza]_ and (2122)" = 2723 


for arbitrary complex numbers 2; and 29. 


1.2 Some initial definitions and results 


ay 


Figure 1.6 The Argand diagram 


Example 1.15 


243% . Fe 5 
——.. find 2" in Cartesian form. 


Given z = 
1—i 


Solution 


We can convert z into Cartesian form by writing 
243i 1+i 
+ 30 +1 1 (Gi-1); 


RE pS 
1-i 1+i ? 
then we can see that 


2" = -$(1+ 5i). 
Alternatively, we can replace every i by —i in the expression for 2, to obtain 
-_ 2-3 
~ 1+i 


and then 


= is x = —3(1 + 5i) 
It would be equally acceptable to write 
are 1 5i 
— (Berg) or <3 aa a 


Exercise 1.9 


‘ 
= £20 in Cartesian form (where a,b, cand dare real, with 2 +? #0). 
me 


Express = 


Exercise 1.10 


Rewrite the complex number z = (1 + 37)? — (3 — i)? in Cartesian form a + ib, where 
a and b are real. 


Find the modulus and the imaginary part of z. Also find 2*, and calculate 


23 
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Polar form and de Moivre’s theorem 
Using polar coordinates, z = x + iy becomes 
z=r(cos@+isin@), Some authors abbreviate 


cos@ + isin@ to cisé. 
known as the polar form of z, where 
y 


r=2?+y?=|2|20, cosé=————., ae 


On an Argand diagram, the modulus r is the distance of the complex 
number from the origin, and the angle @ is known as its argument (see 
Figure 1.7). 


sin@ = 


Tt can be shown that 
2" =r"(cosné + isinné) 


for any positive or negative integer n, a result which is often known as 


de Moivre’s theorem. Abraham de Moivre (or 
a s Demoivre) was an 
Note that it is possible to define powers of complex numbers 2” for any real eee Halla aud 


number n. In fact, de Moivre’s theorem also holds in this case, as will be a contemporary of Sir 


demonstrated in Exercise 1.13. Isaac Newton, He adopted 
various spellings of his 
name while he lived in 


Exercise |.11 England. 


Use de Moivre’s theorem and the fact that cos? @ + sin? @ = 1 to show that 


cos 50 = cos 0(16 cost @ — 20. cos? @ + 5). 


The argument of a complex number 


The polar form of a complex number is not unique, because its argument 
is not unique. For example, the complex number 1 +7 can be written in 
polar form as V2 (cos } + isin 7), but it would be equally valid to write 
V2 (cos % + isin %), or indeed 2 (cos(% + 2km) + isin(| + 2k7)) for any 
integer k. 


“Y 


Figure 1.7 Two possible values of the argument of a complex number 


An argument of an arbitrary complex number z = x + iy is defined to be an 

angle 0 for which sin @ = y/\/x? + y? and cos@ = x/\/x? + y?. 

If z = r(cos@+ isin@) and —7 < @ <7, then @ is known as the principal Some authors define the 
value of the argument of z. An argument of z is written as arg(z), while principal value to lie in the 
the principal value of the argument is written as Arg(z). interval 0 < @ < 2r. 
Example 1.16 


Express the complex number z = V3 +i in polar form. 


1.2 Some initial definitions and results 


Solution 


If @ represents an argument of 2, then cos@ = V3/2 and sin @ = 1/2, and we 
see that 7 is a suitable value for @. Also, |z| = 2, so 

V3 +i=2(cos%+isinZ). 
Any value of the argument of the form § + 2k, where & is an integer, is 
equally acceptable, but 7 is the principal value of the argument. (When 


determining the argument of a complex number, it is often helpful to plot 
the point on an Argand diagram.) 


We can use de Moivre’s theorem to find roots of complex numbers. For 
example, to find the fourth roots of z = /3+i, first express it in polar 
form, 

= =2(cos(Z + 2k) + isin( + 27k)) , 
using the general value of the argument, i.e. with k as any integer. Then 
use de Moivre’s theorem to evaluate the fourth roots: 

24/4 = 2'/4 (cos (i (F + 2nk)) + isin (7 (F + 2th))) . 
For & = 0.1, 2,3, this gives four distinct complex numbers: 

29 = 2"/"(cos H+ isin), 21 = 2'/*(cos BF + isin YF), 

22 = 2'/*(cos 8 + isin Bt), 23 = 2'/4(cos 4 + isin YF), 


whieh are the four fourth roots of 2. Other values of k simply duplicate these 
roots — which can, of course, be expressed in Cartesian form if required. 


Exercise 1.12 


fa) Express the complex number z = 2 — 2i in polar form, 


(b) Find the cube roots of z = 2 — 2i in Cartesian form, 


The complex exponential function 


It is possible to extend the notion of Taylor series to complex numbers, and 
to define 
f i9 (6)? ig)” 
exp(is) =e =14 4, +: {i6) 


2! n! 
7 eo 8. oy 
a Wamu et a ica al al 


= cos@ + isin. 


This result is known as Euler’s formula. 


Any complex number in the form 2 = r(cos@+ isin@) can then be written 


in the equivalent: form z = re’. 


We can also use Euler’s formula to show that 
el 4 @-i0 cid _ --i0 


6 = —— in @ = 
cos 2 and sin i 


(which reinforces the close relationship between the trigonometric and hy- 
perbolic functions). 


Exercise 1.13 


Use Euler's formula to prove de Moivre’s theorem. 


25 


You should know that a 
complex number has at 
most n distinct nth roots. 


26 Chapter | Mathematical methods and differential equations 


Example 1.17 
Express the following numbers in polar form. 
1+iv3 
@int wae Gee 
l+i 
Solution 
(a) 1+%= V2 (cos F + isin ¥) In each case a multiple of 
2m could be added to the 
(b) 1+%V3 = 2 (cos § + isin ¥) argument. 
(c) Since (using the results of parts (a) and (b)) 1+i/3 = 2e'*/3 and 
1+i= V2e'"/4, it follows that 
1+iv3 _ 24¢i7/8 i Ae 
+i Ren v2ei"/!2 = V2 (cos +ising). # 
Complex numbers are often useful in calculations that would otherwise be 
rather tedious. Suppose, for example, that we wish to evaluate 
[et sion at, 
where a and bare real. This integral can certainly be found using the method 
of integration by parts (twice), but the following approach is probably easier. 
Consider the integral (}- © « 
(A) ee 
\ a= ss sn 
= f= dt ; Bteabh “+s 2 = enddd p esalet) 
=—-e2 +2 
it Ark 7 a 
= —— lorie AE cos (LEW SACLE 
ape + atn’* (25 (GAs Cb8)) 
Is isin(bt ; 
= at cos(bt) + Aging ) G-% 
a+ib 
— pat cos(bt) + dein (be) yee ib Ee 
a+ib a—ib 
at 
= a ((acos(bt) + bsin(bt)) + é (asin(bt) — bcos(bt))) + C. 


We can also write 
‘ q . toe 
4) 0S if ettelt dt = f e' (cos(bt) + isin(bt)) dt, = ftSeccoe) i (a aaCtet) 


so, equating real and imaginary parts of these alternative representations 


of J, we obtain , 


eat ; 
feo cos(bt) dt = 2aR (acos(bt) + bsin(bt)) + Cy Maik ha vou 
and fram clard 
oe ett . _ \€ 
fe sin(bt) dt = Sp (asin(bt) — eos(bt)) + C2 | La keg 


(where C = C, + iC). 
So we have found not only the required integral, but also f e™ cos(bt) dt 
with no extra effort. 


In a later chapter we shall meet complex numbers again in the context of 
Fourier series, and the following observations will prove useful. 


1.2 Some initial definitions and results 


For any non-zero integer k, 


[t= 2 [yt - 


sin(kr)=0, (1.6) 


ca 
ne 
while for k = 0 we have 
i ed9= | d0=2n. (1.7) 
=" —« 


To give you some impression of how these results will be used later, consider 
the function 


f(0) = Aye" + Ane”? + Age? + Aye? 


With the aid of the stated results, it is very easy to express the coefficients 
in this expression in terms of the function f. For example, to determine A3, 
multiply by e~* and integrate from —z to 7, to obtain 


a : a a 
f 1(0)e-* ao = Ay f a0 + Ae f edo 
x -n 


+s [ ao A [ 
ie 


. 
phichiciven tty = a af f(0)e- do. In general, 
us 7 


1g nid 
Ape 0 
An = 5 fi Sema 


for n = 1,2,3,4. 


End-of-section Exercises 


Exercise 1.14 

Prove the triangle inequality, that |a + b| < |a| + |b| for any two real numbers a 
and b. 

Exercise 1.15 


Show that 

d 
de 
Hence show that ¢ = Asinhwa + B cosh wa, where A and B are arbitrary constants. 


(sinh a) =coshz and (cosh x) = sinh x. 
dx 


+ . s 5 Pps 2 
is a solution of the differential equation ov? 
de’ 


Exercise 1.16 


. 
in partial fractions, and hence show that > 


we 
n(n +1) 


Express 


= 
n(n +1) 
Exercise 1.17 
Find the Taylor series for sinh(2«) near «= 0. 
Exercise 1.18 


Show that = 


=a+ S29 +O(2'). 


Exercise 1.19 


Show that f(x) = = is odd. Hence evaluate / ree dr. 


i+z a 


27 
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Exercise 1.20 


1 
=F and sketch the graph of y = sinh! x. 
+x 


Show that 2+ ¥1+2?>0ifxeR. 


Sew thay (entict a 
cE 


Calculate S (ne + V1+2)), and hence show that sinh™' x = In(x + V1 +22). 
di 


Exercise 1.21 


4 2 
4a* — 9: 4 
in partial fractions, and hence evaluate f° ==>" 5 ae. 


9r +4 


4, 
xpress ae-D@—) 


Exercise 1.22 


2 
' : 1 4 
For what values of the constant a is the integral [ - dt convergent? 


Exercise 1.23 This exercise is optional, 
20 oo 
Show that if i |f(t)| dt is convergent, then [ f(t) dt is also convergent. Applied mathematics is 
a to dependent on many results 


i a from pure mathematics. 
[Ean Comers Ce) LAC] This exercise is intended to 


give you some indication of 
the steps that would be 
required to justify some of 
our assumed results. 


: sint 
Show that [| “> d¢ is convergent. 


Exercise 1.24 


Express the following in Cartesian form. 


143i 2-3: 1 1 
t ‘ 
ies “Tam faa TH 


Exercise 1.25 


Express the following in polar form. 
1 2 mis 
(a) —— () (1+iv3)? (©) (1+ iv)" 


Exercise 1.26 


Show that cosh(i#) = cos@ and sinh(i#) = ésin@. 


Exercise 1.27 


Show that for any two complex mumbers 2 and w, 


jz+w)? + [2 — wl? = 2}z|? + 2)w)?. 


1.3 Differential equations 


A brief history 


The study of differential equations has a venerable history, and the subject is 
still central to modern mathematics because many physical systems can be 
described by differential equations. As early as 1676, Newton solved a differ- 
ential equation by means of an infinite series, although his results were not 
published until 1693, the year in which a differential equation appeared for 
the first time in the work of Leibniz. Shortly afterwards, Johann Bernoulli 
developed the method of ‘separating the variables’, and he also showed how 
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a differential equation of the form dy/dx = f(y/x) could be reduced to one 
in which the variables were separable. He and his brother Jakob succeeded 
in reducing a large number of differential equations to soluble forms. 


Developments in the eighteenth and nineteenth centuries were rapid and 
involved some of the greatest names in mathematics: Euler, d’Alembert, 
Lagrange, Laplace, Cauchy, Jacobi, Ampére, Frobenius and Lie, to mention 
but a few. 


From these passing references to its history, it should be evident that the 
study of differential equations is a vast subject, and here we can do little 
more than scratch the surface. 


Introduction 


While it is important to be able to solve various differential equations, it is of 
equal importance to be able to recognise that a given equation is of a certain 
type, for then you can refer to the method of solution in one of the many 
standard texts. You may have seen some of these differential equations 
before in other courses, in which case you can regard the corresponding 
exercises as revision, but they are included here for completeness. 


While the variables, constants and parameters in this section are generally 
real, much of this material applies also when they are complex. However, 
certain complications arise in integration of complex functions, so for the 
moment you may assume that everything is real. 


1.3.1 Initial classification of differential equations 


Essentially, the study of differential equations is a matter of assigning them 
to various categories and discovering as much as possible about each such 
category. Here we shall make a start by dividing them into two major groups: 
ordinary differential equations and partial differential equations. 


Ordinary differential equations arise when we are considering the variation 
of a physical quantity that depends on just one independent variable. If we 
denote the independent and dependent variables by x and y, respectively, 
then an ordinary differential equation may involve functions of x and y, and 
the derivatives of y with respect to x. Thus 


Pg d. 
Ce + aft +y=0, -e + cy = 0 (where c is a constant) 
dy 44 F 
and —> =—k’*siny (where k is a constant) 
da? 


are all examples of ordinary differential equations. (Quantities such as ¢ 
and k in these illustrations are sometimes known as parameters; while the 
solution may depend on their values, they are generally assumed to be fixed.) 
Sometimes we omit the word ‘ordinary’ when it is clear that we are referring 
to an ordinary differential equation. 


On the other hand, partial differential equations arise when the depen 
variable is a function of two or more independent variables. If we dehote 
the independent variables by x and t, and the dependent variable by y, 
then a partial differential equation may involve partial derivatives of y with 
, Py _ 92 
respect to « and ¢, as well as functions of x and t; for example, a = oss 
ar 
(where ¢ is a constant). We shall be concerned with such equations and their 
solutions later in the course. 
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Here we are primarily concerned with the process of finding solutions, not 
with the difficulties that may arise in the definitions of the functions and 
their domains. We assume throughout that unless otherwise stated, the 
functions involved are well behaved: for example, we assume that they are 
differentiable as many times as we please on an interval of the real line. (You 
may take this to mean that the functions are ‘sufficiently smooth’.) 


Solution of a differential equation 


A function that satisfies a given differential equation is known as a solution 
of the differential equation, and, while it may not be obvious how one might 
find such a function, it is easy to verify that a given suggestion is correct. 
For example, we can easily verify that the function y = 3? + Ae* satisfies 


ly 
the differential equation a —y = 6x — 3x? for any value of the constant A, 
ve bid 
by writing 


dy 


—y = (6x + Ae*) — (3a? + Ae*) = 62 — 327. 


0 
A solution may be required to satisfy certain additional constraints, for 
example that y(0) = 0, and such conditions will further restrict the choice 
of a suitable function. 


It is not always possible to find the required function explicitly; we may have 
to be content with an implicit definition. For example. it is easy to verify 
that if y = y(x) is a function for which x? + zy + y° = ¢, then differentiating 
with respect to x gives 


d, d; 
Qe +y ta ae 5yi =0. 
dx da 
d, +22 F - seat y 
It follows that 2 = —¥ , 80 y = y(x) is a solution of this differential 
dx x+5y! 
equation; but we cannot find an explicit formula for y in terms of x. We may 
also refer to an implicit definition, such as 2? + ay + y° = ¢, as a solution of 


a differential equation. 


Example 1.18 


Show that the function y(x) = 3sin(A.x) is a solution of the differential equa- 


tion oS = —,*y, where ) is a constant. 
ar 


Solution 


In order to establish that this is the case, we need only differentiate y twice. 


to obtain first ay = 3\cos(Ar) and then 
2 
oy = 3)? sin(Ar) = — 


as required. 


?y, 


In the previous example the solution was given as a function, but more often 
such solutions arise implicitly, as in the following case. 


Example 1.19 


: C 5 ' a . @ =e 
Show that zy = e” + e” is a solution of the differential equation oe = y 4 
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Solution 

Differentiating the equation ry = e* + e” with respect to x, we have y+ ry! 
=e"+e%y', Rearranging this equation gives y/(e! — y—e*, which 
establishes the required result. 


Notice that in this case we are not able to rearrange the original equation 
into a simple algebraic expression for y in terms of z, 


Order of a differential equation 


The order of a differential equation is the order of the highest derivative that 


@; 3 
appears in it. For example, the differential equation a + (#) +y=0 
Es 


is of order two. 


Exercise 1.28 


(a) In each case, form a differential equation of first order by eliminating the 
constant a, 


(i) y=e" (ii) y=sin(x +a) 
(b) Form a differential equation of second order by eliminating the constants a and 
wt+a 
b from y = z 
ath 
. = es 
(c) Solve the differential equation qa 
he’ 
d) Show that a function y = y(a) that satisfies the equation ry + cosy = 3 is a 
y=y og 


y 


solution of the differential equation —— = ———. 
dx siny-—2x 


Linear differential equations 


A linear ordinary differential equation of order n is an equation that 
can be written in the form 


d'y 
dx" 
with a,(x) 4 0, for some function f(a). 


m1 


(1.8) 


d 
an(2) +--+ ay(2)SE + ao(e)y = fl), 


+ Qn-1(2) aa 


A solution of a linear differential equation of order n containing n indepen- 
dent constants is known as the general solution of the equation. For linear 
equations it can be shown that every solution can be obtained by assigning 
values to the independent constants in the general solution. This is not the 
case for non-linear differential equations, 


The differential equation (1.8) is said to be homogeneous if f(x) = 0, 
and such an equation has n independent solutions y;(). y2(x),---.Yn(x)- 
Moreover, every solution y(a) of the differential equation may be written in 
the form 


y(x) = Ary (x) + Azye(x) +--+ + Anyn(2), (1.9) 


where Aj, Ag,...,A, are constants. 


In order to find the general solution of equation (1.8) when f(x) 4 0, we 
need only find a single solution without arbitrary constants (known as a 
particular integral), then add a function of the form (1.9). 
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From the chain rule, the 
derivative of e”") with 
respect to x is y’e”. 


With the aid of a computer 
program, it is possible to 
convince oneself that the 
equation wy = e7 + e” 
defines a function y = y() 
with the negative real axis 
as its domain. 


Part (c) illustrates the fact 
that we might expect the 
solution of a second-order 
equation to contain two 
arbitrary constants. 
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In order to solve a practical problem, it is not usually sufficient to determine 
the general solution of a differential equation; often we need to assign specific 
values to the constants that arise in the general solution. Such values are 
generally determined by certain additional restrictions; loosely speaking, 
one would expect that n pieces of information would be needed to uniquely 
define the solution of a differential equation of order n. 


As a simple illustration, y = Ae” + Be” is the general solution of the 
second-order differential equation y’ = y. The conditions y(0) = 1 and 
y'(0) = 0 are sufficient to determine the two constants A and B (in fact 
A= B= }). This is an instance of initial conditions, which we shall discuss 
shortly. 


1.3.2 Elementary solutions of ordinary differential 
equations 


Separable equations 

A first-order differential equation is said to be separable if it can be written 
d; ‘ ‘ 

in the form 4 = g(x)h(y), and the general solution of the equation may be 


ae 
obtained by evaluating the integrals in the equation 


lim" on | o(oyae. 


Example 1.20 
To obtain the solution of the equation y/ = «xy, we rewrite it in the form 
y'/y =x, then integrate to obtain | (y//y)dx = | dy/y = frac. Evalu- 


ating the integrals gives In y = 3a? +c, where ¢ is an arbitrary constant. In 
fact, the solution is made simpler by choosing the arbitrary constant to be 
c= A, where A > 0 is also arbitrary, for then Iny = $a? +n A, and the 
solution is y= Aexp (32). (This is a solution of the original differential 
equation for any choice of the constant A, i.e. it is the general solution.) 


Example 1.21 


3 . Te: e 
The general solution of the differential equation — may be found 
dx 1+y? 


by evaluating / (1 +y*) dy = fe dx to give y+ 4y* =e" +c. This is an 
implicit equation for y. = 


The solution of a differential equation that satisfies certain additional con- 
straints is known as a particular solution, 


Example 1.22 

To find a particular solution of the differential equation vt +a =0 for 
which y(0) = 3, we first recognise that the equation is separable, then find 
the general solution from the equation ti ydy = — i adx. Thus the gen- 


eral solution is x? + y? =, and it only remains to find the value of the 
constant c. Putting y= 3 when rx =0, we have c= 0? +3? = 9, so the 
particular solution is 2? + y? = 9. 
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Alternatively, we may notice that the solution may be written in terms of 
definite integrals so that the initial conditions are automatically satisfied. Note that in an integral 


7 7 u 
We have uh ydy = — f adx, and evaluating the integrals gives the re- such as [ f(y) dy, we use 
s 8 the same variable for the 
integrand and the upper 
limit. This expression 
Example 1.23 should be interpreted as 


quired result directly. This is a useful trick that is worth remembering. 


A boating pond is 10m wide and 20m long. A boy stands in one corner; his ( i's f(x) ar) 
toy boat, attached to a string 10m long, is at an adjacent corner, as shown 
in Figure 1.8. He walks along the long side of the pond and tows the boat, 
keeping the string taut. Find the equation of the path followed by the boat. 


=u 
2=3 
The ‘shorthand’ is 
frequently useful because it 
saves on substitution of 
variables. 


10m 


Figure 1.8 The path of a toy boat 


Solution 


We use axes as shown in the figure. When the boat is at the point (x,y), the 
boy has moved a distance x + \/100— y?, and the slope of the hypotenuse 
of the shaded triangle in Figure 1.8 is —y/\/10? — y?. Thus the gradient of 
the path followed by the boat is 

dy y 


dx /i00—y?" 
—y 
_fa- [Ba 


Be) 


The integral on the right-hand side can be evaluated using the substitution The integration after the 


u? = 100 —y? for u > 0, which gives substitution uses the 
= rec be » method of partial fractions 
/100 — y? uU u uw d as discussed in 
Sari dy = ale, du) = | a9 Subsection 1.2.5. 
Now notice that y = 10, and therefore u = 0, when x = 0. So, using partial ~s p+ aes 


fractions and the trick mentioned in the previous example, we have gr-lec 


ie u w 7 1° i 
a A pe = 2) ee 
[aw i 2-100 I (: 5(goutwes)) a 


Thus 
=x = [u+5In(10 — u) —51n(10 + w)]% =u+sin(it> “). 
10+u 
which gives 
r=5in( 32) —u ~sn( 24 aa) — ¥ 100 
10-—u 10 — \/100 — y? 


In this example we cannot rearrange the equation to give y = y(x), but we 
could use a computer program to plot the graph of x = x(y) as shown in 
Figure 1.9. 


34 Chapter | Mathematical methods and differential equations 


10 
8 
6 
4 
2 
0 3 10 is Pe 


Figure 1.9 The graph of « = 2(y) (with the independent variable y on the vertical 
axis) 


Homogeneous equations of first order 


The word ‘homogeneous’ is commonly used in mathematics in several dif- 
ferent senses. Here we define a first-order differential equation to be homo- 


A d 
geneous if it is an equation of the form a at s(2). 


In order to solve such an equation, we write y = xv and endeavour to find a 


differential equation in v and x that is separable. Differentiating with respect 


. dy dv eae a + 
to x, we obtain a re + v, and eliminating y from the original equation 
dx da: = 
dv _ f(v)—v 


2 dv 5 Li 
gives c— +v = f(v), which can be rearranged into the form — 
dx dx xr 


which is separable. 


Example 1.24 


In order to solve the equation 
dy _ «+y? 
de oxy 

we first notice that it is a homogeneous first-order differential equation, since 

it can be written as 
dy = 
eae ey. 
a” ff = x 

So we make the substitution y = xv, giving y/ = wv’ + v as above. Elimi- 

nating y gives 


and the differential equation becomes 
dv ar 
r= S sO frw= fF. 
dx ov E 
Hence ju? =Inx +A for «x > 0, and finally replacing v by y/x gives 


(4)’ =2me+c, 


where C = 2A is an arbitrary constant. 


Exercise 1.29 


In each case, show that the equation is homogeneous, and solve it. 
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2 opty Sh 
a) Y= b) Bry? = oF + y3 
(a) ¢ (b) 3ay?—* = a3 + y/ 


Exercise 1.30 
Show that the differential equation 


dy —t Ve? +? 
dx y 


is homogeneous, and use the substitution y = rv to solve it, [Hint: Try the further 


substitution u= V1 + v?.} 


Linear equations of first order 
The general linear equation of first order takes the form 


dy 
ay(x) +ao(x)y = f(x). 
For example, the equation 


3dy 9 5 
3 met 5 
a 4+ 3rry = 2 
dx a 
is of this form, and in this case it is easy to find the solution if we notice 


that 


ody 23 da 
E ie + 3a°y = ie? y)- 
Then the differential equation becomes 
d 7.3.) 
Fale y= 


and we have only to integrate this expression to obtain the solution 

ay = da® +e. 
However, it is not immediately clear how we would proceed if the same equa- 
tion were to arise in a slightly different form, for example, “ + 2ry = Gr. 


Fortunately, there is a simple method that will transform such an equation 
into a form that can be integrated. 


The integrating factor method 


Any linear differential equation of first order of the form 


d; 
& + gla)y = h(e) (1.10) 
may be solved by multiplying both sides of the equation by the integrating 


factor 


P(x) = ex( f a(e)ar) a 


The differential equation then becomes 


a 
ay (yP(#)) = h@) PG), 


which may be integrated to give the general solution 


1 
y= Pa) [beyPCeyae. 
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Example 1.25 
Find the general solution of the differential equation “ + 2ry = Gr. 


Solution 


This equation is of the same form as equation (1.10) (i.e. first-order linear), 
with g(x) = 2x and h(x) = 6x. The integrating factor is 


P(x) = eo(f 2a ar) = exp(z”). 
So the general solution is 
(oe sy ie ee 2 
y= oe) [ srexote )dr = ae) (exp(x?) + A) 
=3+3Aexp(-z?). @ 


Example 1.26 


Find the particular solution of the differential equation 


dq q Differential equations of 
Ra a= E(t), this kind arise in the 
i ; theory of electrical circuits, 

where R = 10, C = 10-4 and E(t) = 20sin(100zt), that satisfies g(0) = 0. 
Solution 
The differential equation becomes 

10% + 10%q = 20sin(100zt), 
which simplifies to 

dq 


ay + 100a = 2sin(100zt). 
This equation is first-order linear (see equation (1.10)), with g = 100 and 
h = 2sin(100zt). So the integrating factor is P(t) = ex( f 1004) = elt, 


and multiplying the differential equation by the integrating factor gives 


£ (6!) = 2e! sin(100zt). 
Integrating, we obtain We discussed integrals of 
1008 the form t sin(bt) dt 
e1%tq — —“ ____(sin(100nt) — cos(100zt)) + C. sh a 
50(1 + 7) page 26. 


Using the condition g = 0 when t = 0 gives C = 7/50(1+ 7°), and finally 


1 
=———.- (sinf = 1000) 
q +a) (sin(1007t) — mcos(100xt) + me ). a 


Exercise 1.31 


d 

a) Find the particular solution of the differential equation &Y + ycotx = 3e°* 
is dx 

that satisfies y = —1 when x = 3. 


‘ind the general solution of the differential equation (2 — 1)(a2—2)—+y 
b). Find the general solttion of the differential equisti wy 


(c) Show that the transformation v = y'~”, where n # 1 is a constant, reduces 


the differential equation 2 + yF(a) = y"G(zr) to linear form. 
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Linear equations with constant coefficients 


Linear differential equations with constant coefficients are particularly easy 
to solve. We shall concentrate on equations of order two, but higher-order 
equations can be solved in a similar fashion. (For order greater than four, 
the auxiliary equation cannot generally be solved algebraically.) 


We consider differential equations of the form 


@y 


are +00 +6 = f(z); 


and we assume that you have seen such equations in a previous course; these 
brief notes are intended primarily as revision. 


where a, b and ¢ are constants, 


Procedure 


The above differential equation is said to be homogeneous when f(x) = 
In order to solve the homogeneous equation 
_ + ott +cy=0 

where a, b and care real constants and a # 0, we use the following procedure. 

e Write down the auxiliary equation am? + bm +c = 0, and find its roots 
my, and mo. 

e Ifthe auxiliary equation has two distinct roots, then the general solution 
is y= Ae™* + Be™*®, 

e If the auxiliary equation has two equal roots, i.e. m, = mg, then the 
general solution is y = (A+ Br)e™*. 

e If the auxiliary equation has a pair of complex conjugate roots, given 
by m =a+i8 and mg =a —i, then the general solution is 
y = e®" (Acos Sx + Bsin 3x). 


(In each case, A and B are arbitrary constants.) 


In order to find the general solution of the second-order inhomogeneous 
linear equation with constant coefficients 


dPy dy 
ag +0 + oy = f(x), (1.11) 
da’ 
we first use the procedure above to find the general solution y-(x) of the 
P 
associated homogeneous equation om + oft +cy=0, Then we find a 


particular solution yp(a) of equation (1.11). Here y.(x) is known as the 
complementary function, and y)() as a particular integral. 


The general solution of equation (1.11) is then 
y(x) = ye(x) + yp(2). 


For an nth-order differential equation, y will be a solution containing n 
arbitrary constants, so it will indeed be the general solution. 


There are various techniques for finding particular integrals yp, and we refer 
you to any of the standard texts on differential equations for the details. 
Here we illustrate one of these techniques, known as the method of un- 
determined coefficients. In order to apply this method, we need to be 
able to guess the form of a solution and then determine the constants in our 
proposed form by equating coefficients. The following example illustrates 
the method. 
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Note that the term 
‘homogeneous’ as used here 
differs from its application 
to first-order equations on 
page 34. 


If the coefficients a, b and ¢ 
are real, and the auxiliary 
equation has complex 
roots, then these roots will 
always occur in complex 
conjugate pairs. 
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Example 1.27 
To find the general solution of the differential equation 


Py dy 

m 5 oi + Gy = 10, 
we consider first the auxiliary equation m? — 5m +6 = 0, which has roots 
m, = 2 and mg = 3. Hence the general solution of the homogeneous equa- 
tion (i.e. the complementary function) is y. = Ae?” + Be**. In order to find 
a particular integral, we guess that a constant function yp(2) = k may be 
suitable, and substituting this into the original differential equation gives 
It follows that the general solution of the differential 


equation is 
Y = Ye + Yp = Ac™ + Be™ + 3. 
This is the general solution of the second-order differential equation because 


it satisfies the equation (this can be checked by substitution), and it contains 
two arbitrary constants. M 


If the number 10 in the previous example were to be replaced by 10 +32 +2, 
then we might guess that a+ bx +exr? is a suitable form for a particular 
integral. 


The method may also be extended to differential equations of higher order. 


Example 1.28 

1 
Find the general solution of the equation a -y=0. 
Solution 


This equation is homogeneous, so it is not necessary to find a particular 
integral. We consider the auxiliary equation m!—1=0: factorising the 
left-hand side gives 


m! — 1 = (m? —1)(m? +1) = (m—1)(m +4 1)(m—i)(m +i), 
so the general solution of the differential equation is 
y = Ae” + Be~® + Ccosx + Dsinx. 


This is the general solution of the fourth-order differential equation, because 
it satisfies the equation and contains four arbitrary constants. Ml 


Exercise 1.32 


(a) Find the general solution of the differential equation 


dy dy 
=e 60: 
ae? Fay UO 
(b) Given that y; =—ze* is a solution of the differential equation 
@y dy = 
da ae 


find the general solution. 


In each of Exercises 1.33-1.35, a particular integral may be found by using 
the method of undetermined coefficients. Try to guess the form of a solution, 
then equate coefficients to find a particular integral. 


1.3 Differential equations 


Exercise 1.33 


Find the general solution of the differential equation 


ey 
da? 


dy 
dx 


Exercise 1.34 


Find the general solution of the differential equation 


+5— +4y = 3-22. 


@y 
da? 


+ 4y = 2x. 


Exercise 1.35 

Find the general solution of the differential equation 
fy . ody . 
2-5 +2 + 3y= Qn 1. 
a tq, +3 2+ 22 — 1. 


Exercise 1.36 
Find the general solution of the differential equation 


dy 
= ty=0. 
drt 7 


Initial conditions and boundary conditions 


We have previously discussed cases where a particular solution of a differen- 
tial equation is defined by the value of the function at a specific point, but 
sometimes we require more than one additional piece of information. As an 
illustration, consider the equation defining simple harmonic motion 

& 

a =-wy, (1.12) 
where w is a specified real positive constant, y denotes the displacement of 
an object, and the independent variable t represents time. The general so- 
lution of this equation is y = Acoswt + Bsinwt, and we require two further 
conditions if we are to determine the values of the constants A and B. Of 
ten such conditions arise in one of two forms: initial conditions, in which 
the value of the function and its derivatives are specified at some chosen 
value of t, or boundary conditions, in which the function is defined at 
two different values of t. Boundary conditions will be dealt with later in the 
course. 


Suppose that an object performs simple harmonic motion and satisfies equa- 
tion (1.12), and that we are given the position yo and velocity vp of the object 
at time ¢ = 0; in other words, we are given initial conditions for the motion 
of the object. We can e: determine that A = yo and B = vp/w, and this 
provides us with the particular solution 


y(t) = yo coswt + (vp /w) sinwt. 


Here initial conditions are sufficient to determine a unique solution for any 
choice of yo and vp. This is generally the case; provided that we are given 
an appropriate number of initial conditions, the solution is always uniquely 
determined. An nth-order equation requires n initial conditions. 
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End-of-section Exercises 


Exercise 1.37 


Solve the equation # = e*+¥, given that y(0) =0. 
ke 
Exercise 1.38 
dj 
Solve the equation z = —2rtany, given that y(0) = 5. 
Exercise 1.39 
di 
Solve the equation = = 5; given that y(1) = 1. [Hint: Show that the 


equation is homogeneous, and use the method of partial fractions.) 
Exercise 1.40 


d ne iP 
Solve the equation «<4 = y(1+Iny—Inz). [Hint: Show that the equation is 


homogeneous, and solve the resulting integral by a suitable change of variable.] 


Exercise 1.41 


ly \* d 
Solve the equation x (#) = ay -r=0. [in First solve a quadratic equa- 
ae a 
tion) then tememb eerily eS (ernie) ae 
ion, then re that “(5 eS | 
dx Vvi+a2 


Exercise 1.42 


dy 


ae + 12y = 12x? + 10x — 11. 


& 
Find the general solution of the equation ——> 
d 


Exercise 1.43 Harder exercise 


(a) Find the family of curves for which the length of the part of the tangent between 
the point of contact (x,y), in the first quadrant, and the y-axis is equal to the 
y-coordinate of the intercept of this tangent with the y-axis. (In other words, 
find curves where the lengths OA and AP shown in Figure 1.10 are equal.) 


Figure 1.10 


(b) Find the equation of the orthogonal family of curves (i-e. the curves such that 
the product of the gradients at the point of intersection with each member of 
the previous family is equal to —1). 


Exercise 1.44 Harder exercise 


A chemical compound C is being formed by the reaction of two substances, A 
and B, During the reaction, substance A combines with substance B in the ratio 
2:1 by mass to form substance C, The rate of formation of C is proportional to 
the product of the masses of the quantities of A and B that remain uncombined. 
Initially, there are 200 grams of A, 100 grams of B, and 0 grams of C in the mixture. 
Given that 100 grams of substance C are formed after 30 seconds, find the amount 
2(t) of substance C formed after time t. 


1.4 Outcomes 


1.4 Outcomes 


After studying this chapter you should be able to: 

e understand the terms continuous, differentiable and order when applied 
to functions, and apply the O notation; 
understand what is meant by the convergence of an infinite series; 
differentiate expressions involving elementary functions: 
evaluate the Taylor series for a function, and understand what is meant 
by its radius of convergence; 

e integrate elementary functions; 

e understand what is meant by the convergence of an integral over’ an 
infinite range; 

e express rational functions as partial fractions, and apply this technique 
to the integration of rational functions and the expansion of rational 
functions in seri 
manipulate expressions involving complex numbers; 
transform between the Cartesian, polar and exponential forms of a com- 
plex number; 

e calculate powers and roots of complex numbers using de Moivre’s theo- 
rem; 

e understand the terms ordinary, partial, linear, order, general solution 
and particular solution when applied to differential equations; 

e solve first-order ordinary differential equations which are separable, ho- 
mogeneous or linear; 

e solve higher-order ordinary differential equations which are linear with 
constant coefficients; 

e apply initial conditions to general solutions, to fix the arbitrary con- 
stants. 
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Solutions to Exercises in Chapter | 


Solution 1.1 
We have 
2 
R,(a) = Yak =a" +a"! 4g"? 4... =a"(1+ata2+---). 
k=n 
For —1 <a < 1, using equation (1.2) gives 


a” 


R,(a) =a" lim S,(a) = 
hk—+90 1 


—a 


Solution 1.2 


From Equation (1-1) we have 


l+atatt+e 


Replacing x by — 
N 


5 (ci 


k=0 


which converges to the limit Ge” N — ow if |2| < 1, but diverges if |x| > 1. 
ze 


Thus the infinite series $7 (—1)'x** converges if |r| < 1, but diverges if || > 1, so 
) 


0 
the largest interval |.c| < p corresponds to choosing p = 1. 


This exercise is one of the very rare occasions when we are actually able to determine 
the function f in a simple form from the coefficients {a}, for in this case 


1+ (—1) 8224+? 
1 +27 


fle) = in| 


Solution 1.3 


From the Taylor series we have sina: = x + O(c) and In(1 +x) = + O(a). It 
follows that 


sina In(1 +2) = [a + O(2*)][a + O(a?)] = 27 + O(2*). 


Solution 1.4 


Solution 1.5 


*y 1 : 
(0) Ine'= i 4 dt and + <1 for ¢ > 1. From the stated result we therefore have 
fi 


if i dt<z—1<z. 


(b) The formula for integration by parts may be written in the form 


Ina = 


[weye'@ae= u(x) v(x) — [ weyea)ar. 


Solutions to Exercises in Chapter | 


(i) Choosing u(a) = Inz and v(x) = Inz, we have 
Ine ge = (nz)? [ax 
J & x 
which can be rearranged to give 
In? ge = }(nz)?. 
© 


Alternatively, make the substitution s = Inz; then 


[ae = [sas = 4s? =3(ne)?. 


T 
Ing 
f = te = a(n). 


(ii) Choosing u(a) = In and v(x) = —r~', we have 


T Ina nT 1 1 nT 
i eiges: (3-7) -cy-1-3-. 
(c) (i) Using the result of part (b)(i), we have 
z 
ta Lie sanT)?, 
ats 
which increases without bound as T increases, so the integral is divergent. 


(ii) From part (b)(ii), 
T Ina 1 MT 
a dp ax Pais I 
f[ Be-1-3 
1 
so the integral is bounded. Since also =F > 0, the integral is convergent. 


Solution 1.6 


For 1 <a < T, we have 0 < = < ze-*", 80 using the stated result gives 


ri Tr 
[vas f ae de =| 
1 1 


© 0 
Since this integral is bounded, and e-*” > 0, we see that the integral [ e* dr 
1 


oa 
is convergent. Using the hint, the convergence of [ e-* dr follows immediately. 
0 


Solution 1.7 
(a) We write 
lia+5 = A a B 
(Qe —1)(3a+2)° 2r-1  3r+2° 
Then putting x = 0 gives —5/2 = —A+ B/2, and putting x = 1 gives 16/5 = 


A+ B/5. Solving this pair of equations for A and B, we obtain A = 3 and 
B=1,s0 


ee 
(Qz—1)(8@+2) 2-1 ° 324+2° 
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(b) We write 
Ac+YQr+7) __A_, B Cc 
(@@—D(e+Det+3) x-1° r+2° x+3 


By substituting suitable values for x, or multiplying both sides by the de- 
nominator (a — 1)(# +2)(a +3) and equating coefficients of powers of x, we 
find 

A=3, Be2, C==1. 


(1.13) 


There is a third method which is often useful. Suppose that we multiply both 
sides of expression (1.13) by « — 1, to obtain 


Yet Qe+7) _ 4, B@e~1) | Cle—1) 
(@+2)(@+3) — z+2 +3” 
and then put x= 1. We obtain A= tet = 3 immediately. Similarly, 


multiplying (1.13) by 2 + 2 and putting « = —2. and multiplying by 7 +3 and 
putting « = —3, gives ‘B and C at once. 


(c) You may have noticed that in the given expression, the numerator is not of 
lower degree than the denominator, so our first step must be to overcome this 
difficulty. This can be achieved by writing 


a3 +20 (a(e?-1) +2) +20 


2-1 wl 
_ @(e? —1) +32 
~~ ie=1) 

3a 


=*+ Ger) 
Proceeding as before, we find 
w+2r 3 5 3 
a—1 °° &e-1) © etl) 


Solution 1.8 
r 3 2 


@) Bote o=8 22 
(b) If 
ee oer ered! 
@-17 ~** @@)’ 
then 
Po)? 1) Be 1 2 u 
Q@) ~ @=1? > =F @=1F 1 
sO 
ee es 
@-1p '* @-) * @-1 
(@-1? 4 1 
©) G-ap@=2) =o" 2-2 
(a) a—-2 _ _2+2 3 1 


Q+2)1+2? 32422) 30+2) (+2)? 


Solution 1.9 
a+bi _a+bi ~ ca di _ 
ce+di ctdi c—di 


3 15 ((ac-+ bd) + (be — ad)i) 
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Solution 1.10 

= = —16 + 12i (and notice that it would be equally acceptable to write z = 12i — 16). 
le] = (16)? + 12? = 442 $3? = 20 and Im(z) = 12 (not 12i). 

z* =—16 — 12i and z2* = |z|? = 400. 


Solution 1.11 


From de Moivre's theorem we have (cos @ + isin @)° = cos50@ + isin 5@. So, from the 
binomial theorem and using the fact that cos? @ + sin? @ = 1, we have 


cos 5 = Re((cos @ + isin 8)°) 

= cos’ @ — 10.cos’ @sin? @ + 5cos@sin* @ 
cos 0 [cost @ ~ 10 cos @(1 — cos? ) + 5(1 — cos* @)?] 
cos @ (16 cos* 9 — 20 cos” +5) . 


Solution 1.12 


(a) First we calculate |2| = /2? +2? = 2V2, then if we plot the point z = 2 —2i 
on an Argand diagram, we can easily see that Arg(z) = —7. Thus, in polar 
form, we have 


2 = 2v2 [cos(—%) + isin(—3)] . 


Any value of the argument of the form —} + 2k7, where k is an integer, is 
equally acceptable. 


(b) The general polar form is z = 2V2 (cos(—¥ + 2k) +isin(—} + 2km)), for k 
an integer. Hence, using de Moivre’s theorem, we obtain 


21/3 — 81/6 (cos (1(—% + 2km)) + isin (+(—% +2kz))) . 
For k = 0,1, 2. this gives three distinct cube roots: 
29 = V2 (cos(— 4) + isin(—)), 
21 = V2 (cos 43 + isin 7) , 
29 = V2 (cos + isin 5) . 
The remaining values of k simply repeat these roots. 
Solution 1.13 


For any real number n, 


2” = (re”#)" = r"el™ = r"(cosnd + isin nd). 


Solution 1.14 
We have 
Ja +b]? = (a+ b)* = a? + + 2ab < Jal? + [b|? + 2Jal|b| = (ja) + {o))? ; 


and the required result follows. 


Solution 1.15 


at wee 
4 (cinh.x) = g (- = )=£ +e =coshz, 


d: dx 2 2 
d d (@+e* emer. 
qq (cosh) cir ( 2 ) = 2 = sinh x. 


If d= Asinhwx + Bcoshwa. then # = Awcoshwr + Bwsinhwr, so 
ao 


—; = Au® sinhwr + Bu* coshwx = wo. 
dx? 
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Solution 1.16 


Weemente: ~~ a4 — 2 Aetna 
ROED) no EL 
a 
4 n(n +1) NT n+1 


This follows because each of the negative terms, except the last, cancels with the 
first term in the next bracket. Thus 


© 1 N 1 
DD, n(n+ 1) stim, Tee ay (: 


n=l n=1 


Solution 1.17 


We know that 


2n+1 


et siieet 
Qn+Hr 


gent 
“Fara 


2 e- 
sinha = 


sinh(2) 


Solution 1.18 


From the Taylor series we have 


a 
e-—+0(r°) and =1+27+ O(c"). 


sing a 


Hence 
sine 
1—a4 


= [x — 42° + O(2°)] [1 +27 + O(z")] 
=r—a+2°+O(2*) 
=2+ $27+0(2*). 
Solution 1.19 
We have 

j= 2S 
ioe 1+(-«)? ~ 


so the function is odd. Thus 


[ : fa)ae= f " fle) de + i: : fla)de 


1 Oo 
= r = =z) d: 
ff seac— fi s-ayae 
! 


For any odd function f(a), and any real number a, it is obvious from the graph 
of f that if F(x) dx =0. 
~a 


Solution 1.20 
If w=sinh~! x, then « =sinhw. Differentiating with respect to 7, we obtain 


1=coshw—. 


dx: 


Here we have substituted x 
for —x in the second 
integral. 
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It follows that 
dwt 1 eat 
dx coshw” /j4sinh?w Vite? 
The graph of y=sinh™'« can be obtained easily by reflecting the graph of 
y =sinhr in the line y = x. 


Figure 1.11 The graph of sinh”! x 


We have V1+ 2? > Vr? > —x for alle € R,s0r+ V1+22 > 0. 
Now 

(in («+ 1+a?)) = : % 
der Vi+a2 


So we see that a = A (In (2 + V1-+22)), therefore w and In (x + V+?) must 
differ by a constant, giving 


sinh! x = In (x +V1+ 7) +c. 


However, sinh! x = In (x + V1 +22) =0 when x = 0, so ¢ = 0, and this estab- 
lishes the required result. 


Solution 1.21 
We can write 
4c7—-9r+4 _ 2 
wx(e@—1)(e—2) 2 
so for « > 2 we have 
4a? — 92 +4 2 
lea” eit 
= 2Inx + In(x — 1) + In(w — 2) +e 
= In (x?(a — 1)(2 —2)) +e. 


Tt follows that 


4 2 
ri ene de = [In (x*(x —1)(2—2))]} = (16/3). 


Solution 1.22 


For a # 1 we have 


- 
[ ze- 1 fay? r*_ yy, 
1 


te l-a 1 1=a 
Se ip Easel oyhnle Sofa <1 thes indy cine notiexset: 
T—00 J, t a-1 


os 
1 
We also showed in a previous example that i i dt is divergent. 
1 


nea | 
Tt follows that [ = dt is convergent for a > 1 and divergent for a < 1. 
1 
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Solution 1.23 
First we notice that 0 < f(t) + |F(t)| < 21f(t)|, so 


iT ‘T oo 
f cosines f yoase [sole 


T 
Thus the integral - (f(t) +|F(t)|) dt is non-decreasing as T increases, and it is 


a 
bounded above, so the limit 


T 
sim f(s +\n00)) a 


T 
exists. But we are given that jim, i |f(t)| dt exists, so 
00 Jig 


T <a ‘T 
im [ f(O)at = im [/ (1) +isto) ae f lo 


20 
exists. In other words, [ S(t) dt is convergent. 
a 


_ 1 sac) er ; * |sint 
Since < BE (for t # 0) and 2 dt is convergent, it follows that roe dt 
1 1 
P sint .. 
is convergent, therefore = dt is also convergent. 
1 
Solution 1.24 
(a) i-1 (b) $+2 (ce) 1 
Solution 1.25 
(a) 1-4 = VBe-""/4, go b= 1 eters, The polar forms of these 
1-4 2 complex numbers can be 


wy 0/9" = er? te a a ar 
(c) (1+iv3) = (2ei"/3)'* = Qldglbin/3 — 9l5g5im — gl5eix 


Solution 1.26 

e+e — (cos + isin®) + (cos@—isind) _ 
=: 2 * 
ec” _ (cos + isind) — (cos ~isind) _ 


ree 
sinh(i9) = © 5 5 =isind. 


Solution 1.27 

jz + wl? +|2—wl? = (2 + w)(z + w)* + (2 -—w)(2—w)” 
=(z+w)(z* +") + (2 —w)(z* —w") 

+ 2wu* 

= 2I2|* + 2|w|” 


cos @, 


cosh (id) = 


Solution 1.28 


(a) (i) If y=e**, then ny = ax, and differentiating with respect to « gives 
a=y'/y. Substituting y//y for a, and using the fact that y # 0, we 
obtain the first-order equation 

ylny=-y'. 
(ii) Differentiating gives y’ = cos( +a); using the identity cos? @+sin? @ = 1, 
we have 


GP +y? =1. 
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(b) 


(c) 


(d) 


Since (a + b)y = 2 +a, differentiating with respect to x gives (x + b)y’ +y = 1. 
Differentiating again, we have (x + b)y +2y’ = 0. Eliminating x +b from these 
two equations gives 


ay)? =(y—1)y”. 


a = 1 twice, we obtain first Ed =x+A, then 
y = $2 + Ax + B, where A and B are arbitrary constants. 


Integrating the equation 


Differentiating the equation xy + cos y = 3 gives ae +y—sin yt = 0, which 


dx 
can be rearranged to give dy a, 
dx siny—x 


Solution 1.29 


(a) 


(b) 


The equation can be written as 


dy 1+#% 
dx 1-2" 


so it is homogeneous. 


Putting y = av as before, and eliminating y, we obtain the equation 


which can be rearranged to give 


[Ga-a)e -{F 


Integrating gives 
tan“! vy — $In(1 +0?) =Inw + A 


for « > 0, and replacing v by y/2 gives 


2 2 
2tan (4) — in (2 a ) =2ine+B, 


2 
that is, 
2tan—"(y/a) —In(x? + y*) = B. 


The equation can be written as 


2 
ae hey, 
dz 3\y 3\r 

so it is homogeneous. 


Putting y = av, the equation becomes 


dv a! 
eee 
which can be rearranged to give 
ot a2 — 2v 


dx 30 
Rearranging again and integrating, we have 


i - [*%. 
1-2” 
hence 5 fe —2v3) = —Ina +} In A ( for x > 0, 20° < 1), which can be rewrit- 


ten more simply as 1 — 2v? = A/z?, so the required solution is 
a? — 2y3 = Az. 


49 


The choice of the arbitrary 
constant in the form }1n A 
simplifies the final answer. 
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Solution 1.30 
The equation can be written as 
dy -1+ VI1+ WP 


dc (y/) 
so it is homogeneous. Putting y = rv, we have 


dv a ritvit+e 
v ” 


which can be rearranged to give 


_ —v?-14 V+? 


and therefore 


f= ie / dv 

a J -+u%)4Vi¢e 

The integral on the right can be evaluated by making the substitution u = V1 + 0, 
and we then have 


Integrating, we obtain 
Ing = —In(—u+1)+InC, 


sO 
c Cx 


[Hayter | aye 
‘Therefore +a — \/x? + y? = C, which can be rearranged to give 
y? = C(C F 2x). 


u 


Solution 1.31 


(a) The equation is first-order linear, and the integrating factor is 
P(x) = exp (/ cot rds) = exp (In(sinx)) = sin. 
Multiplying the differential equation by the integrating factor gives 


4 OD ee eee 
sing de tycst= a sin) = 3e sins, 


and integrating gives 
ysinx = [seorsinzar =e + AL 


Since y = —1 when « = §, we have (—1)sin = —3e"* 3 + A, which gives 
A= 2. so the required solution is 


ysina = —3e°* + 2. 


(b) In order to solve the differential equation Notice that we must ensure 
e that the coefficient of 
(w —1)(e- 2) + y=2?, dy/dz is 1 before using the 
de formula for the integrating 
we first rearrange it into the form factor. It is a common 
2 error to forget to divide 


dy + y= x r the equation 
dx (x—1)(x — 2) (a —1)(x — 2) dy 
- i ay(x)— + ao(x)y = f(a) 
This is of the same form as equation (1.10) on page 35 (i.e. first-order linear), da: 
with by a(2). 


1 1 
a) = G—a)@=2) 2-3 
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so the integrating factor is 


P(a) = exp (fae ar) 
1 
=e f (5-s4)# 


exp (In(a — 2) — In(a — 1)) = 


2-1 
Now we obtain a further rearrangement of the differential equation by multi- 
plying both sides by P(x) to give 


Et ie ae 

z-1) de * @—-1)"" @-1)” 
which can be rewritten in the form 

d (2-2 weet 2: 1 

de \r—1" e=I?" 


where we have ude the pee side in terms of partial fractions. 
Integrating gives the general solution 


a—2 1 
py =2+2n(e—-1I)->T +e 


which can be rearranged to give 


z 1 
y= 25 (e+ 2me- y- r+). 


(c) Differentiating v = y'~" with respect to x gives 
=p 
da’ 
and the differential equation becomes 
Es 
—n . 
Pes can be rearranged into the form 


& + (1—n)y'"" F(x) = (1—n)G(x), 


ae nytt 


2 + yF (x) = y"G(a). 


which then becomes 
dv 
ae (1 —n)vF(x) = (1 —n)G(2), 
which is of linear form (and can be solved by the previous method). 


Solution 1.32 


(a) The auxiliary equation m? +m —6 =0 has roots m; = 2 and mz = —3, so the 
general solution is y= Ae** + Be~**. 


(b) The auxiliary equation m? — 3m +2 = 0 has roots m; = 1 and mg = 2, so the 
general solution is y= Ae* + Be** — xe*. 
Solution 1.33 


‘The auxiliary equation m? + 5m +4 =0 has roots m; = —1 and m2 = —4, so the 
complementary function is y. = Ae~* + Be~**. 


The form of the right-hand side suggests that we should try a particular integral of 
the form yp =a + br. This gives 


& 
or +5 + yp = 5b-+ 4(a + ber) = (da + 5b) + dba = 32x. 
Equating coefficients of x and constant terms gives a = + and b= —}. Thus the 


general solution is 


Y=Yet+Yp = Ae * + Be * +2 —- 
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Solution 1.34 


The auxiliary equation m? + 4 = 0 has roots m; = 2i and mz = —2i, so the general 
solution of the associated homogeneous equation is y.(a) = Acos(2r) + Bsin(2r). 


We now try a particular integral of the form yp) = ax + b, which gives 4(ax +b) = 2x, 
from which we see that a = } and b= 0. Thus the general solution is 


y(x) = Acos(2r) + Bsin(2r) + $x. 


Solution 1.35 


The auxiliary equation 2m? + 2m +3 =0 has complex roots m; = —} + $V5i and 
mg = —4 — 4,//5i, so the general solution of the associated homogeneous equation 
u z7% 

is 


yela) = e~*/? (Acos(} Viz) + Bsin(}v5z)). 
We now try a particular integral of the form yp = ax? + bx + ¢, which gives 

da + 2(2ax + b) + 3(ax? + ba +c) = 2? + 2x —-1. 
Equating coefficients of powers of x gives a = 4, b = 3 and c = —3, so the general 
solution is 


) 


u(x) = e*/? (Acos(} Vix) + Bsin(} Vz) + 42? + 3x ~ B. 


Solution 1.36 


We consider the auxiliary equation m* + 1 = 0, for which the solutions are the four 
possible values of (—1)'/4. Writing —1 in polar form, 
1 = cos(2k + 1)m + isin(2k +1)m, 


for any integer k, and using de Moivre’s theorem, we have 
(-1)'/4 = (cos(2k + 1)m-+isin(2k + 1)x)'/* 
= cos (+(2k + 1)m) + isin ($(2k + 1)r). 
Choosing k = 0,1,2,3, we obtain the four values 
m, = cos % + isin} = a" mp = cos 8 + isin 32 =! 
mg = cos 2% + isin 92 = aie A mg = cos 72 + isin 7= = 


v2 


(It is easy to verify that each of these is a solution of the equation m* +1 = 0.) 


‘The general solution of the differential equation is therefore 
“= (Pos(/v3) + Qsin(x/V2)) etlv2 
+ (Reos(x/v2) + Ssin(x/V2)) e/V2, 
Solution 1.37 


The equation is separable, and we have 


fora [erae. 


Using the initial condition y(0) = 0, we have 


y 2 
feow-f e** da, 
0 0 


so [-e~¥]), = } [e?*]¢, which gives 1—e-¥ = }(e?* — 1), thus 


2 
v=n(=z). 
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Solution 1.38 


‘The equation is separable, and we have 


[cova = -2 [ ras. 


Using the initial condition y(0) = 5, we have 


y 2 
cot y dy = -2 f wdx, 

n/2 0 

so [In(sin y))” = [2?]§. which gives In(siny) = —z?, thus 


y =sin™' (exp(—2*)) . 


Solution 1.39 

Dividing the top and bottom of the right-hand side of the equation by ry, we obtain 
dy 2 
dx 3a/y—2y/x" 


so the Eas Ponce ee Putting y = av gives y/ = rv' + v, and the equa- 


tion becomes «— + v = Rearranging this, we obtain 


da = 


which is separable. Therefore, using the initial condition v(1) = 1, we have 


dx 3 — 20? 
fe -[ wat 


The integral on the right can be evaluated by writing the integrand in partial 
fractions as 


3—2v? v2 v2 3 


v(2v? —1) — ea v2u+1 0 


nh a a 
= =j (A a es i *) Hs 

= [in(v3v — 1) +In(v2v +1) =sinx], 

zs [n (4 ah 

ws (25!) =u (Be) =m (72) 


So the general solution is y* = 2y? — 2?. 


thus 


Solution 1.40 
‘The equation can be rewritten as 
dy 
m2 C*e(G)): 
where the right-hand side is clearly a function of y/z. 
du 
Now put y = xv, so y’ = av’ + v and the equation becomes os +uv=v(1+Inv), 


which reduces to of =vlnv. This equation is separable and gives 


de _ dv 
vine’ 
The ee on the right can be evaluated using the substitution u = Inv, and we 


obtain Inz = In(Inv) +In@. Thus x = C Inv or v = e®* (changing the constant). 
Finally, replacing v by y/x. we have the solution y = xe?*. 
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Solution 1.41 


First we solve the quadratic equation 
dy\* di 
2(#) = —2=0 (1.14) 


to obtain the two possible values 


dy _ yt Var ty? 
ar a z 
This can be rewritten as 


t-te 


so it is homogeneous. Putting y = rv gives 


di 
ro tusvt Vite, 


dx 
so 
def _dv 
r Vi+e2- 


Thus Ine +Ine=+sinh7! v or v = +sinh (In(ex)). Hence 


~In(ex) 


pin(cr) 
y_.e 
Yo 


which gives 
Qey = + (ex? ~ 1). 


It is interesting to note that equation (1.14) also has the solution y = tix, which 
cannot be obtained from the above solution. This is because it is a non-linear 
differential equation. 


Solution 1.42 


The complementary function (i.e. the solution of the homogeneous equation) is 
found by solving the quadratic equation 


m? —7m+12=0. 


The roots of this equation are m, =3 and mg 
is ye = Ae™ + Be!*, 


4, so the complementary function 


In order to find a particular integral, we try a solution of the form y = ax? + br + ¢, 
and when we substitute this into the differential equation and equate coefficients, 
we find that a= 1, b=2ande= 4. 


‘Thus the required general solution is 


y = Ae™* + Be +27 420+ 4. 


Solution 1.43 


d 
(a) Let oa be the gradient at a fixed point (x,y) on the curve, with x >0 and 


y = 0. Then the equation of the tangent line at this point is 
dy iy 
Y-y= ae wr), 


where X and Y vary along the line. When X = 0, we have 
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This equation can be rearranged to give 
dy 
dx” 


2 


y —2ry 


and then 


Qry 
This equation is homogeneous, so we substitute y = .rv to obtain 


ef 
“Je ji+e™ 


Integrating, we obtain — Ina +In(2C) = In(1 + v?), which gives x? + y? = 2Cx 
or 


(e-CP+y=Cc? 
(a semi-circle of arbitrary radius C with centre at (C,0)). 
(b) In order to find the orthogonal family of curves, we must solve the equation 


ay 
fae 
Substituting y = rv gives 


Integrating gives 


dar 1-2? 1 _ 2 
[*-/[ae-[G-)« 


= sine Deh tiainnlitias eee = 
so Ina —1n(2A) = Inv —In(1 + v?), which simplifies to WA Tew and there- 
fore 
x? +(y— A)? = A? 
(a semi-circle of radius A with centre at (0, A)). 


‘The two families of curves are shown in Figure 1.12. 


Figure 1.12 Each curve of the second family is orthogonal to each member of 
the first family 
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Two lines with gradients T’ 
and —1/T will always be 
orthogonal. 
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Solution 1.44 


Suppose that at time ¢ there are z(t) grams of substance C, which must have been 
formed from 22(t)/3 grams of A and 2(t)/3 grams of B. Then 


dz " 22 z 
= =k (200 - =) (100 - :), 


where k is a constant which dictates the reaction rate. This can be rearranged to 


give 
fae = = f ee eee 
9 ~ 9 J (z—300)? 300-2 


Since 2(0) = 0, we see that c = —1/300, thus 
a re 
9 ~~ 300-2 300° 
We can deduce the reaction rate k from the fact that 2(30) = 100, which gives 
CORE INE he ak 
9 — 200 300 600° 
so k = 1/4000, It follows that 
3008 
0 00 


CHAPTER 2 


Multivariable and vector 
calculus 


2.1 Introduction 


As in the previous chapter, our intention here is to provide some of the basic 
tools of applied mathematics. and the emphasis is on methods rather than 
justification and proof. 


ion, but 
The reasons for introducing some of 


Some of the material may be familiar, so you can treat it as re 
there is much that you may find new. 
the techniques may not be entirely apparent at this stage, but their relevance 
will become evident as the course develops. For the moment you should 
concentrate on acquiring the skills that will be needed later, and attempt as 


many of the exer 


ises as possible. 


Some id 
component of every mathematician’s background, but we shall not ask de- 
tailed questions on their use. Other notions, such as partial differentiation, 
will be used throughout the course and you will need to gain some expertise 
before you move on to later chapters. Generally, you will be able to judge 
the importance of the topics by the number of e» ises we ask you to 
tempt. However, the Gauss divergence theorem is an e3 
of crucial importance later. 


such as continuity, are included because they are an essential 


ption and will be 


Subsections 2.2.1 and 2.2.2 are mainly background reading. Subsections 2.2.3 


and 2.2.4 discuss various techniques that will be essential later on. 


material which may be familiar, but they 


ue in Block III. 


Subsections 2.2.5 to 2.2.8 discus 
lead on to ideas that we shall pu 


Subsection 2.3.1 covers material that should be familiar, and, while the 
remaining subsections discuss techniques that are of universal appli 
they lead on to the Gauss divergence theorem, which is the most important 


result of this chapter. 


ation, 


Section 


introduces techniques that will be used throughout the course. 
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2.2 Multivariable calculus 


Functions of two variables 


Most of the techniques that we discuss in this section can be extended to 
functions of three or more variables, but we concentrate on functions of 
two variables for two reasons. First, while the complexity of the results 
generally increases with the number of independent variables. functions of 
two variables often provide a useful insight into the more general cases. 
Secondly, functions of two variables are relatively easy to visualise because 
they can often be represented graphically in a three-dimensional coordinate 
system. 


A real function of two real variables 2 = F (x,y) is a mapping from a set of 
points D in the Cartesian plane to the real numbers. The set D is known 
as the domain of the function, and each point (x,y) of D maps to a unique 
value z. The variables x and y are known as the independent variables, and 
z is the dependent variable. 


Consider the real function F’ defined by F(x,y) = /1—(z2+y2). For Remember that the square 
F(«,y) to be real, we require x? + y? <1, so we take the disc defined by Toot symbol indicates the 
this inequality to be the domain of F. ROSIE aUBED ACIS: 

Here, and generally, we 
assume that the domain of 
a function is the largest set 
on which the definition is 
asible. 


Figure 2.1. The hemisphere 2 = \/1 — (x? + y?) 


If we put z= \/1 — (a? + y?) and let (#,y.2) correspond to a point in a We use either (a, y) or 
three-dimensional Cartesian diagram, then the function F can be repre- (*,y-0) when referring to 


sented as a hemisphere, as shown in Figure 2.1. points that lie in the 
(x, y)-plane. 


Generally, a function w= F(x,y) can be represented as a surface, and 
(x,y. F(x, y)) represents an arbitrary point on the surface. 


2.2.1 Level curves 


It is not always easy to visualise a surface in three dimensions, and there 
is an alternative pictorial representation of a function of two variables that 
is sometimes helpful. On an ordnance survey map you may have seen lines 
connecting points of equal height above sea level. Such lines are known as 
contours, and they give some indication of the shape of the surface. The 
same idea can help us to understand the behaviour of functions. 


If we are given a function z = F(x,y), then we may be able to plot the 
contours or level curves corresponding to the equations F(x.y) = A for 
various values of the constant . 


2.2 Multivariable calculus 


For the function F(x, y) = V1—-(@2 +9) we consider the curves 
V1 (2? +9?) = A, where0<A<1. 

In this case it is convenient to rearrange the equation into the form 
ety=1-), 


and we can then see that the level curves are circles of radius \/1—?, as 
shown in Figure 2.2. Such a diagram can often give some insight into the 
behaviour of an unfamiliar function. 


2 


Figure 2.2. ‘The level curves x? + y 


Exercise 2.1 


Sketch the level curves of the functions F(x, y) = 3 — 2x + y* and G(x, y) = y/x?. 


2.2.2 Continuity 


Perhaps the best way to appreciate the meaning of continuity is to examine a 
function that fails to be continuous; we begin with a function of one variable. 


Example 2.1 


Sketch the graph of the function f(x) = 


Solution 


Figure 2.3 f(x) =1if x >Oand f(x) =—-lifx<0 
The function is defined for all real values of « except x = 0, with 
2 to > O, 
ne)={_4 if <0. 
The graph is shown in Figure 2.3 with a break at the origin. The essential 


point here is that whatever value we assign to f(0), the resulting function 
will be discontinuous at the origin. 
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You saw in Chapter 1 that we define a function f(x) of one variable to 
be continuous at x =a if lim f(x) = f(a). Clearly, in the above example 
za 


lim, f(x) does not exist because f(a) approaches —1 when x approaches 0 
x 
from below, and f(x) approaches 1 when « approaches 0 from above. 


Now we move on to discuss functions of two variables. 


Example 2.2 
© 


Sketch the surface that represents the function F(x, y) = eal 


Solution 


peas 
\zy| 


Figure 2.4 The function F(.x, y) = 


We see that F(x, y) = 1 when « and y have the same sign, and F(«,y) = —1 
when they do not. The corresponding surface is shown in Figure 2.4, and 
it is clear that the function behaves badly on the axes. First, the function 
is not defined on the a- and y-axes, but, even if it were to be defined there, 
the sudden changes in value from —1 to 1 would ensure that the function is 
discontinuous. 


A geometric description of continuity is a useful aid to understanding, but 
it cannot be extended easily to functions of more than two variables. A 
slightly different approach, based on comparing the values of the function 
at neighbouring points in the domain, is more productive. 


For a function of one variable, we define f to be continuous at x =a if 
jim f(x) = f(a). In other words, we can make f(x) as close to f(a) as we 
please by choosing zx sufficiently close to a. It is this notion that we wish to 
extend to functions of two variables. 


We say that a function F is continuous at a point (a,b) if F(x, y) converges 
to F(a,b) as (x,y) approaches (a,b) from any direction. (We shall make this 
statement a little more precise when we discuss functions of many variables 
later.) 


In this course the difficulties that arise from discontinuous functions are not 
our main concern, so you may assume that any given function is continuous 
unless we say otherwise. 


Many of the results that come later in the course are crucially dependent on 
the fact that the relevant functions are continuous, but we shall be avoiding 
these complications. 
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2.2.3. Partial derivatives 


In addition to continuity, we often require that a function should be ‘smooth’, 
which in intuitive terms means that there are no sudden dents or sharp 
corners on the surface that represents the function. This is equivalent to 
the requirement that it is possible to construct a flat surface, known as the 
tangent plane, that touches the original surface at the point in question. 
The existence of such a tangent plane is intimately connected to the existence 
of the partial derivatives of the function, so first we introduce the notion 
of partial derivatives and then we use them to define the tangent plane. 
The partial derivative of a function F(x, y) with respect to x is obtained 
by differentiating the function with respect to « while keeping the variable 
y fixed, and similarly for the partial derivative with respect to y. These 
partial derivatives are often denoted by 

OF a OF 

—— Jee —* 

Ox Oy 
although alternative notations can be useful, e.g. F.(«,y) and Fy(x,y), or 
just F, and Fy. 

OF OF 

The partial derivatives — and by when evaluated at « =a and y = 6, 
represent the slope of the surface z = F(x,y) at the point (a,b, F(a,6)) in 
the x and y directions, respectively. 


Example 2.3 


Given F(x, y) = x? + y?, calculate a and - and hence deduce the slope 
in the x and y directions at the point (3,2). 


Solution 


OF - x 
To calculate a’ we differentiate with respect to x keeping y constant: 
OF _ 
se 


OF ‘ ‘ . 
To calculate a, differentiate with respect to y keeping x constant: 
y 


2z. 


Alternatively, we might have written 2 =? + y°, in which case — 


On 
3y?. 


= 22 


Oy 
Hence the slopes in the a and y directions at the point (3,2) are F,,(3,2) = 6 
and F,(3,2) = 12, respectively. 1 


Exercise 2.2 


Given F(a, y) = x?y? + sin(xe”), calculate GE od oes 
Ox Oy 
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Example 2.4 

Given an equation such as x? + y? — 2? = 1, we need to clarify which are 
the dependent and independent variables. 

(a) Supposing that z is the dependent variable, while x and y are the inde- 


Oz 


pendent variables, evaluate oe and =. 
Ox Oy 


(b) Alternatively, supposing that y is the dependent variable, while x and 
a 
z are independent, evaluate oe and ae 
Solution 


(a) Rearranging the equation, we obtain z = +./x? + y? — 1, which is real 
if a* +y* > 1. Choosing the positive square root, and differentiating 
first with respect to « and then with respect to y, we obtain 


Oz = aa dz = y 


ac ~ V(x? + y*) -1 Oy f(a? +y?)-1 


(with similar results if we choose the negative root). 


Alternatively, we may differentiate the equation x? + y? — 2* = 1 di- 
rectly, First we differentiate partially with respect to .r, and notice that 
the partial derivatives of x? and y* with respect to « are just 2x and 0, 
respectively, while the partial derivative of z* (a function of x and y) is 
Oz : 
care (keeping y fixed and using the chain rule). Thus we obtain 
Oz 
22 — 2z—- =0. 
Ox 


Similarly, differentiating partially with respect to y, we obtain 


Rearranging these equations, we have 
Oz Oz 
On and ay = a where z # 0. 
(b) Differentiating the equation x? + y* — 2? = 1 partially with respect to 
@ and = now gives 


Oy _ eee ae 
2a + Qua, =0 and Wa, —2z=0, 
and rearranging gives 


=—= and a 


On 5 = oD wherey #0. 
Example 2.5 
df [ide ‘a 1 
Show that: a (/ ora) -[ a (ca) dx. 
Solution 
We have 


s(f gee \ 2.8 sl te fie ae 
dt\Jo (+t?) dt| at+t], dt\t 14+t) (+t? #2 


and 


1a a) a= [| = al-a eo 
Lele *<f @+o® |@re|, G+ & 


Note that we use partial 
derivatives (0) to 
differentiate the integrand, 
and total derivatives (d) to 
differentiate the integral. 
This is because the 
integrand is a function of x 
and t, while the integral is 
just a function of t. 


2.2 Multivariable calculus 


Generally, we shall assume that we are justified in ‘differentiating through 
the integral sign’, but if necessary we may appeal to the following result. 


Leibniz’s rule 


Let F(x, t) be a continuous function with a continuous partial derivative 
OF /0t on the rectangular region a < « < b, ty <t < te of the (x. t)- 


plane. Then, for t; < t < fa, 


b 


d 
all. F(z,t)dz = 


This result will be used in Block II. 


Exercise 2.3 


1 1 
Evaluate i sin(ax) da, then use Leibniz’s rule to evaluate xcos(ax) dx. 
0 0 


(It is not necessary to verify the conditions for Leibniz’s rule.) 


2.2.4 Linear approximation 


It is often useful to obtain an approximation to a given function near a 
given point. For a function of one variable, y= f(a) Say, we may use a 
linear approximation f(a) ~ f(a) + f/(a)(a — a), which is equivalent to ap- 
proximating the function by the tangent line at 2 =a. 


For a function of two variables we can obtain a similar linear ¢qtimate in 
terms of the partial derivatives, and in geometric terms this estimasacan be 
interpreted as approximation by the tangent plane. 


The essential idea is that we find a plane that touches the correspondin® , 
surface at a given point, and then we approximate values on the surface by 
values on this nearby plane. To make the ideas more precise, we need the 
equation of a plane. 


yoy) 


Figure 2.5 The surface representing F near the point P 


In Figure 2.5 the surface corresponding to the function z = F(x,y) passes 
through the point P with coordinates (a,b, F(a,b)), and our aim is to esti- 
mate the value of F(x,y) when (x,y) is close to (a,b). 
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The equation of a plane 


The equation z = A+ B(x—a)+C(y—b) represents a plane passing 
through the point (a,b, A). Differentiating partially with respect to x and 
then with respect to y, we obtain 

dz Oz 

= =B = =C. 

Bz and dy Cc. 
In other words, the slope of the plane in the x direction is B, while the slope 
in the y direction is C. 


The tangent plane to the surface z = F(x,y) at the point (a,b, F(a,b)) is 
defined as having the same slope as the surface in the x and y directions, so 
its equation is 


z= F(a,b) +(e— aoe +(y- Oe 


where the partial derivatives are evaluated at (a,b). Thus our linear approx- 
imation for the given function is 


OF 
F(z,y) = F(a,b) + («7 — a) 5 +-9. 


This result is often written in terms of the delta notation for small incre- 
ments, in which case it becomes 
OF OF 
~ ——dar + ——dy, 21 

OF had + Oy iy (2.1) 
where 6F = F(a,y) — F(a,b), 6a =a —a and dy = y—b. Equation (2.1) 
tells us approximately how much the function F changes when we move 
from a point (a,b) to a neighbouring point (a + dz, b+ dy). This approximate 
change in the function is depicted in the inset of Figure 2.5. 


2.2.5 Local extrema 


You will probably recall that the equation (x — a)* + (y—)? = & defines a 
circle of radius 6 with its centre at (a,b). The interior of this circle is defined 
by the inequality (a — a)? + (y— 6)? < &. 


A point (a,b) is said to be a local maximum of a function F(x, y) if there 
is a region (« —a)? + (y—b)? < & such that F(a,y) < F(a,b) at all points 
(x,y) in the region. 


A point (a,b) is said to be a local minimum of a function F(x. y) if there 
is a region (x — a)? + (y —b)? < & such that F(x, y) > F(a,b) at all points 
(x,y) in the region. 


Local maxima and minima are known collectively as local extrema. 


Stationary points 


A point (a,b) is a stationary point of the function F(2,y) if the partial 
derivatives e and a are both zero when evaluated at (a,b); in other 
words, F;,(a,b) = F,(a,b) = 0. 


Stationary points of functions of two variables play much the same role as 
such points for functions of one variable. This stems from the fact that it is 
possible to show that (for a well-behaved function) every extremum must be 
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a stationary point. Thus, if we wish to find a local maximum or minimum 
of a given function, we may restrict our search to its stationary points. On 
the other hand, not all stationary points are local maxima or minima. A 
stationary point that is not an extremum is called a saddle point. 


Figure 2.6 shows examples of typical stationary points. 


wmgy 


local maximum local minimum saddle points 


Figure 2.6 Typical behaviour of a function near a stationary point 


Exercise 2.4 


P) P 
(a) Find the partial derivatives ce and oe of the function 
On” By 
F(a,y) =e" +2? +siny. 
(b) Given /x + /y+ V2 = 1, and regarding 2 as the dependent variable, calculate 
Oz Oz 


(c) Given cos(cyz) = 1, and regarding = as the dependent variable, calculate 
Oz and 22 
= and =. 
Ox Oy 


(d) (i) Show that the surface 
through the point (1, 2 


z= F(x,y) =1+ (x —2y+1)(y — 2x + 1) passes 
—1), and find the tangent plane at that point. 


(ii) Show that F(z, y) has a stationary point at (1,1). Show that the tangent 
plane at the point (1, 1, 1) is parallel to the (, y)-plane. Putting s =x —1 
and t = y—1, or otherwise, find an expression for F(x.y) — F(1.1) in 
terms of powers of « — 1 and y — 1. 


(iii) The lines y —wx and y = x intersect at the point (1,1). What can you 
say about the sign of F(x,y) — F(1,1) on the line y= 2—- and on the 
line y = x? What do you conclude about the stationary point at (1,1)? 


(e) Show that the function F(x. y) = (x? + y*)e” has two stationary points, and 
that the stationary point at the origin is an extremum. 


2.2.6 Higher-order partial derivatives 


From Exercise 2.4, you may have observed that the partial derivatives of a 
function F(«,y) are themselves functions of two variables. This means that 
they too can be differentiated partially with respect to x and y. 
For example, if F(x,y) = 2" siny, then 

OF OF 2 


an =2rsiny and Oy = cosy. 
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Differentiating again, 

aos (5) = 2siny, 42 (=) = 2rcosy, 

Ox \ Ox Oy \ Ox 
2 (F) =2xrcosy and x (=) = —x? sin y. 
These are known as the second-order partial derivatives or second partial 
derivatives of the function F(x,y) and, respectively, they are commonly 
abbreviated to 

PF OF OF OF 

Oat Bydx? Ody “ye? 


or alternatively Fy:, Fyre, Fry and Fyy. 


ia al this is generall, 
Sa os IS 15 general 
Dyor  drdy 2 ¥ 
true provided that the function F(x, y) is sufficiently smooth for the second- 
order partial derivatives to exist and to be continuous. (This result is often 


known as the mixed derivative theorem: we shall not prove it here.) 


It is no accident that in this example we have 


It is, of course, possible to take still higher derivatives, e.g. Frrx, Fyex. 
The mixed derivative theorem then states that the order in which these 
derivatives occur does not matter; so, for example, Fyrxx = Feyrx = Frryx = 
Fraey: 


Taylor approximations 


It is often the case that the linear approximation to a function is not suf- 
ficiently accurate to provide the information that we require: we may need 
a higher-order approximation. Suppose that near (0,0) it is possible to 
approximate the function F(,y) by a quadratic expression so that 

F(x,y) ~ A+ (Bax + Cy) + (Px? + 2Qry + Ry”), 
where A, B, C, P, Q and R are constants. 
This is known as a quadratic approximation for F(x, y) near (0,0). 
More precisely, we suppose that 

F(a,y) = A+ (Bu + Cy) + (Px? + 2Qry + Ry?) 

+ (higher powers of x and y). 

Putting « = y = 0 on both sides of this equation gives A = F(0.0). 
Differentiating both sides of the equation partially with respect to «x gives 

a = B+ (2Px + 2Qy) + (higher-order terms); 
then putting x = y =0 gives B = F,(0,0). Similarly, differentiating par- 
tially with respect to y and putting 2 = y = 0 gives C = F,(0.0). 
Differentiating both sides of the approximation twice with respect to x gives 

as = 2P + (higher-order terms): 
then putting « = y =0 gives P = $F rx (0, 0). Similarly, Q = 3 Fry(0. 0) and 
R= 4F,,(0.0). 
Thus, for a function F(x,y) that is sufficiently smooth near (0,0). the 
quadratic approximation to F(x, y) near (0,0) becomes 

F(x. y) = F(0.0) + F,(0,0)2 + F,(0,0)y 

+4 (For(0,0)x? + 2F,y(0.0)ry + Fyy(0.0)y?) . (2.2) 
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This is known as the Taylor approximation of order two for F(x. y) 
near (0,0). 

(Notice that the linear terms in this expression are identical to the linear 
approximation that we discussed earlier.) 


It is possible to extend this process and obtain an infinite series for F(x. y) 
in ascending powers of x and y. Such a series is known as the Taylor series 
for a function of two variables. 
Example 2.6 
Find the Taylor approximation of order two for F(x,y) = (a + 2y + le” 
near (0,0). 
Solution 
We have 

F, =e + (a+ 2y+1)ye™ and Fy =2e + (a+ 2y+1)zre™. 
Then 

Fre = (2+y + ay + 2y")ye™, 

Foy = (1+ 20 + 4y + ay + cy + Qay?)e™, 

Fy = (44+ a+ 2ry + 2)re™, 
Thus F(0,0) = 1, F,(0,0) = 1, F,(0,0) = 2, Fre(0,0) = 0, Fry(0,0) = 1 and 
Fy,(0,0) = 0. Substituting these values into the general expression (2.2), the 
Taylor approximation is 

(a + 2y + L)e™” ~ 1+ (x + 2y) + $(2ay) 

=1l+at+2y+czy. 

Alternatively, we could obtain the same result more easily by noticing that 
eV = 1+ Hay) + day)? +--+, 80 

(a +2y + Le™ = (w+ 2y +1) (1+ 2y + pay) +) 

=1+x+2y+ zy + (higher powers of « andy). @ 

A similar argument applies to approximations near an arbitrary point (a,b). 


For a function F(x,y) that is sufficiently smooth near (a,b), the Taylor 
approximation of order two for F(ax.y) near (a,b) is 


F(x,y) = F(a,b) + F,(a,b)(a — a) + F(a, b)(y — 6) 
+4 (Frx(a,b)( — a)? + 2Fyy(a,b)(a — a)(y — b) 
+ Fiy(a,b)(y — b)*) « 
(Again, the linear terms in this expression are identical to the linear approx- 
imation to F(«,y) near (a,b).) 
The corresponding infinite series in powers of z — a and y — b is known as 
the Taylor series of F at (a,b). 


We shall not require terms of the Taylor series beyond those corresponding 
to the Taylor approximation of order two. However, if we let dx = x — a, 


by = y —b, po = (a,b), dp = (5, by), and define an operator U = aa + ius. 


then we can write the Taylor series at (a, 6) in the neat form 
= 
1 
F (po + 6p) = F(mo) + > GUE 
k=1 


(where all the partial derivatives are evaluated at po). The first three terms of 
this series give the Taylor approximation of order two for F near (a,b). 
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In this course we shall not 
generally require 
approximations of order 
higher than two. 
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The classification of stationary points 


Suppose that F(x, y) has a stationary point at the origin, so that the first- 
order partial derivatives are zero at (0,0). Then the Taylor approximation 
of order two for F(x,y) near (0,0) becomes 


F(x,y) ~ F(0,0) + $ (Fee(0,0)2? + 2Fry(0, O)ary + Fyy(0,0)y?) - 


Near an extremum at (0,0), it is the sign of F(x, y) — F (0,0) that determines 
if the point is a local maximum or a local minimum, so we are led to consider 
the sign of expressions of the form 


Px? + 2Qry + Ry? 


(where the constants P, Q and R are not all zero). It turns out that the 
sign of this expression is dependent on the sign of PR — Q?. 


Suppose that PR —Q? > 0. Then P cannot be zero, so we may write 


2 _~ Q%)y2 
Pa? + 2Qay + Ry? = Gar ener (R= Oy aaa ew 


where the numerator is positive (except at the origin). 


Thus, at points other than the origin, Px? + 2Qry + Ry? takes the same 
sign as P. 


Now putting P = F,.(0,0), Q = Fr,(0,0) and R= F,,(0,0), we see that 
the function F(a,y) has an extremum at the origin if 


F,(0,0) = F,(0,0)=0 and PR—Q?>0, 


and this extremum is a local maximum or minimum according as P (or R) 
is negative or positive, respectively. 


The above result can be easily extended to a function F(x,y) that has a 
stationary point at a point (a,b) other than the origin. In this case we need 
only observe that the function F(x +a,y +) has a stationary point at the 
origin, and, applying the previous result, we obtain the following. 


A test for extrema 
The function F(x, y) has an extremum at the point (a,b) if 
F,(a,b) = F,(a,b)=0 and PR—Q?>0, 
where P = F,,(a,b), Q = Fyy(a.b) and R = F,,(a,b). 
This extremum is a local maximum or minimum according as P (or R) 
is negative or positive, respectively. 
If PR—Q® <0, it can be shown that F(x.y) does not have an ex- 
tremum at the point (a,b). 
A stationary point that is not an extremum is known as a saddle point. 


Near an extremum at (a,b), F(x,y) — F(a,b) does not change sign, but 
near a saddle point there are always points for which F(x. y) — F(a, 6) 
is positive and points for which it is negative. 


If PR—Q? =0, then the above test for an extremum fails. The point in 
question may or may not be an extremum, and this may be determined by 
examining a higher-order approximation, although we shall not discuss the 
general process here. 
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Example 2.7 
Find and classify the stationary points of the function w = «* + y’ + 3zy. 
Solution 
We have 

Ow 4 Ow 

=- = 32? d — =3;7 ; 

On ac" + Sy ane Oy 3y° + 32. 
so at a stationary point, x? + y = 0 and y* + x = 0. Solving this pair of equa- 
tions, we obtain « = 0 and y = 0, or « = —1 and y = —1, giving stationary 
points at (0,0) and (—1,—1). We also have 

Ow Pw Pw 

= 62, =3 = by, 

pat" Deby SEE Cape ee 

so 


PR-Q= Pw Pw (= 


Ox? Oy? 

At (0,0) we have PR — Q? = —9 < 0, so the stationary point is a saddle 

point.-At (—1,—1) we have PR — Q? = 27 > 0, so the stationary point is an 
Pay 

extremum. Also at this point we have P = a = 6r evaluated at (—1,—1), 


giving P = —6 < 0, so the extremum is a local maximum. @ 
In the functions in the following example, PR — Q? = 0, so the test fails. 


However, in these cases, we can use elementary analysis to determine the 
nature of the stationary point. 


Example 2.8 


(a) Given F(a,y) = 24y!, show that PR —Q? =0 and that the function 
has a local minimum at (0,0). 


(b) Given F(a,y) = xy, show that PR —Q? = 0 and that the function 
has a saddle point at (0,0). 


Solution 
(a) Since F,(0,0) = F,(0,0) = 0, (0,0) is a stationary point. We have 
Ox? Oy* Oxdy 
so at (0,0) we have PR—Q?=0. However, F (x.y) = aly! >0= 
F(0,0), so there is a local minimum at (0,0). 
Since F,(0,0) = F,(0,0) = 0, (0,0) is a stationary point. We have 
Ox? Oy? Oxdy 
so at (0,0) we have PR — Q? = 0. However, F(x,y) — F(0.0) = 2%y', 
which is positive when a and y have the same sign, and negative when 


they have opposite signs. Thus there cannot be a local extremum at the 
origin, and the point is therefore a saddle point. Mf 


2 
) = (1222y4)(1224y2) — (1629y2)?, 


= 


2 
) = (6ry$)(625y) — (92°y?)2, 


Exercise 2.5 


(a) Find the Taylor approximation of order two to the function w = e*Y near (1,2). 
(b) Find and classify the stationary points of the function w = sin cosh y. 
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Later we shall discuss partial differential equations, and this is an opportune 
moment to introduce a solution of such an equation. 


Exercise 2.6 
Show that any function of the form F(x,y) = f(x) + g(y), where f and g are 


differentiable, satisfies the equation 


catia 
Foy 


2.2.7 Functions of many variables 


Suppose that we wish to investigate the behaviour of a function F' of the n 
variables 2, 2r2,...,2,. It is often convenient to write @ = (21,22,...,2n) 
and to regard F as a function of the vector a. 


Let @ = (@,2,...,2,) and @ = (a), a2, mn). Then we define a func- 
tion F(a) to be continuous at a@ if lim F(x) = F(a), and note that 2 may 
za 


approach a from any direction. 


The linear approximation is easily extended to functions of many variables. 


If a = (a), a2,..., a), we write 6F = F(x) — F(a) and dx; = 2; — a; for 
i=1,2,...,n. We then have 
OF OF OF 
6F ~ —6: 6. s+ ——6tn, 
OF On ry Saas ng + ee in 
where all the partial derivatives are evaluated at a. 


It is also possible to extend Taylor approximations of higher order to fune- 
tions of many variables, although, as one might expect, the complexity of 
the expressions increases with the number of variables and with the order 
of the approximation, 


Exercise 2.7 
Given w = (1r),a72,2r3) and 


F(x) = A+ (Bix, + Borg + Bgx3) 
+ (Criat + Coor} + Cygx3 + Corie + Cogr2r3 + Cyyr97}), 


calculate the coefficients in this expression in terms of partial derivatives of F’. 


Taylor approximation of order two 


For an arbitrary function F(a) of three variables, the Taylor approximation 
of order two near a = a is 


F(a) > F(a) + (x1 — a1) + (x2- a) + (23 — a) 
+5 (Cera + (ea — an) 8 + on 09) 
+ 2(x) — ay) (x2 — on) gor 
+ 2(a2 —a2)(z3 — 03) go 


+ 2(a3 —a3)(a1 — apna) 4 


Recall that a partial 
differential equation is any 
equation involving partial 
derivatives. 


We shall discuss vectors in 
the second part of this 
chapter. 


Exercise 2.7 illustrates 
such an approximation 
when a= 0. 
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where all the partial derivatives are evaluated at a =a. It is easy to see 
how this result can be extended to more than three variables. 


Stationary points 


The notions of stationary points and extrema can be easily extended to 
functions of more than two variables. 


Near the point @ = (a;,@2,....@n), it is the nature of the difference 
6F = F(a) — F(a) that characterises the behaviour of the function F. If 
this difference is positive (or zero) for all a close to a, then the function will 
have a local minimum at a; alternatively, if the difference is negative (or 
zero), then the function will have a local maximum at a. However, if in every 
neighbourhood of a there are points at which 6F = F(x) — F(a) is positive 
and points at which it is negative, then there cannot be an extremum at a. 


There is an alternative approach to stationary points which will be relevant 
in later chapters. The idea is to examine the behaviour of a function as we 
approach a stationary point from a nearby point. In Block III we shall see 
that this alternative approach can be extended to a situation in which there 
is no direct analogue of the notion of a partial derivative. 


Since we are concerned only with points close to a, it is convenient to write 
z=a+e€, where € = (€.6,-.-.€,), with @+G+---+@ =1, and e 
is small. (Essentially, this means that we approach the point @ along the 
straight line joining a to a+.) We expect dF to be small when is small. 
We might expect that 6F = F(a +e) — F(a) is approximately proportional 
to ©, or, more specifically, that 
F(a + ¢€) — F(a) 
€ 

is bounded (i.e. not infinite) when ¢ is small. For well-behaved functions, 
and for an arbitrary point a, this is generally the case, and we have 


OF OF OF 
Flateé)— Fla) =e (Se + eat + 5E&) = 000) 


Example 2.9 
Given F(x) = xyx2 + 273 + 232) and a = (1, —1,2), find F(a + <€) — F(a) 
and show that F(a+<€) — F(a) = O(¢) for all €. 
Solution 
We have 
F(a +e€) = (1 + €))(—1 + €2) + (—1 + €€))(2 + €€3) 
+ (2+ e€5)(1 + €€)) 
= —1+e(-€ + & — &3 + 2& + €3 +28) 
+ £7 (Ei£> + fobs + £361) 
= —1 + £(& + 3Ep) + €7(Exb> + bobs + E361) 
and F(a) = —1. so 
F(a + 2€) — F(a) =e(€; + 36) + 7°(8:& + £28 + 6361) = Ol). 
The previous example illustrates what we expect to happen for an arbitrarily 


chosen point a. If < is small, then <? is smaller, so the difference F(a + c€) — 
F(a) is approximately proportional to <. 


7I 
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However, for a given function F’, there may be exceptional points a for which 
F(a+c€) — F(a) =O(e2) for all € 
(in other words, 6F is smaller than we might expect). These are, of course, 
the points for which 
Or On. Ore 
6m 822 Om * 
so that the linear approximation for 6F is zero. These exceptional points 
are in fact the stationary points of F. 


Example 2.10 


Given F(a) = 222 + ror3 + xgx; and a = (0,0,0), find F(a + e€) — F(a) 
and show that F(a +<€) — F(a) = O(e?), so F has a stationary point at 
the origin. 


Solution 
In this case we have 
F(a +e€) — F(a) = £7(E,£ + &€3 + &36:) = O(€*), 
and we conclude that the point a = (0,0,0) is a stationary point. 


We could, of course, have reached the same conclusion by calculating the 
partial derivatives 


OF _ OF _ OF 
Ox, Oxo ~ arg 
is a stationary point. 


We can see that = 0 at (0,0,0), which confirms that it 


We could, therefore, have defined the stationary points of a function F to 
be those points a for which F(a +¢&) — F(a) = O(<*), and indeed it is this 
formulation that is most easily extended to the cases we shall consider later 
in the course. 


Generally, the classification of stationary points (into maxima, minima and 
saddle points) involves the determination of the sign of a quadratic expres- 
sion in many variables, but we shall not discuss that process for functions 
of more than two variables. 


2.2.8 The chain rule 


Often we encounter situations where we wish to transfer between coordinate 
systems. For example, we may find it convenient to transfer from Cartesian 
coordinates (x,y) to polar coordinates (r,@), in which case we may wish 
to represent partial derivatives with respect to x and y in terms of partial 
derivatives with respect to r and @ (or vice versa). The device which makes 
this possible is known as the chain rule. 


Suppose that we are given a variable w which is dependent on two variables, 
a and y. Then the linear approximation gives 
bw ~ ox + Sy, (2.3) 
ae 


where the accuracy of the approximation increases as dx and dy both ap- 
proach zero. But suppose now that the variables « and y are themselves 
dependent on a parameter, ¢ say. For example, we could be given w = e™ 
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where x = cost and y = sint, and wish to examine the behaviour of w as 
t varies. (In this case the point (a,y) would lie on a circle, so we would 
examine the behaviour of w = e*” for points on this circle.) 


Suppose that dt is the increment in ¢ that gives rise to the increments dx 
and dy. Then (dividing both sides of equation (2.3) by dt) 

dw | Owda | Ow dy 

ot ~~ Ox dt Oy dt” 
and in the limit as 6¢ approaches zero, the approximation becomes equality 
and we have 

dw Owda 4 Ow dy 

dt = dx dt * dy dt 
This result is known as the chain rule for functions of two variables. 
(It is usually abbreviated to the chain rule when there is no danger of 
confusion.) 


Suppose now that «2 and y are dependent on two variables, u and v say 
(rather than a single variable t). Then the previous argument applies with 
one minor difference. Keeping v fixed, the derivatives in the chain rule now 
become partial derivatives, and we have 


Ow  Owdxr . Ow dy 
Ou Ox du” dy du" 
Similarly, if we keep w fixed, then we have 
Ow Owdr Ow dy 
Ov” Ox av * By du’ 


These are the relationships that enable us to change from one coordinate 
system to another. In Section 2.4, we demonstrate how to switch between 
derivatives involving Cartesian and polar coordinates. 


Exercise 2.8 


dw 
(a) Given w = e*, where x = sint and y = cost, calculate a first by substituting 
for x and y, and then by using the chain rule. 
(b) Given a function w = w(«,y), we make the substitutions «= au + bv and 
y = cu + dv (where a, b, c and d are constants) to express w as a function of 
u and v. Express the partial derivatives 
dw dw Ow Pw Pw 
Ou’ Av’ Au?" Audv Ov? 
in terms of partial derivatives with respect to x and y. 


and 


2.2.9 Multiple integrals 


In Chapter 1, we introduced the idea of definite integration in a single vari- 
able as the limit of a sum. In this section we generalise this idea to integra- 
tion over two or more variables. We begin with the case of integration over 
two variables. 


Integrals over plane surfaces 


We consider a function f(r) = f(x,y) defined on a region S of the Carte- 
sian plane. We first divide the region into a large number, N, of small 
elements 5A,,. as shown in Figure 2.7(a), so that the sum of these elements 
approximates the region S. 
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‘A pair of functions 
(a(t). y(t)) defines a 
parametric curve. 


Note that we use partial 
derivatives (0) when 
differentiating w with 
respect to x and y, but 
total derivatives (d) when 
differentiating « and y 
with respect to t. This is 
because w(z, y) is a 
function of two variables, 
but «(t) and y(t) are 
functions of a single 
variable. 


These small elements are 
not necessarily rectangular, 
although it is convenient to 
define them by a grid as in 
Figure 2.7. 
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(a) (b) 
Figure 2.7 Regions divided into rectangular subregions 


For each such region we choose a vector r,, corresponding to a point lying 
in JA,,. Then we define the integral of f over S to be the limit of the sum 


N 
DY frm) 64m 


m=1 


as N tends to infinity (and the size of the elements shrinks), so that 


N 
: L(r)dA = lim > f(tm) 6Am- 
s N—00 
‘ m=1 
It is not essential to form the small regions using a Cartesian grid as in 
Figure 2.7, although this may enable us to write the integral in a more 
recognisable form. 


The simplest case is when the region S is a rectangle, as shown in Fig- 
ure 2.7(b). In this case we have JA; = 62; dyj, and the sum over m can 
be found by first keeping j fixed and summing along a horizontal strip (i.e. 
over i), then summing over j. This process corresponds to first integrating 
{(a.y) with respect to 2 while keeping y fixed, and then integrating the 
result with respect to y. Thus the integral becomes 


[saa = i ([ sae) dy 


(summing first across an arbitrary horizontal strip, and then vertically). 
Alternatively, we could reverse the order of summation to obtain the same Changing the order of 


value: summation can produce a 
b d different value, but only in 
1A = r, y) dy ) d: very exceptional cases, 
[soe | (/ Fu) v) <t which will not be 


5 ‘ * r encountered in this course. 
(summing first up an arbitrary vertical strip, and then horizontally). These 


forms are known as repeated integrals. 


In our first example we shall consider a particularly simple case, namely 
where f(x,y) = 1. 


Example 2.11 
Ly pe 
Evaluate the repeated integral [ ( fi ay) dx. 
0 Jat 


Solution 
We have 
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Note that in later chapters we use the notation / dx f(x) rather than 


[re dx, because writing fe fate y) in place ot [ ( [ste y) ay) dx 
avoids the need for brackets. 


Example 2.12 
Evaluate the integral | f(r)dA where f(r) = f(x.y) = 27y and S is the 


s 
rectangle with its sides parallel to the axes for which one vertex is at the 
origin and the diagonally opposite vertex is at the point (1,2), as shown in 
Figure 2.8. 

Solution 
We have 


A slight adaptation of the method allows us to integrate over more compli- 
cated regions, as in the following example. 


Example 2.13 
Evaluate the integral [teaa where f(r) = f(x,y) = x7y and S is the 


8 
region bounded by the curves y = x? and y = 1— 2, for « > 0, and the 
y-axis (so S is defined by the inequalities y > x, y < 1 — a? and x > 0; see 
Figure 2,9). 


Solution 


Figure 2.9 An elemental ‘rectangular’ area in the region S 


Writing the integral as a repeated integral, where we integrate first along a 
strip from one curve to the other, as illustrated in Figure 2.9, we have 


1/2 1—z? 
[saa a? (/ vi) 
Js 0 2 
1/v2 
-[° ? vs 
0 
1/ v2 one 
if 2*((1— 27)? — 24) de 
0 


o 
im * a2(1— 222) de = V3/60. . 


0 


2 


dx 


U1 
tole 
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Figure 2.8 


It is often useful to draw a. 
diagram of the integration 
region, as this helps in the 
construction of the 
integration limits. 
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You may wish to repeat this calculation with the order of integration re- 
versed, i.e. perform the x integral first, and then the y integral. However, be 
warned that it will be necessary to alter the boundary conditions. It may 
be instructive to first examine the solution to Exercise 2.9, which solves a 
similar double integral in two ways. 


Exercise 2.9 


Evaluate the integral if a? dA where S is the region in the first quadrant bounded 


by the hyperbola xy = 16 and the lines y = 7, y =O and 2 =8 (soy <2,0< <8, 
y >Oand y < 16/2). 


Integrals involving polar coordinates 


Recall that a point (x,y) in the plane can also be labelled by polar coordi- See page 24 of Chapter 1. 
nates (7,0), where 2 = rcos@ and y =rsin0. 


The preceding method can be extended to integrals over areas defined in 
terms of polar coordinates. 


Example 2.14 


Figure 2.10 The cardioid r = 2(1 — cos@) 


The curve shown in Figure 2.10 is known as a cardioid, and is defined in terms 
of polar coordinates by the equation r = 2(1 — cos@) where —z < 6 < 7. 


Find the area of the region bounded by the curve. 


Solution 


(b) 


Figure 2.11 The region bounded by the cardioid divided into elements that are 
approximately rectangular 


In order to find the area S bounded by the curve, we divide the region into 
small elements as shown in Figure 2.11(a). The element shown is approx- 
imately rectangular, so its area is approximately r4@dr. Summing over all 
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such elements produces an integral of the form "fi if r dr d@, and writing this 
as a repeated integral we have 2 


7 2(1—cos @) 
[frerae- | rdr ) d0 
on 


-2/" Er] 2] Leas A) 


=f 4(1—cos0)*d0=67. 


Jo 


The preceding example illustrates a general rule which allows us to transform 
between integrals of functions using Cartesian and polar coordinates. For 
any function f and integration region S, 


[ff semaeay = ff r5¢rcos0,rsin®) doar, 


where the region S$ and the function f are expressed on the left-hand side 
in Cartesian coordinates, and on the right-hand side in polar coordinates. 
Transforming between coordinates is often essential for performing an inte- 
gral, as the following example illustrates. 


Example 2.15 


Show that for a > 0, the integral of the Gaussian function satisfies 
This result was quoted in 


a 
a oe ere ey kis Chapter 1 as 
a) = [exo eM vz equation (1.4) 


Solution 


We first make the change of variable u = Vax. so 


o 
I(a) = =f exp(—u’) du = 4a T(1). 
90 


So it is necessary to evaluate /(1). We do this by considering its square: 


wo (rere) foes 
=f (fee ae) av 


We transform this integral to polar coordinates (r,@), remembering that 
r? = x? +y? and that the region of integration is the whole two-dimensional 
plane 0 < 0 < 27,0<r<oo: 


2 0 an a oC . 
(1(1)) =| (f ee rad) dr= a [ e" rdr 
0 0 0 


=—1 [en all Rss 
Hence we have shown that I(1) = 7, so I(a) = VE as required. Wf 
Exercise 2.10 


Evaluate the integral [ r? dr dO, where S is the region defined by 0 <r <1 and 
Ss 
0<0<7/2. 
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Triple integrals and multiple integrals in general 


The notion of a double integral can be generalised easily to integrals of 

functions of three variables. We consider a function f(r) = f(x, y,z) where We assume that you have 
r = xi+yj+ zk is defined on a region R of three-dimensional space. We met vectors before; we 
then divide the region R into a large number, N say, of small regions, each shall refresh your memory 
of volume AVj, and the triple integral is defined as a limit of a sum: a 


[ff tensa =, 


Such integrals are often known as volume integrals. 


N 
im SO fis de 24) AVi- 
=I 


In this course multiple integrals generally arise from theoretical considera- 
tions, and we rarely need to evaluate them except in cases where the sym- 
metry, or other factors, makes the calculation particularly easy. However, 
as an aid to understanding, it may help to see how we might undertake such 
a calculation. For this purpose we shall choose a very simple region. 


Example 2.16 


Figure 2.12 The rectangular box R divided into elemental subvolumes 


In this example we suppose that R is the rectangular box defined by the 
inequalities 

O<a<a, O<y<sbh O<zK<ce 
(as shown in Figure 2.12), and that 

S (2,2) = (2? + ye. 
We then have 


ri om y,z)dV = he (f ( [a +y)e* 7) ay) dx 
= i (fw +y)0- ay) dx 


a b 
=(1-e~) ( (2? + vay) dx 
0 0 
=(1- @9f [xy + 1p dx 
=(1- ef (bx? + 3b3) der 
0 


= 4(1—e™*)ab(a? +6"). 


Integrals over surfaces that are not planar 


We define integrals over surfaces that are not planar in a similar fashion, 
first dividing the surface into a large number of small elements, then forming 
a sum, and finally taking the limit. In practice we try to choose a coordinate 
system which allows us to write the integral as a repeated integral, as in the 
following example. 
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Example 2.17 
Given f(r) = 32 — \/x? + y?, evaluate the integral [seas where S$ is 
s 


the curved surface of a cylinder of height h/2 and radius R with its axis on 
the z-axis and its base on the plane z = h/2. 


Solution 


Figure 2.13 The curved surface of a cylinder divided into elemental rectangles 


In this case it is helpful to choose cylindrical polar coordinates (r, 0, z) where 
a=rcos@ and y=rsin@. Then f(r) = 32 — /2?+y? =32—r, and an 
area element on the surface is 6A = R 6062 (see Figure 2.13). 


On the curved surface of the cylinder we have r = R, so 


h ‘Qn 
[teaa= Of (a: — R)Rae) dz 
s n/2 \Jo 


he Qn 
= [ (f 7) (32R — R®) dz 
Snj2 \Jo 


32R — R*) dz 


= hx(9h—4R)R. @ 
Most of the integrals of this kind that will arise in this course are particularly 


simple. For example, if f(r) = 1 and S is a sphere of radius p with its centre 
at the origin, then 


[soaa 
Ss 
Similarly, if f(7) 


I firydd= (pet 2) [ dA = 4mp?(p +2). 
s s 


dA = 4rp". 


r|* +2, then 


End-of-section Exercises 
Exercise 2.11 


Find the partial derivatives 


Pw Pw. 
Sand <2 if w= — ’ 
Oe Oe Vee 


Exercise 2.12 


Verify the mixed derivative theorem (see page 66) for the function 
fd 


Fev) = ae 


719 


Recall that a sphere of 
radius p has surface area 
4np?. 
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Exercise 2.13 


Show that each of the following functions is a solution of 
This equation will be 
discussed briefly in 


Section 2.4. 

the two-dimensional form of Laplace’s equation in Cartesian coordinates. 
(a) w=e*cosy — (b) w Bary? 
Exercise 2.14 
Find and classify the stationary points of the function 

F(a,y) = 1 — 22x — 12y + 10x? + 2ary + 5y?. 
Exercise 2.15 
Find the point on the plane 2x — y + 2z = 16 that is nearest the origin. 
Exercise 2.16 
Given w = a? + y? + 2?, calculate 

Pw ow ow 

oe aes ds SS 

Ordydz' dz0xdy Dyd=dx 
Exercise 2.17 

=f pcosd 
Evaluate the repeated integral [ ( i psn) dé. 
Jo \Jo 

Exercise 2.18 This exercise is optional. 


The region S in the (a, y)-plane is bounded by the three lines « = 3, y = x/3 and 
y =0. Write the integral of the function exp(2?) over the region S$ as a repeated 
integral in two ways, namely with the «x integral performed first, and then with the 


y integral performed first. Use one of these to evaluate exp(a*) dA. 
s 


2.3 Vector calculus 


In this section we shall be concerned mainly with the behaviour of func- 
tions whose domains are sets of vectors, so we now revise some elementary 
properties of vectors. 


2.3.1 Vector algebra 


We assume that you are familiar with the fundamental notions of vector 
algebra. This is essentially the material that you are likely to have met in 
previous Open University courses, but for students who have gained their 
knowledge of the subject from a different source, we provide a very brief 
summary and explain our notation. 


A brief summary 


Scalars are quantities that are defined by their magnitude, such as mass 
and length, and we use the term to distinguish them from vector quantities, 
which have both magnitude and direction, such as velocity. Vectors are often 
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denoted in printed material by bold type. for example v or F, or by under- 
lining in handwritten material, v or F (although several other conventions 
are in common use). They are often specified by pairs of numbers in two 
dimensions, and triples of numbers in three dimensions: for example, (x,y) 
in two-dimensional Cartesian coordinates, or (x,y,z) in three dimensions. 


The following statements are made in the context of three-dimensional vec- 
tors, but similar results apply in two dimensions. The vector (x,y) in two 
dimensions is often equated with the vector (x.y,0) in three dimensions. 


A vector r = (x,y, 2) may be multiplied by a scalar a to give 

ar = (ax, ay. az), 
and two vectors in Cartesian form may be added to give 

ry +12 = (1,491.21) + (@2, yo. 22) = (1 + 22,41 + Y2, 21 + 22)- 
The zero vector 0 corresponds to (0,0.0) in Cartesian coordinates. 
The magnitude or length of a vector r = (x,y,z) is denoted by |r| and 
defined to be |r| = /2? +y*+2?. A unit vector is a vector of magni- 
tude 1. The unit vectors in the directions of the Cartesian coordinate axes 
are labelled i, j and k, so a three-dimensional Cartesian vector r = (x,y, 2) 


may also be written in the form r = xi+yj+ zk. The scalars x, y and z 
are known as the Cartesian components of the vector r. 


We write a ‘hat’ symbol over a vector to denote the fact that it is a unit 
vector. Thus a unit vector in the direction of a given vector r is denoted by 
= Sug 

T, where F = —. 


{rl 
The scalar product 


The scalar (or dot) product of two vectors a and b is 

a-b= |a||b| cos 0, (2.4) 
where @ is the angle between the vectors. Using a unit vector 6, the scalar 
product a- 6 represents the length L shown in Figure 2.14; this is known as 


the projection of the vector a in the direction of b. Note that the scalar 
product of two vectors is a sealar. 


Figure 2.14 The projection of a vector a in the direction of 6 
In Cartesian components, the scalar product of vectors a = ayi + aaj + agk 
and b = byi + boj + bgk can be written as 

a-b=ayb; + agb2 + agb3. (2.5) 


Note that it is not immediately obvious that equations (2.4) and (2.5) are 
equivalent. However, you should be prepared to use either definition as the 
need arises. 


The following list includes some useful properties of the scalar product. 


Strictly speaking, (x,y, =) 
is a representation of the 
vector r rather than the 
vector itself, but we shall 
not be concerned with such 
fine distinctions. 
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Some properties of vectors 


i-i=j-j=k-k=1, i-j=j-k= 
a-b= |a||b|cos@ 
= (a,i + aaj + agk) - (bi + boj + bak) 
= ab; + agbe + agb3 
=b-a, 
lal = aj +a} +43 = Va-a, 
\Aa| = |A||a| for any scalar A, 
(a+b)-c=a-c+b-ec, 
(A+ p)a = Aa + pa for any scalars \ and pr. 


Example 2.18 


For a fixed vector v, and any unit vector n, show that v-7 has maximum 
value |v|, when 7 is chosen to be parallel to v. 


Solution 
Using equation (2.4), 
vi = |v| [7] cos 0 = |v| cosd. 


For any given v, this has a maximum value |v| when cos@ = 1. This corre- 
sponds to @ = 0, i.e. to choosing % to be parallel tov. 


Exercise 2.19 


(a) Find the cosine of the angle @ between the vectors a = (1,2.—1) and 
b=(2,1,1). 


(b) Find the magnitude of the vector ¢ = 3i + 2j — 5k. 


2.3.2 Scalar fields 


In applied mathematics, real functions of two and three real variables are 
often known as scalar fields. The term ‘scalar’ is used to distinguish them 
from the ‘vector fields’ that we shall discuss shortly. We begin with a physical 
situation which may help you to understand the general cases that come 
later. 


If a drop of ink falls onto a large sheet of wet blotting paper. then the stain 
will spread. The density F of the stain (i.e. the mass of ink per unit area) 
at a particular instant of time might be represented quite well by a function 


F(x,y) = Aexp[-a(x? + y’)] . 


where A and a are positive constants (and the blotting paper is assumed to 
lie in the (x, y)-plane with the centre of the stain at the origin). In this case 
F(x,y) is an example of a two-dimensional scalar field. 


Suppose that time is kept fixed, and imagine that we follow a curve drawn 
on the paper through the ink stain: we shall notice a change in the ink 
density as we traverse the path. At any particular point on the path, we 
intend to determine the direction in which the ink density increases most 
rapidly. 
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If we represent F(,y) as a surface, as in Figure 2.15(b), then the height of 
the surface above the point (a. y) represents the density of the stain at that 
point. 


(a) 


Figure 2.15 (a) The ink stain in the (x, y)-plane; (b) the surface representing 
F(x. y) = Aexp[—a(x* + y")] 


The gradient 


At most points on the surface representing F we should be able to deter- 
mine both the direction in which the slope is greatest and the magnitude of 
this greatest slope. We can then associate a vector having this magnitude 
and direction with each point (x,y). Since the magnitude and direction of 
greatest slope will vary from point to point, the magnitude and direction 
of the vector representing this will also vary from point to point. In other 
words, this vector will be a function of position (x,y). This vector func- 
tion is known as the gradient of F', denoted by VF; we sometimes write Read V as ‘grad’. 


VF (x,y) to emphasise its positional dependence. Some authors write grad /” 
i instead of VF. 
For our ink stain illustration, VF always points in the direction for which 


the ink density increases most rapidly. In Figure 2.15 we show VF at two 

points, (x,y) and (2’,y’), by the darker arrows. If we were to draw the level The lighter arrows point 
curves for the surface representing F, then we would obtain a diagram very uP the lines of greatest 
similar to Figure 2.2 (page 59), and we would see that the gradient at each slope, and the darker 


st i et SON arrows are their projections 
Hat is pepeeealss to the tangent to the level curve that passes through Onto the (x, y)-plane. 
hat point. 


We now formally define the gradient of a function, and state some properties 
associated with it. 


The gradient 


The gradient of any function F(x, y) is defined as 


OF, OF 
F = —i+—j. 
x Ox ue Oy 2 
At each point (x,y), VF is a vector which points in the direction of 
maximum slope (uphill) of F; its magnitude |VF| is the value of the 
maximum slope at that point. 


Additionally, if 7% is any unit vector, then VF'- 7 gives the slope of F 
in the direction of 7. 


The proof of these properties is presented as an optional example. 
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Example 2.19 This example is optional. 
Prove the aforementioned properties of VF. 
Solution 


We first need to find the slope of the surface representing F in an arbitrary 
direction, say along the direction of the unit vector m = cos @i + sin @j. 


We consider the change in the value of F between two adjacent points on the 
surface, P and Q say, with coordinates (x,y,z) and (x + dx, y + dy,2 +42), 
respectively, as shown in Figure 2.16. 


Figure 2.16 The points P and Q on the surface representing F 
We now keep P fixed and consider the slope of the line PQ. If we let 
6s = \/6x? + dy?, then we see from Figure 2.16 that the slope of the line 
joining P and Q is given by Ha We know (from equation (2.1)) that 

bz > F,dx + F, dy, 
where the partial derivatives are evaluated at P. Therefore 


ox oy 
at ry, 


This result holds in the limit as ds — 0, so 


slope of PQ = "= ~ Fy = F, cos@ + F,sind. 


slope in the direction of A = (F,i+ F,j)-ma= VF +n. 


So, keeping P fixed so that F, and F, are constant, we can calculate the slope Of course, this argument 
in any direction simply by varying the direction of m. From Example 2.18, _ fails if F, = F = 0, since 
we see that the maximum slope has value |VF|, and this occurs when VF then there is no direction 
is parallel to”. of maximum slope. 


We mentioned earlier that at every point VF is perpendicular to the tan- 
gents of the level curves of F, but we shall not prove this here. 


In the case of the ink stain example, we have F(x, y) = A exp[—a(zx? + v)], 


so 
tae = —2a Ax exp[—a(«? + y)) 
Ox 
and 
OF 
oy = —2aAy exp[—a(x? + y)| o 
therefore 


VF = —2aAexp[—a(zx? + y”)] (xi + yj). 
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This is a vector in the same direction as a vector from the point (x,y) to 
the origin, which is what we would expect from the symmetry shown in 
Figure 2.15. 


Similar ideas can also be applied to a three-dimensional scalar field F(x, y. 2). 
for which we define 

OF OF, OF. 

VF =—i+ —i+;-k. 

ioe Bye Oz i 
It is sometimes useful to regard V = ig +45, + wo as an operator that 
can be applied to functions, in which case it is known as a vector differ- 
ential operator. 


Example 2.20 
Calculate VF for the function F(x, y,z) = x?ysin z. 
Solution 
We have 

OF Soar eiis OF Btn» dd Ls) ele ee 

De = 2eysinz, i sinz and 53> =2'ycosz, 
3 OF, , OF. , OF 

VF= Ont + Dy! + a 

= (2ay sin z)i + (x sin z)j + (x?ycosz)k. 

Example 2.21 


Calculate the gradient of the function F(a, y) = a? — xy + 3y? at the point 
(2;—1). 


Solution 


The point (2, —1) refers to a point in the domain of F (i.e. the point (2,—1,9) 
on the corresponding surface). F,(x,y) = 2a —y and Fy(x.y) = —x + 6y, 
so F,(2,—1) =5 and F,(2,—1) = —8. It follows that VF (2, —1) = 5i — 8). 


A schematic representation of VF, together with the level curves of F, is 
shown in Figure 2.17. The arrow at each point represents the vector VF at 
that point. Note that at each point, VF is perpendicular to the tangent to 
the level curve of F. 


= 
- 
-_ 
-_ 


Figure 2.17 The level curves and gradient of F(x.y) =2*—ay+3y? 


V is known as del or 
nabla in some texts. 
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Exercise 2.20 


Calculate the gradient of the function F(x.y) = (x? + 2y”)sin|(2 + y)z] at the 
point (1,1). 


2.3.3 Vector fields 


A vector-valued function F(x, y) = Fy(x.y)i+ Fo(x.y)j in two dimensions, 
or F(x, y,2) = Fy(x,y.z)i+ Fo(a.y, 2)j + F3(x.y, z)k in three dimensions, 
is often known as a vector field, and we commonly make abbreviations like 
F = F\i+ Faj when the choice of independent variables is clear. 


We saw in the previous subsection that the gradient of a function is a vec- 
tor field; but not all vector fields arise as gradients. Physical examples of 
vector fields include fluid velocity (where the vector at each point is equal 
in magnitude and direction to the fluid velocity), and the force on a charged 
particle in an electric field (where the vector at each point is equal to the 
force on the particle at that point). 


The spatial variation of a vector field is inevitably more complicated than 
that of a scalar field: for example. in three dimensions, it involves the nine 
partial derivatives 

OF, OF, OF, OF, OF, OF, OF3 OF; OF 

Oz’ Oy’ O2'’ Ox’ Oy’ dz" O2" Oy' dz- 
These partial derivatives completely describe the spatial variation of F, 
and certain combinations of them occur quite often. The combination that 
interests us here corresponds to the scalar field known as the divergence of 
F, defined as the scalar product of the operator V with the vector field F; 


Oe ein e. ins Some authors use the 
5 (a +5,)+ ak) + (Fil + Faj + Fak) notation div F in place of 
5 a VF. 
— 2H Om, OF 
Be By: Oa, 


Note that although F is a vector field, V - F is a scalar field. The importance 
of V - F will become clearer when we discuss the Gauss divergence theorem; 
however, this theorem involves integration over surfaces and volumes, so first 
we need to say a little more about integration. 


Exercise 2.21 


(a) Given F(x, y, 2) = wi+ yj +k), calculate V - F. 


ya ( 


(b) Given f(x, y,z) =x? + y? + 27, calculate V+ (VS). 


2.3.4 Integration 


We now return to the topic of integration. We begin with the problem of 
how to determine the length of a curve, before moving on to discuss the 
integration of vector-valued functions. 


2.3 Vector calculus 


The length of a curve 


We continue with the idea of integration being the limit of a sum. 


o+- 
ay 


' 
' 
' 
[a 

Figure 2.18 A curve divided into short subintervals 
Suppose that we wish to determine the length L of the curve y = f(x) 
between « =a and x = 6, as shown in Figure 2.18. One approach might be 
to divide the curve into a large number (say \) of small segments that are 
almost linear. The sum of the approximate lengths of these small segments 
would then give us an approximation to the total length of the curve. If 
we increase the number of such subdivisions, we might expect to obtain an 
increasingly accurate approximation to the total length. In other words, the 


length of the curve could be represented as the limit of a sum, i.e. as an 
integral. 


For the length 6s, of the kth segment, we have ds? ~ dx? + dy?. So 


N N N 2 
L~ > bax = > dx} + OZ = DY dey + (**) ; 
k=1 k=l k=1 


and in the limit as N — oo we have 


v= f+) « 


Example 2.22 


Calculate the length of the graph of the function y = x°/? from the origin 
to the point (1,1). 


Solution 
The required length is 


1 dy 2 1 Tame 
L= f y+ (2) de = [ yil+ (32!?) de 


1 
-1/ V4+ 9a dx 
0 


1 
oe 3/2 
[4 +92) iF 
+ (13V13—8). 


i] 


The same idea can be extended to a curve that is not the graph of a function 
(ie. not of the form y = f(x)), for example the curve shown in Figure 2.19. 
Such curves often arise in the parametrised form xr = x(t), y = y(t) with 
ty, <t < ta, and are called parametric curves. A smooth curve is defined to 
be a curve for which the functions x(t) and y(t) have continuous derivatives. 
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Figure 2.19 The curve a(t) = sin 2t, y(t) =sint with —4% <t< % 


The curve in Figure 2.19 can be specified in parametrised form by putting 
x(t) = sin 2t and y(t) =sint with —3* <¢< %. 

In general, if s denotes the distance along the curve from any fixed point, 
we have ds* ~ dx? + dy? as before, therefore 


(5) =) +8). 


and in the limit as dt — 0 we have 
ds\* _ (de)? | (dy 
dt) \dt dt} ~ 


ds dx\?  (dy\? : E 
Thus ae (¢) + (F) , 80 for the length L of the curve shown in 


Figure 2.19 we have 


Bar) “\2 2 
wan 
—3n/4 dt dt 
‘3n/4 
= Vv 4cos? 2t + cos? t dt ~ 6.917 84. The integral was evaluated 


80/4 by numerical methods. 
Such line integrals arise in a number of areas of applied mathematics. For 
example, it is relatively easy to calculate the work done when an object 
is pushed up a straight path against a constant force F, as shown in Fig- 
ure 2.20. We simply multiply the length of the path d by the component F’ 
of the force along the line of the path. 


Figure 2.20 A constant force acting along a straight path 


Imagine now that we push an object up a winding cliff path, and that we 
wish to find the total work done. In theory we could divide the path into a 
large number of almost linear sections, then calculate the work done when 
the object is pushed along the many short straight paths. and thus obtain 
an approximation to the total work done. We expect the approximation to 
improve as the number of subdivisions increases. 


In the following subsection we shall make this idea a little more precise. 
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Line integrals of vector fields 


Suppose that we are given a parametric curve C defined by a vector r(t), 
so the points on the path are determined by the values of the parameter t. 
Now suppose that a particle is propelled along the curve by a force F(r), 
and that we wish to calculate the work done by the force. If a constant force 
F acts on a particle and displaces it through a vector distance ds, then the 
work done by the force is F - 6s. So in order to calculate the total work 
done by the variable force F(r) as it moves along a curve, we approximate 
the path by a large number of straight-line subpaths on which the force and 
the direction are approximately constant (see Figure 2.21). We then have 
N 
total work done ~ aS F(r,) + orp. 
k=1 


nt) 


AY, 


1% + Ory 


Figure 2.21 A subpath ér, on the curve r(t) 


This approximation will improve as N increases (and the maximum length 
of the subpaths decreases), leading us to define a line integral as 


N 


= im » F(tn) + orn, 
= 


which we shall write in a more recognisable form shortly. 


Such integrals can also be used to represent quantities other than the work 
done. 


Suppose now that C is a curve in three dimensions, and that it can be 
parametrised by writing 


r(t) = x(t)i+y(t)j + 2()k. 


with the beginning and end points of C conrenpendling tot =a pea te b, 
pect = bol + byl + dk, 00 © = Py yy 8 
respectively. On C we have dr = dxi + dyj + dzk. so fe re i+ a + i k, 


where 6t is the increment in ¢ that gives rise to the increment ér, therefore 
dy dz 
ory (Gi+ eens *) ot = (2'(t)i+y'(t)j +2/(Ok) ot. 


If F(r) = F\(ax.y.z)i+ Fo(a. y, 2)j + Fa(a.y,z)k, so that on C the compo- 
nents of F are functions of t, then 


N 

SEF (rn) orn 

nes 

=S5 (Filtn)i + Fa(tn)i + Fa(tn)k) + (2"(tn)it y'(tn)5 + 2'(tn)k) 6th. 
n=1 


where ft; = a and ty = b, so Nét = b—a. 
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In two dimensions, 
r(t) = (x(t), y(t); 

in three dimensions, 
r(t) = (x(t), y(t), 2(t)). 


The i component of F is 
Fy (x(t). y(t). 2(t)) = F(t). 
and similarly for the other 
components. 
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In the limit as N — 00, 


J dx dy dz % dr 
r= f[ (ring mod +mG) a= | F.(G)a 


(which is just a ‘normal’ integral with respect to ¢). 


The line integral J is usually written as 


[Fo -dr, 


where it is understood that the curve C has been parametrised so that 
r = a(t)i+ y(t)j + 2(t)k is a function of a parameter t. 


The following example (in two dimensions) may help to clarify the notion 
of a line integral. 
Example 2.23 


Evaluate the line integral - F(r)+dr where F(r) = 27i + x7y?j and C is 


the part of a parabola defined by the equation y = 1 — x? from A, the point 
(0,1), to B, the point (1.0), as shown in Figure 2,22. 


Figure 2.22 Part of the parabola y = 1 


Solution 


First we parametrise the curve by writing x(t) =¢ and y(t) = 1—¢?, with 
0<t<1. Then t = 0 and t= 1 correspond to the points A and B, respec- 


d. 
tively, Also notice that a =1 and 4 = —2t. We then have 
a a 


5 + £71 — £7)74 


F(r) =27i+27yj= 
and 


dr 
dt 


= £ (e(oi-+ w(0s) =i-2tj. 


50 
1 
[Fo-a= [ (Pi + P(1— #7)2)) « (i — 29) dt 
Cc 0 


1 
= [ (P-280-P))dt=}. © 
40 


The evaluation of such integrals in three dimensions is very similar. 


Exercise 2.22 


(a) Calculate the length of the curve 2 = ¢?, y =¢* from t= 0 tot = 4. 
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(b) At time ¢ seconds, the position of a particle is given by x =t, y = t?/V2, 
t#/3. How far does the particle travel in the two seconds from ¢ = 0 to 
2? 


(c) Evaluate if F-dr when F(x,y,2) = x°yi+ 2j + ry*k and C is the straight 


G 
line joining the origin to the point (1, 2.3). 


2.3.5 The Gauss divergence theorem 


We can now move on to discuss the Gauss divergence theorem, which in- 
volves an integral over a volume and an integral over the corresponding 
surface. 


The flux of a vector field across a surface 


The concept of the flux of a vector field a surface arises in many 
areas of applied mathematics, but perhaps the easiest to visualise is the 
flow of a fluid or gas across a surface. Imagine the flow of a fluid in three 
dimensions, where at each point P with position vector r = (zr, y, z) the fluid 
has velocity F(r); i.e. the fluid velocity is described by the vector field F(r). 
Now imagine that we immerse a wire mesh in the fluid, and that we wish to 
estimate the yolume of fluid passing through the mesh in one second. The 
situation is depicted in Figure 2.23(a), where the mesh is shaped rather like 
the visor of a motoreycle helmet. 


Figure 2.23 A surface approximated by a ‘wire mesh’, and the volume of fluid 
passing through an aperture in one second 


First we consider just one small aperture in the mesh, and calculate the 
volume of fluid passing through it. Suppose that the aperture is at position P 
with area AA, and that fluid passing through this aperture has velocity F’. 
Also suppose that the vector 7 is a unit normal to the surface at P. 


A vector perpendicular to the tangent plane at a point P on a surface is known 
as the normal to the surface at P. At every point on a smooth surface 
there are two normals, one ‘into the surface’ and one ‘out of the surface’. If 
the former is denoted by 7, then the latter will be —m. For a closed surface 
(ie. a surface without a boundary) we usually take the outward normal to be 
positive. 


In general, the fluid velocity and the unit normal will not be parallel, and 
in Figure 2.23(b) the shaded volume represents the fluid passing through 
the area AA in one second. This is a parallelepiped with base area AA and 
length F (see Figure 2.24). 
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is not a real wire 
but rather a 
mathematical construct 
with infinitesimally thin 
and ultimately with 
an infinitely fine mesh. 


If the aperture is small 
enough, then the fluid 
velocity will not vary 
across it. 
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Figure 2.24 


From the diagram we see that the volume of this parallelepiped is equal to 
the volume of a box with base area AA and height F-m, where 7 is the 
unit normal at P. Thus it is the component of velocity perpendicular to the 
mesh that is relevant, and the volume of fluid passing through the aperture 
in one second is approximately F(r) +n AA. 


Tf we sum over all apertures, and form the limit of the sum as the size of the 
mesh shrinks to zero and the number of subdivisions increases, we obtain 
an integral 


kL F(r)-idA. 
This integral is the flux of the vector field F across the surface S, which in 
this case represents the volume of liquid crossing the surface in one second. 
It is often convenient to abbreviate A dA to dA, in which case the integral 
is written | F(r)+dA. 
s 


Example 2.24 
If the velocity of a fluid is given by the vector field 
Fa op rit a? y2?j + 272k, 


calculate the volume of fluid per second passing through the rectangular 
area S given by r= 1,0< y<1,0<2<1. This situation is depicted in 
Figure 2.25, 


Figure 2.25 


Solution 


The unit normal to S is given by 7 =i, so dA = idA =idydz. Hence the 
volume of fluid passing through S is 


v= [rans -dA 
Ss 
1 


1 
= if (f (y?27i + yo?j + y?2k)- idy) dz 
JO 0 
1 1 1 
=| (f say) a= f idz= 3. 
0 0 0 


We say that the flux of fluid across the surface S is 3- If the fluid velocity is 
in metres per second, and the coordinates along the x- and y-axes represent 
metres, then $ has area 1m?, and }m*s~! of fluid passes through S. 
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The vector dA = ndA has 
magnitude dA and 
direction normal to the 
surface. Since surfaces are 
in general curved, the 
direction of @ and dA will 
depend on r. 


Note that we could have 
chosen 7% = —i, in which 
case we would be 
calculating the volume of 
fluid passing through S in 
the —i direction, which 
would simply give a change 
of sign in the final result. 
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More generally, consider the flux of an arbitrary vector field F across the 
surface S of a small rectangular box with sides of length éx, dy and dz, as 
shown in Figure 2.26. 


ae 


(X,Y =) 


Figure 2.26 An elemental rectangular box 


The flux across the right-hand shaded area of this box is approximately 
F(x + 6x, y, z) « idy dz, whereas the flux across the left-hand shaded area is 
approximately —F'(a,y, 2) + idy éz (because here the outward normal is —i). 


Applying a similar argument to the other two pairs of parallel surfaces, we 
obtain an approximation for the flux across the entire surface S: 
if F-dA~ (F(x +6z,y,2) — F(x,y,2)) -idy dz 
a + (F(a,y + dy.z) — F(x. y.2)) + jér dz 
+ (F(x. y,2 + 62) — F(a,y,2)) + kda dy 


, (Fila + 6x, y, z) — Fi(a.y.2) 
= 6V 
bx 
Fo(x,y + by, z) — Fo(x,y.2) 
pe EE 
by 
F3(x,y,z + dz) — wang) 
‘errr |e 
bz 
where dV = da dy dz. If we now take the limit as the volume shrinks to zero, 
we obtain 
n 1 OF, , OFh | OF 
a . = — + — + —=V-F, 2 
a av fe Co eo eee 2s) 


where V is the vector differential operator discussed on page 85. This result 
has an interesting consequence when F describes the velocity of a fluid. 


‘The integral [ F -dA represents the volume of fluid crossing the surface in 
s 


one second (both into and out of the region). For an incompressible fluid, 
where the amount of fluid entering the region equals the amount exiting, we 
would expect this volume to be zero (unless there is a source of fluid inside 
the region or a sink allowing some fluid to escape). In other words, for F 
to represent the velocity of an incompressible fluid, we require V - F = 0 
throughout the region. 


The Gauss divergence theorem 


Equation (2.6) is valid for an infinitesimally small volume element dV. The 
Gauss divergence theorem generalises this result to a finite volume R. 
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Gauss divergence theorem 


For any (sufficiently smooth) vector field F, the volume integral of 
V -F over a region R is equal to the (outward) flux of F across the 
surface S$ enclosing R. This statement can be expressed as 


[v-rav = [ P-aa. 
JR JS 


Figure 2.27 Schematic representation of the Gauss divergence theorem 


This result is easily established from equation (2.6). First we use a Cartesian 
grid to divide the region R into a large number of box-shaped volumes 6V, 
with corresponding surfaces 6S;. Then 


[vera = D0 - PM, 
R i 


where (V+ F); means V - F evaluated at some point within the small vol- 
ume 6Vj. Using equation (2.6), we obtain 


v-Fav ~ >> ( F-dA.), 
Fi 6S, 


where F'-dA; represents the flux of F out of the box-shaped region 6Vj. 
5S, 


YR 


Figure 2.28 ‘Two adjacent elemental boxes 


Now, adjacent boxes have a rectangular face in common, and the outward 
normals on this common face have opposite signs (as shown in Figure 2.28). 
The flux of F out of this face from one box is equal in magnitude to the flux 
of F into the adjacent box, but these fluxes have opposite signs, so their 
contributions to the sum cancel. This cancellation occurs for all the faces 
except those at the boundary of the region R, so we have 


[vera [ P-aa. 
R s 


As the number of elements in the subdivision increases, and their size de- 
creases, this approximation becomes more accurate, and in the limit we 
obtain the required result. 


The Gauss divergence theorem is a very important result. It relates an 
integral over a region to an integral over the boundary of that region. In 
this sense it is a generalisation of the well-known result in one dimension 
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b 
dj 
# ae = f(b) — f(a). In fact, the theorem can be generalised to any 

number of dimensions. At a practical level, it enables us to calculate the 

flux of a vector field by calculating a volume integral of a scalar field. In 
general, this is much simpler than calculating integrals of vector fields over 
surfaces. 


Example 2.25 
Verify the Gauss divergence theorem for the vector field 
F = ayi+ y2j+ 2k 
and the region R given by the unit cube 0 <2 <1,0<y<10<2<1. 


Solution 
First we integrate V - F over the region R. Since 
Died) Atare a 

Stew), Olye) Slee) — 


es Ox Oy Oz 


ytete, 


we have 


[vere = i (/ ( [ori dr) a) aa 
‘s : (f G@ ++ :)av) dz 


1 
= [ G+de=$. 
Next we integrate the field F over the surface of R. Since F is a cube, 
we need to evaluate six integrals. one for each of the faces at x = 0, x = 1, 
y=0, y= 1,2=0, 2=1. 
For the face at z = 0, the unit normal out of the cube is given by n = —i, 
so the surface area element is dA = —idydz. Hence F(0.y,2)+dA = 0, and 
we obtain 
S=| F(0.y2)-dA=0. 

r=0 
For the face at « = 1, the unit normal out of the cube is given by 7 = i, so 
the surface area element is dA = idydz. Hence F(1,y,2)-dA = y? dydz, 
and we obtain 


1 1 
Sy = F(1.y.2)-da = [ (f yay) de = 
=i io \Jo 


For the face at y = 0, we have i = —j, so dA = —jdrdz. Hence we have 
F(x.0,2)+dA = 0, and we obtain 


S3= | F(a,0,2)-dA=0. 
y=0 


For the face at y = 1, we have dA = jdxdz, so F(x,1,z)-dA = zdxrdz, 
and we obtain 


Sim fFleas)-aa= [ ([ =«) z=}. 


For the face at z = 0, we have dA = —kdzdy, so F(x.y,0)-dA =0. and 
we obtain 


x= f F(x.y.0)-dA=0. 
==0 
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Here, the points a and b 
form the ‘boundary’ of the 
line interval [a,b]. 


96 Chapter 2 Multivariable and vector calculus 


For the face at z = 1, we have dA = kdrdy, so F(x,y,1)-dA = xdxdy, 
and we obtain 


1 1 
s='] F(x,y,1)-dA= (f zd) dy =}. 
2 0 SO 
Adding these surface integrals gives 


Si) + So + $3 + Si +55 + Se = 4. 


the same result as obtained by finding the volume integral. 1 


Exercise 2.23 


Verify the Gauss divergence theorem for the vector field F = x sin(xy)i+ xzj over 
the region 0<2<1,0<y<1,0<2<1. 


2.4 Changing coordinates 


In Subsection 2.2.8, we showed that if the Cartesian coordinates « and y are 
functions of two independent variables u and v, then any function w(x, y) 
can be differentiated with respect to u and v, using the chain rule. First, 
keeping v fixed, we have 


Ow Oy 
Oy Ou" 

Now keeping wu fixed, we have 
Ow _ OwOx | Ow dy 


=a + < 
Ov Ox Ov Oy Av 
In this section we use these results to transform between derivatives in polar 
and Cartesian coordinates. 


Polar to Cartesian coordinates 


Recall that in polar coordinates, points (a,y) in the plane are labelled by 
(r.0) where 2 = rcos@ and y = rsin@. 


Now suppose that we have a function w(x,y) and we wish to express its 
partial derivatives with respect to r and @ in terms of its partial derivatives 
with respect to x and y. 


First we find 


Ox Ox : Oy ‘ Oy 
=cos0, =-—rsin8, =sin@, =rcos6, 
ar cos 6. 30 rsin or sin 0. 20 rcos 8. 
and then, using the chain rule, 
Ow _ Owdx , Owdy _ du nu 


; v 
— cos 9 + sin@ 


oy 


Or dx dr" Oy dr Ox 


Pe & Ow y Ow 
~ Vata oe” Stee Oy 
and 
Ow 3 Ow Ox ne Ow Oy apa Ow fe reos@ ow = 7! dhe “file 
00 «Ox 00 Oy OO Ox Oy Oy Ox 
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Similarly, we can calculate the second-order partial derivatives with respect The algebra may be 
to r and 6. unpleasant, but the process 
is straightforward. 


Example 2.26 


For an arbitrary function w(«,y), calculate the second derivative with re- 
spect to @ in terms of the second derivatives with respect to x and y. 


Solution 


We have just shown that 
Ow Ow Ow 
= es =a 
since w is arbitrary, it is convenient to consider the operator 
) 0 oO 
30 Oy - YOe: 
This allows us to express a in terms of partial derivatives with respect 


to x and y in the form 
Pw dO (dw a) a Ow Ow 
ae 00 (3) > (5 -1x) (Fe - vz). 
It remains to simplify the expression on the right: 
Ow a a) Ow Ow 
ar (5 i Uys) (ox * rd) 
oO ( dw Oo ( dw a ( dw Oo ( dw 
(=F) ~y (vz) Yn (5) Yon (vz) 
é é Pw Ow Pw 2 w 
Viziy) ~¥ (35 **deay) *¥ act 
20?w Ow Ow 


Von "de oy ™ 


Example 2.27 

Given z = «? + y?, calculate oe ee and a in terms of r and @. 
Or’ 00 Orde 

Solution 

We have 


Oz  Oz0x , Oz Oy 
PE ECE: | SREY = Seen 5 
or x Or * Oy Or scapes 
= 2r(cos? @ + sin? @) = 2r, 
Oz Oz0r | Oz0y 
eS ee hee el ret ros 6 
30 a2 00 3 dy 00 mes rsin 0) + 2y(r cos 6) 
= 2r*(—cos@sin @ + sin @cos @) = 0, 

2 =p 

Orde 
In this example we could have noticed that z = x? + y? =r, from which 
the results follow immediately. 


Cartesian to polar coordinates 


Again we have the relationships x = rcos@ and y = rsin@, but while pre- 
viously we considered x and y to be functions of r and @, now we suppose 
that r and @ are functions of x and y. 
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. Pw . , ey . 
Suppose that we wish to express >— in terms of partial derivatives with 


respect to r and @. One might suspect that it would be useful to write 


r=Ve2+y? and @=tan'(y/x) for >0, 
and proceed as before. However, using this approach the algebra becomes 
unpleasant, and there is a better way. 


We begin with the equations = rcos@ and y = rsin@. First we differenti- 
ate them partially with respect to x, to obtain 


or Ne aoe or. - 
1= 5, cosd—rsind 2 and O= 7B sind + rcosd 


and these equations can be solved to give 


or 00 sind Note that en and kel are 
> = cos? and ~~ =———. Ox Ox 
Ox Ox . not simply the inverses of 
Similarly, differentiating the equations x =rcos@ and y = rsin# partially dx nd Ox 
with respect to y, and solving the resulting equations, we obtain Or 7° 06° 
OF. 00 cosé 
ay sin@ and i re 
These relationships are exactly what we need in order to find expressions 
for — and ae in terms of — and ne From the chain rule, we have 
of Oe Oy OF 00° a 
dw _ Ow r | Ow 20 _ Ow _ sind Ow 
Ox Or Ax * OO Ox or r 00 
and 
Ow _ Ow dr - Ow 00 _ sin Qo fa cos 0 Ow 
Oy Ordy Ody Or + OO” 
If we now replace w by a in the first of these expressions, we obtain 
Pw eae 9 (dw) _ sind dO (dw 
Ox? Or \ Ox r 00\ 0x )* 
Ow 
and using the above expression for 2’ this becomes 
Pw 30 a oad Ow _ sinddw\ _ sind 0 pow _ sind dw 
BaF FO” Or 80 Be" Orr 8B 


Finally, using the rule for differentiating a product, we have 


Pw _ ; ay ow _ 2sin@cos@ Pw e sin? 0 Pw 
ae 8 Or? r 0ra0” POP 
sin?@dw , 2sin@cos0 dw 
Tot 88 27) 


In the following example and exercise we assume that a function 2(r,@) be- 
comes a function 2(:,y) after a change from polar to Cartesian coordinates. 
Thus the independent variables change from r and @ to x and y, while the 
dependent variable remains as 2. 


Example 2.28 


a. 
Given z = rsin 20, show that r— + y—— ==. 


Ox Oy 
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Solution 

We have 
Oz Oz0r 0200. : ae sind 
Du Bat DO ae an eos + Or cos20 ( = ) 


= sin 20cos @ — 2cos 26sin@ 


Oz0r 0200. ‘ cos 4 
Or By + 20 dy > sin 26sin # + 2r cos 20 =a, 
= sin 20sin @ + 2 cos 20 cos, 


ot ae = rcos 6 (sin 20. cos 0 — 2. cos 20 sin @) 
+ rsin 0 (sin 20 sin @ + 2cos 20 cos@) 
=rsin20(cos*@+sin?0)= 2. 


Exercise 2.24 
Pw 


. » . ; w, 
Given a function w(2,y), find an expression for yz im terms of r and 0, where 
A Pw, Pw 
a =rcos@ and y=rsin@. Hence find an expression for ra + rd 
ss paneer ! ed 
partial derivatives with respect to r and @. y 


in terms of 


Laplace’s equation in polar coordinates 


The operator 

& Gees wi 

aa? * oF OF 
is known as the Laplacian operator; you will meet it again in Block I. 
The partial differential equation 


Pw Pw Pw _ 0 

On? Oy On 
is the three-dimensional Cartesian form of Laplace’s equation, which 
arises in many areas of applied mathematics. For various reasons, it is 
often useful to transform this equation into spherical or cylindrical polar 
coordinates, and to do that we require methods similar to those developed 
above. Applying the results of Exercise 2.24, the two-dimensional form of 

: . Ow Pw 

Laplace's equation a + => = 0 becomes 


Oy? 

Pw i: 1 
Or? * 1? ag? 
(This form of Laplace’s equation may appear to be more complicated than 


the original. In fact, as we shall see later when we discuss separation of 
variables in Block I, this form can be advantageous.) 


Exercise 2.25 


Transform the partial differential equation 


dw , Ow 
ax +b =0, 
dx Oy 
where a and b are constants, by making the change of coordinates u = br — ay and 
v = ax + by, so the new independent variables are u and v. 
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Exercise 2.26 


Transform the partial differential equation 

Pw Pw 

a? OF 
by making the change of coordinates u = x+y and v = x — y, so the new indepen- 
dent variables are u and v. 


End-of-section Exercises 


Exercise 2.27 


Find the cosine of the angle @ between the vectors a = 2i—j+3k and b=i+2j—k. 


Exercise 2.28 
Given a = i+ 2j+3k, b = 2i+2j —k and e=i—j+ 2k, calculate 
(a-c)b—(a-b)ec. 


Exercise 2.29 


(a) Show that the curve y? = daz, where a is a constant, can be defined in param- 
etrised form by writing 2 = at? and y = 2at. Find the length of the curve 
from the origin to the point (a, 2a). [Hint: Express the length as an integral 
in terms of f, then use integration by parts.] 


(b) Show that the surface 


can be specified in terms of two parameters u and v in the form « = 2(u+v), 
3(u—v), 2 =4uv, and show that the point P with coordinates x 
|, = 0 lies on the surface and corresponds to choosing u = 1 and v = 0. 
Find the equation of the tangent plane to the surface at P, and find a vector 
perpendicular to this plane. 


Exercise 2.30 
(a) Find the equation of the tangent plane to the surface 
z= f(x,y) = (e+ 2y + 1)(e—y +2) 


at the point (0,0,2) on the surface. Show that the point corresponding to a 
vector of the form 


r= 2k + Mi +3k) + u(j +3k), 


where \ and j are arbitrary scalars, lies on this tangent plane. 


= 


Given @(«.y, 2) = z— f(a, y), evaluate Vo at the point (0,0, 2). Show that any 
vector r — 2k (which is parallel to the tangent plane at (0,0, 2)) is perpendicular 
to V¢@ evaluated at the point (0,0.2). 


Exercise 2.31 


(a) Calculate the gradient of the scalar field F(x, y, 2) 
at the point (1,—1,2). 


Bry + y> + Qz+2 


(b) Calculate the divergence of the vector field F = sinh(xy)i + cosh(yz)j + e*"k. 


2.5 Outcomes 


Exercise 2.32 

Find the work done by the force F = (x + yz)i+ (y+a2)j+(2+zy)k in moving 
a particle from the origin to the point P with coordinates (1,1, 1): 

(a) along the straight line joining the origin to P; 


(b) along the curve r =t. y =f? i: 


Exercise 2.33 


(a) The function @ is defined by o(2. y, z) iz, and the vector field F is given 
by F = V¢. A point P with coordinates (cost, sint,sin 2t), 0 < t < 2z, is an 
arbitrary point on a curve C. Show that the point P; corresponding to t = 0 
coincides with the point P corresponding to t = 27, and show that 


[ F(r)-dr =0. 
ic 


[in Calculate S (cos? tsint sin 2t). 


(b) Generalise the result of part (a) by considering an arbitrary vector field F = 
V¢ for some function (x,y,z), where P, with coordinates (x(t). y(t), z(t), 
ty <t < te, is an arbitrary point on a curve C. If the point P; corresponding 
to t; coincides with the point P2 corresponding to tz, show that 


[Fo)-ar=o. 


2.5 Outcomes 


After studying this chapter you should be able to: 

e understand the meaning of continuity for functions of many variables; 
e understand the meaning of and be able to calculate partial derivatives; 
e determine the tangent plane to a surface at a point; 
° 


calculate first- and second-order Taylor approximations for functions of 
several variables; 
calculate and classify stationary points of functions of two variables; 


¢ apply the chain rule for functions of several variables; 

e integrate over two-dimensional planar regions in Cartesian and polar 
coordinates; 

e integrate over non-planar regions when given an appropriate coordinate 
system; 

e understand the concepts of scalar and vector fields; 

e calculate the gradient of a scalar field and the divergence of a vector 
field; 

e calculate the length of a curve specified either by a function or in para- 
metric form; 

e calculate line integrals of vector fields; 

e calculate the integral of a vector field over a surface to determine the 
flux; 

e verify the Gauss divergence theorem in elementary cases; 

e use the chain rule to find derivatives with respect to new coordinates; 

® appreciate the derivation of Laplace’s equation in polar coordinates. 


You may assume that the 
curve C is smooth. 
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Solutions to Exercises in Chapter 2 


Solution 2.1 


The equation y = 2 
as its axis of symmetry. The equation x = 


is recognisable as the equation of a parabola with the y-axis 
y* also represents a parabola, but with 


the z-axis as its axis of symmetry. The equation x = a+ y? represents a translation 
of the latter parabola along the z-axis. These observations allow us to sketch the 
following curves. 


Figure 2.30 The curves y/x? = 


Figure 2.29 The curves $— 24 + 


The level curves in the first diagram are typical of a ‘well-behaved’ function, while 
the fact that the curves meet at the origin in the second figure means that the 
function is not defined there. 


Solution 2.2 


ve 3 + e¥ cos(weY) and uh = 3r7y? + xe” cos(re"). 
Or Oy 

Solution 2.3 

We have 


1 
| said (ti=rcowa)! 
o @ 


1 
ji sin(ar) de 
Jo 


and differentiating partially with respect to a gives 


1 
ww (f sin(ax) de) set (a cosa) 
da \Jo da \a 


But, using Leibniz’ 


ae 
am (f sin(ar) dr) -[%  (cin(ax)) dx 
da \ Js , 


therefore 


1 
[vain 
h 


[- cos(ar) 


a 


rule, 


: 
[ ecuelaz) ite: 


0 


ney 2.4 


Or 
eee = ye™Y + 2x and Py = ae" + cosy. 
Oy 


Solutions to Exercises in Chapter 2 


(b) Differentiating partially with respect to x gives $a7'/? + pez =0, so 

a: & az SE Ox 

‘a > 0). Similarly, —~ = —*= (y > 0). 

fie ar ( ) eae (y>0) 

(ce) Ifcos(ayz) = 1, then xyz = 2k7 for some integer k, so z = 2k / (ay). It follows 
a e ca Oz _ 
that re —2kr/ (xy) and ime Qkr/(ay?). 
(d) (i) Putting « = 1 and y = 2 gives z = —1, thus the point (1,2,—1) does lie 
on the surface. The partial derivatives are 
Oz 
= =(y—22 —2(a — 2; =5y=40— 
Oe (y — 2x + 1) — 2(@ — 2y + 1) = 5y —4e —-1 
and 
Oz 
= = 5a —4y—1. 
Oy a 
Evaluating the partial derivatives when 2 =1 and y = 2, we obtain 
F,(1,2) = 5 and F,,(1,2) = —4. Thus the equation of the tangent plane is 
2=5(x—1)—4(y—2)-1. 

(ii) Putting 2 = 1 and y = 1 gives z = 1, so the point (1,1, 1) does lie on the 
surface. At this point, the partial derivatives are F,(1,1) = F,(1,1) = 0. 
from which we deduce that this is a stationary point. At (1,1), the 
equation of the tangent plane reduces to z = 1, so it is parallel to the 
(a, y)-plane. 

Putting s =a — 1 and t = y — 1, we obtain 
F(a,y) — F(1.1) = (s — 2t)(t — 2s) 
= —2s? + 5st — 2¢? 
= -2(a — 1)? + 5(x — 1)(y—1) — Ay—1)2. 
You may have noticed that this expression contains no linear terms in 
powers of z — 1 and y — 1; this is a consequence of the fact that there is 
a stationary point at (1,1). 

(iii) Putting y = 2-2 gives F(x,y) — F(1,1) = —9(a — 1)? < 0, and putting 

y = gives F(x, y) — F(1,1) = (@ — 1)? 20. 
We conclude that every neighbourhood of the point (1, 1) contains points 
at which F(a,y) — F(1,1) is positive and points at which it is negative. 
Thus the stationary point at (1,1) cannot be an extremum. 
(ce) We have 
= = (a? + pP)e* + 2xe* = ((x +1)? +y?—1)e* and > = eA 
so stationary points occur at the intersections of the circle (a + 1)? +y? =1 
with the line y = 0, namely at (0,0) and (—2,0). 
Since F(«r.y) — F(0,0) = (x? + y2)e* > 0, the stationary point at the origin is 
a local minimum. 
Solution 2.5 
(a) we = ye™ and w, = ze™, so wre = ye, wry = E™ + zye™ and wy, = 


ae", The quadratic approximati is therefore 
w(a,y) + w(1,2) + w,(1,2)(2 — 1) + wy (1,2)(y — 2) 
+ 3 (wee (1, 2)(@ — 1)? + 2wey(1,2)(— 1)(y — 2) 
+ wyy(1.2)(y — 2)") 
=e" (1+ 2(a —1) + (y—2) + 2(x — 1)? +. 3(x — 1)(y — 2) 


+4(y—2)?). 
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Shortly we shall see a 
simple criterion that will 
allow us to reach this 
conclusion more easily. 


You can find these by 

substituting y = 0 into the 
equation (2 + 1)? +4? =1 
to obtain 2 =0 or x = —2. 
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(b) Since 
the solutions of the pair of equations cos x cosh y = 0 and sin rsinh y = 0. Since 
cosh y > 1, the first of these equations gives cos = 0, so 


2n+1 
2 


Ow 
=cosrcoshy and aoe sinzsinh y, the stationary points occur at 


a= §+nn= ( )s for n an integer. 


For these values of x we know that |sin | = 1, so the second equation gives 


y=0. 
The second partial derivatives are 
a - ew oe 
a = -—sinwcoshy, a = cos rsinh y, oe =sinzxcoshy. 


Using the previous notation, at the point (5 +n7,0) we have 


PR—Q? = —(sinxcosh y)? — (cosrsinh y)? = 


if 


so the stationary points are not extrema (but saddle points). 
Solution 2.6 
oy 


(This simple observation will be of interest later, when we discuss the general 
solution of partial differential equations such as F,, = 0.) 


If F(x.) = f(z) +9(y), then & = g'(y), so aa =0. 


Solution 2.7 


Putting @ = 0 (i.e. x; = x2 = rg = 0), we obtain A = F(0). Then differentiating 
partially with respect to each of the variables in turn, and putting #2 = 0, we have 


OF OF OF 
= » b=(— > B=([—— 5 
Oxy Es ; a ie (ae, = 


In order to find the remaining coefficients, we differentiate partially twice and then 
put @ = 0, giving 


<A Poee ae fae 1 (PF 
cum 3 (3a). cu= (33), cu = (Fa) 


PF PP 
ys ae aomenias* a C= he 


Solution 2.8 


(a) We have w = e085", so 
ae = (cos? t — sin? t)esm tort, 
Using the form w = e*”, we have 
Ow dx dy 
= yer, = ae", = cost, =-—sint. 
Dr ye Oy re at cos t a sin 
Hence, using the chain rule, 
dw _ Owe | dwdy 
dt Ox dt " Oy dt 


= ye™ cost — re™ sint 
= (cos? t — sin? te" *et, 
which agrees with our first calculation. 


(b) We have 
Ox Ox Oy _ Oy _ 


ae d. 
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Then, using the chain rule, 


Ow = Ow Ox es Ow Oy re a (ou 

Ou Ox Au Oy Ou Ox Oy 
and 

dw - Ow dx | Ow dy _ ,Ow +qoe 

Ov Ox Av Oy Ov Or Oy” 


These relationships are valid for any function w, and it is convenient to write 
@) @) a a a a 
c: an = + * 


ay Eee iy eae fag a 
Then 

Aw _ (9 0) (dw , aw) _ ee w 9 Ow +2%e, 
fet regulary [lores ular) aaa or Oy * 
Pw a. A\(,dw , dw ‘i ee Pw 
aude = (0g. +5) (052 +05) = 5 art bebad) So + oda 
Pw (a a Ow | Ow a Pw - pw 
a (x +45) (0% +42) =f ra + bd oy +d DF 

Solution 2.9 


Figure 2.31 The region S$ divided into (a) vertical strips and (b) horizontal strips 
There are two alternative approaches. We could divide the region S into two regions, 
J and K, as shown in Figure 2.31(a), where: 


J is bounded by y = 0, 2 =4 and y= 2: 
K is bounded by x = 4, y = 0, x =8 and y = 16/2. 


We then have 


[eae [ears f dA, 
s J K 


where the areas J and K written as repeated integrals are 


ry ca 4 
[eae f 2(f ay) a= f adr = [t2']! = 64 
J 0 i) 0 


and 


Le wdA= fe (fe p) ae [10ede~ [8x7]; = 384. 


Thus 
[ a? dA = 64 + 384 = 448. 
s 
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Alternatively, we could divide the region S into two regions L and M, as shown in 
Figure 2.31(b), where: 


L is bounded by y = 2, y = 16/x and y = x; 
M is bounded by y = 0, x = 8, y = 2 and y= x. 


We then have 


if wdA= [ wdA+ v x? dA, 
s L M 


where the areas L and M written as repeated integrals are 


4 16/y am 16/y 
[@aa- f / ede ay = f Be] dy 
Jt J v 2 


* 4096 y' 
=| (Fr - 5) av = 108 
Bey ; 


= ) dy = 340. 
Thus 
? dA = 108 + 340 = 448, 
as above, 
Solution 2.10 
We have 
m/f 2 1 : m/2 
[Pada (/ Par) a= 3 do = x/6. 
s 0 0 0 
Solution 2.11 
We have 
Ow z Pw 2? — y? 
pirate and = —. 
Ox (a? 4- y?)3/2 Ox? (a? + y?)5/2 
Similarly (and taking advantage of the symmetry), 
Solution 2.12 
We have 
OF y-2? =a OF ____— ry 
de wee Oy eR 
Then 
OF PF _ 2y(3r?-y) 
Oydx Axdy (a? +4?) * 
Solution 2.13 
Zoe fe Ott eh 
(a) Ifw=eF cosy, then 5~ = e* cosy and By 7 sin. 
2 ay ” a Yo , 
s Pw Pw Fw Pw afi 


0 Far = eT eosy and Fa = e* cosy, which gives 55 + Fa 


Solutions to Exercises in Chapter 2 


2 9 Ow ow 
= 23 — 324? ee =-6zy. 
(b) If w = 29 —3ay’, then Oe 3a? — 3y? and Oy Gary. 


Pw aw ye, Ow, Pw _ 
So 5 = Gr and Ret ~62r, which gives 53 oF 
Solution 2.14 
We have 


Fy = —22+20r+2y and F, =—12+ 22+ 10y. 


The single stationary point occurs at the intersection of the lines 10x + y = 11 and 
a + 5y = 6, and thus is at (1,1). Since F,, = 20, Fy, =10 and F,, = 2, it follows 
that (using the standard notation) 


PR-Q? = 200-4 = 196 > 0, 


therefore the stationary point is a local minimum. 


Solution 2.15 


Let P, with coordinates (x,y,z), be an arbitrary point lying in the given plane. 
Then the distance from P to the origin is d= \/2? + y?+2*. The algebra is 
simpler if we minimise w = d? = «? + y? + 2? rather than d, and we note that on 
the plane we have y = 2a + 22 — 16, so w + (2x + 22-16)? +27. As P ranges 
over the plane, there will be a point ch is closest to the origin, and at this point 
the function w(sr, z) will have a mi m. At this stationary point, 

Ow 

on = 2a + 4(2a2 + 2z — 16) = 102 + 82 —-64=0 
and 


Ow 
Je = Ae + 2z — 16) +22 = Be + 102 — 64 =0, 


therefore « = z = 32/9, so y = —16/9. Since there is only one stationary point, we 
know that this must be the required minimum, and it occurs at (32/9. —16/9.32/9). 


We can confirm this conclusion by performing the evaluation 


Ox? Oz? Oxdz 


2 
) = 10x 10-8? = 36 >0, 


2 ae, Paeetsscd . F 
so the stationary point is an extremum, and since isa 10 > 0 this extremum is 
Hie hr 
a minimum, 


Solution 2.16 


We have 


Bw Ow Ow _ 3ary= 


Oxdydz ~ p2drdy OyOzdxr (+ y+ey? 


Solution 2.17 


We have 


<7 jeme " es 
if + psinOdp a= [ [dp?sin]°°*" ao = [cost osinoa, 
> \o 5 0 


Now we can either make the change of variable u = cos@, or note that sm @d@ = 
—d(cos 0): 


=f pcos® 
a ¢ psnede) dé =-} 
0 \Jo 


cos? 6 d(cos) = 2 [—cos* 6] 5 = 4. 
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Note that here we have 

chosen to eliminate y, but 
we could equally well have 
chosen to eliminate x or 2, 
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Solution 2.18 


Figure 2.32 


Dividing the region into horizontal strips, then summing over the strips (see Fig- 
ure 2.32(a)), we obtain 


Ife exp(x*) dA = i. (eso ae) dy, 


which cannot be evaluated directly in terms of the elementary functions. 


Alternatively, dividing the region into vertical strips, then summing over the strips 
(see Figure 2.32(b)), we obtain 


3 pela 3 i 
II exp(x?) dA = [ (/ ext) dx =/ [yexp(a?)]2" ax 
Ss ‘0 \Yo fi 
2 
= [3eexp(z2)] dz 
0 
3 
= 5 ler], 
= t(e?-1). 
‘The point to notice here is that the order in which we perform the integration may 
well affect the difficulty of the calculation. 
Solution 2.19 
(a) We have 


ja] =V1+4+1=V6 and |b} = Va+1+1= V6, 


and 
a-+b=(1,2,—1)+(2,1,1) =3. 
So 
eS 
cos = a =a: 


(b) The magnitude of the vector ¢ is V3? + 22 +5? = V38. 


Solution 2.20 
We have 
Fy = 2xsin{(x + y)r] + (x? + 2y?)rcos{(a + y)7] 
and 
Fy = 4ysin{(x + y)x] + (a? + 2y?)rcos|(x + y)z]. 
It follows that F,(1,1) = F,(1,1) = 37, therefore VF = 3n(i +). 
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Solution 2.21 
(a) We have 
_ OF, , OF , OFs_ ty? +2? 3 
Va a ty te @ryteye @eP tae 
2 
(b) We have 
Vjf= BE 4 OF hy = 2Qri + 2yj + 2zk, 


Oy 
therefore V + (Wf) =6. 


Solution 2.22 


(a) “ =2t and 7 ay = 3t?, so the length of the curve is given by 


[y( OE (#) tae [ Varvara 


-[ (4+ 9t?)'/? at 
0 

s9]* 
=$[a+9r)"] 
= £(37v37 - 1). 


(b) We have 2/(t) = 1, y'(t) = V2t, 2/(t) = #?, so the distance travelled is 
+(2), dt = is Vi+22 + tat 


-[ (1+#)dt 
0 


= [t+ 4e5]3 = 42. 


(c) The straight line joining the origin to the point (1,2,3) can be parametrised 
asc =t, y =2t, z = 3t from t = 0 tot = 1. Then r = (i+ 2j + 3k)t and 
1 
£ F(r)+dr= [ (2¢3i + 3tj + 4t°k) « (i+ 2j + 3k) dt 
fey 0 
1 
-| (6t + 14¢9) dt = 63. 
0 


Solution 2.23 


First we integrate V - F over the region R. Since 


_ Asxsin(zy)) , Awz) _ 
V:.F= he + a = sin(zy). 
we have 
1 1 1 
[v-Fw= | (f sin(ny) dr) ) as 
R 0 


Next we integrate the field F over the surface of R. Since R is a cube. we need to 
evaluate six integrals, one for each of the faces at x = 0,2 = 1, y= 0, y=1,z=0, 
z=1. 
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For the face at « = 0, we have dA = —idy dz, so F(0.y. 2) -dA = 0, and we obtain 


S,= F(0,y,z)-dA =0. 
2=0 
For the face at « = 1, we have dA = idydz, so F(1,y,2)-dA = sin(my) dydz, and 
we obtain 


1 1 
So= | F(l.y,2)+dA -[ (f sin(xy) dy) wa 
r= 0 0 * 


For the face at y = 0, we have dA = —jdadz, so F(x,0,2)-dA = —xzdirdz, and 
we obtain 


a= f {Flno:)-da=— ([ z:42) de =-4. 
kes 


For the face at y = 1, we have dA = jdxdz, so F(x,1,z)-dA = xzdardz, and we 
obtain 


1 1 
Si [ F(e,1.2)-da= [ (f sede) de =}. 
j=1 0 0 


There is no contribution to the surface integral from the faces z = 0 and z= 1, 
since in each case F- dA = 0. 


Adding the surface integrals gives 
2 
$1 + 52+ $3 +54 = 5 


the same result as obtained by finding the volume integral. 


Solution 2.24 


We have 


Fw _ 9 (dw 
Oy? — Ay \ Oy 


_,9f. ,dw , cosddw cos8 0 (. dw | cosddw 
= sind — | sind -— + —-—- ara sind — + ——— 


or or + OO 00 or 00 
os in?” 2sinAcos@ Ow ct cos? 0 Pw ns cos*OOw _ 2sin@cosd dw 
= ere + oreo + OF” + OF POO 


Hence, using equation (2.7), 
Ow i Pw Pw 1 Pw 4 Ld 
Ox? * By? — Ar? * 2 ag * r Or” 


Solution 2.25 

From the chain rule we have 
dw _ dw du Ow Ov _ pow Ps Ow 
Ox Oude” Ovdr Du” Ov 


and 

dw dwdu , dw dv Ow Ow 

By CUO Ou ey = Aa Oa 
therefore 

aoe +0 = (0% +05) +b (-«52 +050) = (a? + woe, 
The given partial differential equation therefore reduces to 

Ow 


maT 
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Solution 2.26 

We have 
du du av av 
te) yh ae By 


It follows from the chain rule that 
dw  Owdu , Owdv dw , dw 
Or dude * Wve du’ Ov 
and 
Ow — Ow du ve: Ow Ov _ Ow _ dw 
Oy Ou dy” Av dy Au Av™ 
These relations are valid for any function w, and it is convenient to write 


a9_ 29 0 

de Ou Ov 
and 

a_a a 

dy Ou Ov" 


in which case we have 


Ow _ ( 2) ee a _@w Pw , ew 
ae - 


da? ~ \Ou du * Ov) ~ du * Bud * Be? 
and 

Pw 0 d@\ (dw dw Pw ,Fw ew 

oF (x x) (F Z i) = Oe “aude * Oe 
It follows that 


Pw Pw Pw 


Ox? Oy? Oude” 
so the given partial differential equation reduces to 


Pw _ 
Oudv 
Solution 2.27 
We have 
cos 0 = a:b _ (2i—j+ 3k) -(i+2j-k) 
fall] 2? + (1 +8 1? +2? + (1)? 
= 252-8 
~ VuVv6 
3 
~~ 2var" 
Solution 2.28 
We have 


a-c=(i+2j+3k)-(i-j+2k)=5 
and 

a-b= (i+2j + 3k) - (21+ 2) -—k) =3, 
so 

(a+ c)b—(a- b)e = 5b — 3c = 5(24 + 2 — k) — 3(i- j + 2k) 
i+ 13j— 11k. 
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Solution 2.29 
(a) Since t = » owe have x = at? =a (ey = cm so y* =4azr. 
2a" 2a 4a 
Since x = at? and y = 2at, we have 
de dys 
a 2at and me 2a. 


Thus the length of the curve is 
1 1 
S= ii v (2at)? + (2a)? dt = 2a f V1+# dt. 
0 0 
1 
Putting J = t V¥1+?# dt and integrating by parts, we have 
0 
aca 1 jy 
T= |t¥1+#?| — dt 
Wel, f ae 


14+? 1 
=v2-— ars [ at 
te lw lp VIFE 


=v2-1+ [sinh-' ¢]?, 


$= 2al =a(v2+In(1 + v2). 


(b) Just as a curve can be parametrised by a single parameter as in part (a), so a 
surface can be parametrised by two parameters, and here we have « = 2(u + v), 
y =3(u—v) and z= duv. Solving the first two of these equations gives 
es Spd reyegee! 
w=ats and = 2u 273 


Then substituting into the third equation gives 


2 2 
x y 
eee eee 
4 9 
showing that the parametrisation in terms of u, v and w does specify the given 


surface. 


Putting u=1 and v =0, we obtain « = 2, y= 3 and z = 0, so the point 
(2,3.0) does lie on the surface. 


2 YP 
Differentiating = = — — ry partially gives 


4 
Oz = Oz Qy 
Crt and Be. So 


so z,(2,3) = 1 and z,(2,3) = -2. The equation of the tangent plane at the 
point (2,3, 0) is therefore 


2=(«—2)-2(y—3), 
which can be written as 3a —2y—3z = 0. 


There are several ways of finding a vector perpendicular to the tangent plane, 
but the following is probably the simplest. First we notice that the plane 
3a — 2y — 32 = 0 passes through the origin, and if we let (x,y. z) be any vector 
from the origin to a point in this plane, we can see that (3, —2,—3) - (x,y, 2) = 
3a — 2y — 32 = 0. So the vectors (x,y,z) and (3, —2.—3) are orthogonal, thus 
the vector 3i — 2j — 3k is perpendicular to the tangent plane at (2.3.0). 


Tt is also worth noting that the equation of the plane may be written as 
f(x,y, 2) = 3x — 2y — 3z = 0, and then the vector Vf = 3i— 2j — 3k is per- 
pendicular to this plane. This is because Vf is always perpendicular to the 
(tangent plane to the) level curves of f. 
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Solution 2.30 
(a) We have 


of =(x+2y+1)+(e@—y+2)=2e+y+3 
and 

of 

oy 
The equation of the tangent plane at (0,0. 2) is therefore 

z= f2(0,0)x + f,(0,0)y + f(0,0) = 3x + 3y + 2. 


The point corresponding to the vector r = 2k + A(i+ 3k) + u(j + 3k) has co- 
ordinates (A, j4,2 + 3A + 3), and these satisfy the equation 2 = 3a + 3y + 2, 
so the point does lie in the plane. 


(b) We have 


= 2(n —y+2)— (24+ 2y+1) =27—4y4+3. 


06; , 90, 0b, __ Of; OF 


Via ae 6, a= Bee ait & 
whose value at (0,0,2) is —3i—3j +k. 
Tt follows that 


Vo (r — 2k) = (—3i — 3 + k) - (Ai + pj + 3(A + 1k) = 0, 


so the vectors are perpendicular. 


In order to clarify what is happening here, note that the level curve (x, y,2) = 
z— f(x,y) =0 corresponds to the surface z = f(x,y). So at a given point on 
the surface, vectors which lie in the tangent plane to z = f(x,y) will always 
be perpendicular to Vo. 


Solution 2.31 
(a) We have 
F, =2r—3y+22, Fy=—3e+3y?, F,=2r+1, 


VF = (2x — 3y + 22)i + (—3x + 8y?)j + (22 + 1)k. 
At the point (1,—1,2), we obtain VF = 9i+ 3k. 


(b) We have 
_ Asinh(xy)) , A(cosh(yz)) | Ale**) 
V-F= Or + Oy + Oz 
= ycosh(ry) + 2sinh(yz) + xe**. 
Solution 2.32 


(a) We can parametrise the straight line joining the origin to the point P by writing 
x=t,y=tandz=t withO<f<1. Then 
F(r) = («+ yz)i+(y + x2)j+ (z+ 2ry)k 
=(t+@)i+ (t+ )j+ (t+P)k 
=(t+#)(i+j+k) 


2 AU) Oe an 
fat at geaititk 


so the work done is 
a 
[Poa [ereirire-G+i+nde 
fo} 0 


1 
=3/ (@+@)at=3 [bP + 18)) = 3. 
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b) The particle moves along the curve x = t, y = #?, z = f° from the origin to the 
is rigin 
point P as ¢ increases from 0 to 1. On this curve, we have 


F(r) = («+ yz)i+ (y+az2)j + (2+ 2y)k = (t+ O)i+ (PF + tj + 26k 


and 


Gr de, dy. dz, oa 
mt ait ait gen it ts k, 


so the work done is 


1 
i F(r)-dr= J (+ O)i+ (2 + tj + 25k) - (1+ 2tj + 307k) de 
IC 0 


1 
= [ (t+0° + 20(t? + ¢4) + 6¢°) dt 
0 


1 
=| (t+2t° +90) dt =§. 
0 
Generally, the value of a line integral depends upon the path of integration; this is 
an example of one that does not. For the general vector field 


F(r) = F\(#,y.2)i + Fala, y. 2)j + Fa(a.y,2)k, 


it can be shown that the integral vi F(r) dr is independent of the path C between 
two given end points if there exists a function (x,y, 2) such that 

Oo a@ =, _ a. 

3a’ i= Dy’ Fy = 55% 


in other words, if F(r) = V@. In this example it is easy to verify that we can use 
= 4 (a? +? + 27) + aye. 


F,= 


Solution 2.33 
(a) We have 
(cos 0, sin 0, sin 0) = (1,0,0) = (cos 2x, sin 27, sin 477), 
so P, and Py are the same point. 
Now 
F=Vo= rei eho +2x= Qryzi + x23 + 27 yk, 


so 


[ F(r)-dr= vis (Qaryzi + x? =j + x*yk) + (dari + dyj + dzk) 
Fa 


" (an(t)y(t)z(t)i + x(t)? 2(t)j + a(t)? y(t)k) + Pe els th dt 
tata 
pe 
= (2costsin t sin 2ti + cos” tsin 2tj + cos*tsintk) + (—sint i+ costj +2cos2tk) dt 
0 
ae . : 
= (—2cost sin? t sin 2t + cos* t sin 2t + 2cos* t sin t cos 2t) dt 
0 


Qn gy 
=[ 77 (cos? t sin tsin 2t) dt 
ot 


= [cos” tsint sin 2¢]." =o 


(b) Using the chain rule, we have 


= 00, 00, _ Ao dx, dy, , dz 
[Fea ff ( Fei + O85 4 x). (Fit qe x x) a 
= Oodx | dody | Oodz 
-{ Ga ans cs 


ta 
| (4) a= [o()]2 =0. since 9(t1) = o(t2). 
th 
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